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1.   INTRODUCTION 

The  Saval  Ranch  Integrated  Model  was  developed  over  the 
period  from  November,  1981  to  January,  1982 .   Two  workshops  were 
held,  with  a  total  of  forty  participants,  who  provided  the  insight 
and  data  necessary  to  construct  and  refine  the  model. 

All  of  the  functional  relationships  and  assumptions  of  the 
model  are  described  in  a  detailed  report  (Sonntag,  et  al.,  1981), 
which  should  be  consulted  when  reading  the  following  documentation.* 

The  model  is  designed  for  the  MTS  operating  system  and  the 
SIMCON  simulation  language.   MTS  is  available  on  computers  at 
several  universities;  SIMCON  may  be  purchased  from  the  University 
of  British  Columbia,  Institute  of  Animal  Resource  Ecology.   The 
model  is  currently  stored  on  tape  at  Wayne  State  University.   Persons 
or  institutions  interested  in  using  the  complete  model  would  need  to 
obtain  permission  from  the  Bureau  of  Land  Management,  and  liaise 
with  ESSA  Ltd.  regarding  transfer  of  the  model  from  Wayne  State. 
For  those  persons  interested  in  modelling  particular  components  of 
the  system  (e.g.,  hydrology,  sage  grouse  population  dynamics),  the 
program  listings  presented  herein  contain  some  useful  techniques 
and  structures. 


*All  references  to  figures  and  tables  in  the  text  of  Section  I 
refer  to  Sonntag,  et  al . ,  1982. 
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2.   VARIABLE  DICTIONARY 

The  following  pages  describe  all  constants  and  variables 
used  in  the  integrated  model.   For  arrays,  these  indices  are  used 
to  explain  frequently  used  array  dimensions: 


ctype 


cattle  type: 

1  =  cows 

2  =  calves 

3  =  yearlings 


ncover  -   cover  types.   Two  breakdowns  are  used: 

(1)  grass,  shrub,  tree,  canopy  and  ground 

(2)  shrub,  bare,  grass,  rock  and  litter 


ngrt 


grouped  range  sites: 

1  =  claypan 

2  =  loamy 

3  =  riparian 


nlyr 


soil  layers: 
1=0-  10" 

2  =  10  -  20" 

3  =  20  -  40" 


nman 


type  of  vegetation  manipulation: 

1  =  ploughing 

2  =  burning 

3  =  spraying 

4  =  mining 


nmead 


1 


meadows: 

1  =  Tremenan 

2  =  Gance 

3  =  Haystack 
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npas    -   pastures: 

1  =  West  Darling 

2  =  East  Darling 

3  =  Lower  Mahala 

4  =  Middle  Mahala 

5  =  Upper  Mahala 

6  =  Lower  Sheep  Creek 

7  =  Upper  Sheep  Creek 

8  =  East  Independence  North 

9  =  East  Independence  South 

nper    -  deer  migration  phases; 

1  =  March  15  -  May  15 

2  =  May  15  -  June  15 

3  =  June  15  -  November  1 

nrt     -   range  sites: 

1  =  loamy  bottom 

2  =  wet  meadow 

3  =  claypan  10  -  12" 

4  =  loamy  8  -  10" 

5  =  claypan  12  -  16" 

6  =  loamy  10  -  12" 

7  =  loamy  slope  10  -  14" 

8  =  loamy  slope  14  -  18" 

9  =  aspen  woodland 
10  =  slope  12  -  14" 

nst     -   identified  sage  grouse  strutting  grounds  (see  Figure  6.5, 
page  104  of  Sonntag,  et  al . ,  1982) 

nveg    -   vegetation  types: 

1  =  bitterbrush 

2  =  sagebrush 

3  =  other  brush 

4  =  forbs 

5  =  perennial  grasses  (decreasers) 

6  =  perennial  grasses  (increasers) 

7  =  cheatgrass 
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stream   -   creek  3: 

1  =  Gance 

2  =  Mahala 

3  =  North  Fork 

weeks    -   the  period  from  March  15  to  November  30  corresponds  to 
weeks  1-38.   December  1  to  March  14  is  "week"  39. 
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2,       VARIABLE  DICTIONARY 


2.   1  Vegetation  Submodel 


VARIABLE 
ACNEW(npas,nrt,nlyr) 

AMEAD(3) 
ARET 

AWC (nrt,npas) 

CVPX(2) 
CVPY(2) 

CWST 

DPB 

DROOT 

ER(nrt) 

GMAX(nveg) 

GMRED (npas , nrt , nveg ) 

GPCWT (npas, nveg) 

GPLANT (npas, nrt, nveg) 

GS0XM(1,5) 
GSOXM(2,5) 
GSOXM(3,5) 
GS0Y(5) 
GTEMPX ( 3 ) 

GTEMPY { 3 ) 


DEFINITION 

available  water  capacity  by  soil  layer  by 
range  site  by  pasture 

area  of  hay  meadows  (a„) 

area  returned  to  pasture  from  seeding 
reserve  (a.) 

old  available  water  capacities  by  range  site 
by  soil  layer 

week  (for  sip  function  CVPY) 

cheatgrass  %  protein  «■  f  (week) 

cow  wastage  as  fraction  of  amount 
consumed 

net  growth  rate  in  above  ground  biomass 
(including  grazing) 

growth  in  root  biomass  (lbs/a.) 

erosion  index 

maximum  plant  growth  rate/week 

fractional  growth  reduction  due  to  soil 
moisture  (0-1) 

%  of  each  plant  type  on  each  pasture  (by 
weight) 

weekly  growth  rate  in  above  ground  biomass 
pregrazing 

claypan  %  moisture  (SLP  function  for  GSOY) 

loamy  %  moisture  (SLP  function  for  GSOY) 

riparian  %  moisture  (SLP  function  for  GSOY) 

growth  rate  reduction  =  F  (soil  moisture) 

temperature  (degrees  F.)  (SLP  function  for 
GTEMPY) 

growth  rate  reduction  =  F  (temperature) 
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VARIABLE 
GTPROD (npas) 

HAYOP(nmead) 

HAYP(nmead) 

HAYY{3,5) 

HMMR 
HX(5) 

HY(5) 

IFDATE 
I  LUMP 

IOC (npas) 

IVM 

NEWRS(ngrt,5) 

NPAS 

NRT 

PAST(npast) 

PC(npas,nrt,nveg) 

PCMAX (ngrt,ncoverl) 

PCMAX2 (ngrt,ncover2) 

PCWT (npas,nrt,nvag) 

PLNTMX(nveg) 


DEFINITION 

total  production  of  forbs  and  perennial  grasses 
by  pasture  at  flowering  time 

level  of  hay  operations;  1  =  poor  irrigation, 
2  =  good  irrigation,  3  ■  fertilization 

hay  production  by  meadow  (1  =  Tremenan ,  2  = 
Gance,  3  =  Haystack) 

SLP  function  for  hay  production  =  F  (level 
of  hay  operations  (HAYOP,  1st  subscript) 
and  streamflow  (HX,  2nd  subscript) ) 

%  moisture  causing  100%  mortality  of 
shrubs 

SLP  function  for  streamflow  (a.  -  ft. /a. 
of  meadow) 

temporary  copy  of  one  row  of  HAYY 

week  of  flowering 

switch  to  lump  old  range  sites  into  new 
ones 

array  used  to  skip  over  pastures  (skip  if 
IOC (npas)  =  0) 

switch  for  vegetation  manipulation;  none 
if  =  0 

array  for  grouping  of  old  range  sites  listed 
in  columns  into  new  range  sites  (rows) 

#  pastures 

#  range  site  types 
pasture  areas  (a.) 

%  cover  of  each  vegetation  type 

%  of  cover  types  (1)  for  each  grouped  range 
site 

%  of  cover  types  (2)  for  each  grouped  range 
site 

%  of  each  plant  type  by  range  site  by  pasture 
(by  weight) 

maximum  vegetation  biomass  (lbs/a.) 
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VARIABLE 
PMR 
PMX ( 2 ) 

PMY(2) 

PREM 

PRTS 

PSTR 

PTOT 

PVPX(2) 
PVPY ( 2 ) 

PWS (nveg) 
PWTMAX(ngrt,2) 

RA(npas?nrt,nveg) 
REQGA 

REQGB 

RLP 

ROOT (npas , nrt , nveg) 

ROOTMX(nveg) 

RPB 

RWS 

SHKILL(nman) 


DEFINITION 

weekly  mortality  rate  of  vegetation 

SLP  function  with  PMY;  vegetation 
weekly  growth  rate 

SLP  function  %  mortality  of 
vegetation  =  F  (growth  rate) 

fraction  of  range  site  removed  by  vegetation 
manipulations 

fraction  of  root  biomass  sent  to  shoots/week 
(before  flowering) 

fraction  of  shoot  biomass  sent  to  roots  (after 
flowering) 

total  weight  of  all  plants  on  one  range  site 
in  one  pasture 

SLP  function  with  PVPY;  week 

SLP  function;  %  protein  of  grasses  and  forbs 
(=  F  (week) ) 

plant  winter  survival 

maximum  %  by  weight  (1)  and  %  cover  (2)  of 
shrubs  and  grasses 

biomass  (lbs/a  .)  of  each  plant  type 

weekly  growth  rate  required  before  flowering 
(if  less,  get  from  roots) 

weekly  growth  rate  required  after  flowering 
(if  more,  pass  to  roots) 

ratio  of  %  cover  of  grass  and  litter  to  % 
cover  of  grass 

root  biomass  (for  all  plant  types  but 
cheatgrass)  by  pasture  and  range  site 

maximum  root  storage  (lbs/a.) 

(roots/above  ground)  at  start 

survival  of  root  potential  biomass  over 
winter 

fraction  of  shrubs  killed  by  type  of 
manipulation 
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VARIABLE 

SPB 

SPBEXP 

SVPX(all) 

SVPY(all) 

THAYP 

TOTPCT 

TPROD (npas,nveg) 

TRED 

VMANIP (npas,nrt,nman) 

VP (nveg) 
VSEED (nveg) 
WMAX 

WPA(l-3) 


DEFINITION 

exponent  for  setting  initial  plant  biomass 

exponent  for  setting  initial  plant  biomass 

SLP  function  with  SVPY;  week 

%  protein  of  shrubs  (=  F  (week) ) 

total  hay  production 

%  canopy  cover  plus  %  ground  cover 

total  plant  production 

fractional  growth  reduction  due  to  temperature 

#  acres  manipulated  in  given  pasture  and 
range  site 

%  protein  content  of  each  plant  type 

plant  biomass  after  seeding 

soil  moisture  (%)  in  wettest  soil  layer 
available  to  plant 

water/acre  for  3  hay  meadows 
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2.2   Hydrology  Submodel 


VARIABLE 

APAR 

ARGNC(nrt) 

ARMAH(nrt) 

AWCNEW (npas , nrt , nlyr ) 

BETOT 

BPAR 

BSTAB 

BSTB1 

BSTB2 

CATFAC 

CFAC 

CN 

CNA(ngrt) 

CNB (ngrt) 

CUFS 

DMTEMP 

DPRECT(npas) 

ER(nrt) 

ETPOT 

ETPLNT 

ETSOIL 

FLFUDG 

FLGNC 

FLMAH 


DEFINITION 

used  in  calculating  curve  number 

areas  of  range  types  in  Gance  basin 

areas  of  range  types  in  Mahala  basin 

available  water  capacity  in  inches 

total  bank  erosion  (tons) 

used  in  calculating  curve  number 

%  stable  streambank 

parameter  for  calculating  bank  stability 

parameter  for  calculating  bank  stability 

coefficient  relating  cattle  to  bank  erosion 

function  of  ground  cover  used  in  erosion 

SCS  curve  number 

parameters  for  calculating  curve  number 

parameters  for  calculating  curve  number 

converts  inch,  acres/day  to  cubic  feet/sec. 

daily  mean  temperature  (°F) 

snowmelt  and  precipitation  in  inches 

erodability 

potential  evapotranspiration  in  inches 

evapotranspiration  by  plants 

evaporation  from  soil  in  inches 

■ — not  used — • 

flow  in  Gance  Creek  in  cubic  feet/sec. 

flow  in  Mahala  Creek  in  cubic  feet/sec. 
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VARIABLE 
FREZDP (npas,nrt) 
FTOT (stream) 
FWTOT (stream) 
LUPTR(npas) 

NDAYS 

NFFLOW (weeks) 

OCN 

ORUN 

OS 

PACAN 

PMECH 

PMULT 

PREC(365) 

PRECIP(npas) 

RADLAN 

RFAC  • 

RFUDG 

RUNOFF (npas , nrt) 

S 

SEDDEL 

SMAX 

SMELT 

SNOPAC (npas, nrt) 

SO (npas, nrt, 3) 

SOFDG 

SOILOF (npas, nrt) 


DEFINITION 

depth  of  frozen  soil  in  inches 

total  flow  to  date  in  current  year  by  stream 

weekly  total  flow  in  each  stream 

pointer  from  pasture  to  hydrologic  soil 
group 

#  of  days  in  current  time  step 

flow  in  North  Fork  by  weeks 

curve  number  for  area  outside  Saval 

runoff  from  area  outside  Saval 

retention  parameter  for  area  outside  Saval 

percent  active  canopy 

factor  related  to  soil  conservation 

multiplier  for  rain  on  USFS  lands 

daily  rainfall  on  BLM  lands  (inches) 

--not  used — 

solar  radiation  in  langleys 

parameter  related  to  storm  intensity 

— not  used — 

water  runoff  in  each  spatial  unit 

soil  retention  parameter 

sediment  delivery  coefficient 

maximum  soil  retention  parameter 

snowmelt  in  inches 

snow  storage  in  inches  of  water 

weekly  mean  soil  water  as  %  AWC 

— not  used — 

erosion  in  tons/acre 

0  n  n  t  ? 


VARIABLE 


DEFINITION 


SOWAT (npas,nrt,3) 

TMAXL(365) 

TMAX(npas) 

TMINL(36  5) 

TMULT 

UEGNC 

UEMAH 

WATRES (nlyr) 

WAT  IN 

WATQAL 

WINDAY 

WQUAL 

WTWOOD 

XLS(nrt) 


inches  of  water  in  each  horizon 

daily  maximum  temperature  on  BLM  lands  (°F) 

weekly  mean  maximum  temperature 

daily  minimum  temperature  on  BLM  lands  (°F) 

multiplier  for  temperature  on  USFS  lands 

upland  erosion  into  Gance  (tons) 

upland  erosion  into  Mahala  (tons) 

depth  of  soil  horizons  in  inches 

infiltration  in  inches 

daily  water  quality 

#  of  days  in  winter  time  step 

weekly  mean  water  quality 

weighted  %  woody  cover  in  riparian 

length  slope  parameter  used  in  erosion 


,'\  f\  O 


J 


0 


2.3   Cattle  and  Economics  Submodel 


VARIABLE 

BEAT (ctype) 

CALF 

CALFS 

CALFW 

CE  (npas,ngrt,nveg) 

CLFB 

CMIN(2) 

COSTS 

COW (ctype) 

COWB 

COWMAX 

COWS 

COWW (ctype) 

CVA 

DAW 

DBH 

DGH 

DPW 

DPWN 

DSH 

FC 

GVO 

HAYB 

HAYR( ctype) 


DEFINITION 

weekly  biomass  intake  of  each  cattle  type 

#  calves  born  before  1st  year  simulation 

#  calves  sold  in  fall 

calf  weight  on  April  15  (lbs) 
initial  consumption  (lbs/acre) 

#  calves  bought  in  spring 

target  spring  weights:  cows,  calves 
total  ranch  costs  ($) 
initial  #  cattle 

#  cows  bought  in  spring 
target  #  cows  on  Saval  Ranch 

#  cows  sold  in  fall 
initial  cattle  weights  (lbs) 
county  value  added  (Equation  5.10) 
$/acre  to  improve  irrigation 
$/ton  to  buy  hay 

$/ton  to  grow  hay 

#  days  per  week 

#  days  in  winter  period 
$/ton  revenue  for  sold  hay 
fixed  ranch  costs 

gross  value  of  county  output  (Equation  5.8) 
tons  hay  bought 

lbs  hay  required  for  each  cattle  type  over 
winter 

nnn  1  a 
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VARIABLE 


DEFINITION 


HAYS 
HPP  ( 3 ) 

IAHB 

IAHS 

IENR 

IGR(npas) 

IGRT 

IPAST (npas) 

IPCNT(ipas) 

IPREF(nveg) 
IPTO (npas) 
IPW(npas) 
IPYN(npas) 
ISTOP 

ISTR 
IUSE 
IWM 

NPMAX 

PAS 

PCALF 

PCS 

PPCLF 

PPCOW 

PPYR 

RA(npas,ngrt,9 


-10) 


tons  hay  sold 

slope  of  winter  growth  versus  intake 
(Figure  5.8) 

hay  buying  switch;  =  1,  needed  hay  bought 

hay  selling  switch;  =  1,  excess  hay  sold 

week  cattle  taken  off  range 

slot  of  pasture  before  simulation  begins 

initial  #  weeks  to  date  on  pasture  NPAST 

— not  used — 

#  years  pasture  has  been  in  present  rotation 
slot 

plant  type  preferences 

transition  array  for  rotation  (Table  5„2a) 

#  weeks  cattle  graze  each  pasture 

#  years  pasture  is  kept  in  same  rotation  slot 

switch  for  computing  fall  and  winter;  =  0 
fall/winter  not  done 

week  that  cattle  put  on  range 

pasture  currently  being  used 

if  =  9/  feed  only  "near"  water,  =  10,  feed 
"distant"  too 

#  pastures  grazed  by  cattle 
proportion  cows  sold  in  fall 
%  calving 

proportion  calves  sold  in  fall 

$/lb  of  calf  bought/sold 

$/lb  of  cow  bought/sold 

$/lb  of  yearling  bought/sold 

proportion  of  range  sites  near  water,  distant 
f\f)f\  i  rfrom  water  (respectively) 
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VARIABLE 


REV 

RNA 

RNP 

SLI 

SLTDN (nveg) 

TGWA 

TGWAO 

THAYG 

THAYR 

TREQ  (C- 

type) 

TTREQ 

VCC 

WGAIN ( 

ctype) 

WGAR ( 8 

,6) 

XCLF(6) 

XH1: 

XM2 

XM3 

XWNG(all) 

YCLF(all) 

YEARS 

YEARW 

YII 


DEFINITION 

total  ranch  revenues  ($) 

ranch  value  added  (Equation  5.9) 

ranch  purchases  (Equation  5.7) 

slope  of  intake  versus  weight  relationship 
(Figure  5.3) 

slope  of  proportion  TDN  versus  proportion 
protein  (Figure  5.5) 

counter  for  acres  added  to  "near"  water  class 

counter  for  acres  added  to  "near"  water  class 

total  lbs  hay  grown 

total  hay  required  over  winter  period 

total  weekly  food  requirements  for  each  cattle 
type  (lbs) 

total  weekly  food  requirements  for  entire  herd 

per  cattle  variable  costs 

weekly  weight  change 

input  matrix  for  Figure  5.6.   The  model 
interpolates  down  column  1  (cattle  weight) , 
then  interpolates  across  to  get  daily  intake 
(lbs  TDN) 

X-array  for  calving  function  (Figure  5.7) 

economic  multiplier  (Equation  5.8) 

economic  multiplier  (Equation  5.10) 

economic  multiplier  (Equation  5.9) 

X-array  for  winter  weight  gain  (Figure  5.8) 

Y-array  for  calving  function  (Figure  5.7) 

#  yearlings  sold  in  fall 

weight 'yearlings  bought 

Y  intercept  of  intake  versus  weight 
relationship  (Figure  5.3) 


VARIABLE  DEFINITION 

YRB  #  yearlings  bought  in  spring 

YTDN(nveg)  Y  intercept  of  proportion  TDN  versus 

proportion  protein  (Figure  5.5) 

YWG(5)  input  Y-array  for  weight  gain  (Figure  5.6) 

YWNG(all)  Y-array  for  winter  weight  gain  (Figure  5.8) 
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2.4      Wildlife  Submodel 
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VARIABLE 
AHEN(nst) 

AMALE(nst) 

ANEED 
ANESMN 

ANX ( 3 ) 
ANY  (3) 

CAL(nveg, 12) 

CHENMN 
CHX ( 3 ) 
CHY ( 3 ) 

CMALE 

DAYH 
DENERF 

DENERA 

DFEM 
DINGW(2) 

DMALE 

FAWN 

FK 


DEFINITION 

initial  condition:;  adult  hens  on  strutting 
grounds  for  5  identified  populations 

initial  condition;  adult  males  on  strutting 
grounds  for  5  identified  populations 

kcal  *  10  *  3  needed/adult/week 

minimum  adult  nesting  success  if  forbs 
density  too  low 

forb  and  grass  biomass  (X-axis  of  Figure  6.6) 

adult  nesting  success  as  function  of  average 
density  of  forbs  and  grasses  (Y-axis  of 
Figure  6.6) 

nutritional  value  for  each  plant  type  by 
month  (103  *  kcal/lb) 

minimum  chicks/nest 

chicks/hen  (X-axis  of  Figure  6.7) 

number  of  chicks  per  hen  as  function 

of  average  forbs  and  grasses  density  (Y-axis 

of  Figure  6.7) 

proportion  of  sage  grouse  that  are  male  (in 
fall) 

#  days/hunter/year 

initial  condition  -  kcal  103  ingested  by 
fawn  deer  in  previous  year 

initial  condition  -  kcal  103  ingested  by 
adult  deer  in  previous  year 

initial  number  of  female  adult  deer 

lbs/week  ingested  by  1)  fawn;  2)  adult  deer 

initial  number  of  male  adult  deer 

initial  number  of  fawns 

half  saturation  parameter  in  extra  mortality 
function  applied  to  fawn  deer 

000  1  n 


VARIABLE 

FLEAVE(nveg) 

FNEED 

FPERC 

FRAC 

HFRAC 

HHX ( 3 ) 

HHY(3) 

ISTRUT(nst) 

NPASD  (riper, npas) 

NPERP(3) 

NPERT ( 3 ) 

NTYPE(nper,nrt) 

PERIN(all,all,all) 

PERMIT 
PREF(nveg,12) 

RRA ( npas , ngrt , nveg ) 

SAGEG 
SHUNT 
SK 


DEFINITION 

amount  of  forage  not  eaten 

kcal  *  10  *  3  needed/fawn/week 

proportion  of  fawns  that  are  male 

fraction  of  deer  population  overwintering 
east  of  Saval  project 

proportion  of  game  management  unit  covered 
by  Saval  project  area 

density  of  bucks  on  range  for  SLP  function 
(X-axis  of  Figure  6.4) 

deer  hunting  success  per  unit  of  effort  - 
hunter  days  (Y-axis  of  Figure  6»4) 

pastures  containing  sage  grouse  strutting 
grounds 

pastures  used  by  deer  in  each  of  3  periods 

number  of  pastures  utilized  by  deer  in  3  periocj 

number  of  range  sites  used  by  deer  in  3  period 

range  sites  used  by  deer  in  each  period 

initial  condition  -  proportion  of  each 
plant  type  in  deer  diet  by  pasture  and  range 
type 

#  deer  hunting  permits  in  game  management  unit 

monthly  preference  parameters  for  each  type 
of  vegetation 

dummy  variable  to  reflect  density  of  plants 
after  removing  cattle  grazing  -  otherwise 
same  as  RA 

total  number  of  sage  grouse  on  Saval  ranch 
March  15 

proportion  of  sage  grouse  population  removed 
by  hunting 

half  saturation  parameter  in  extra  mortality 
function  applied  to  sage  grouse 


r\  o  n  1 


n 


VARIABLE 

TOTHEN 

TOTMAL 

WE (npas ,ngrt ,nveg) 

WHENS 

WMALES 

XFAW(all) 

XFEM(all) 

XFNY(all) 

XHFAC(4) 

XHSUR ( 3 ) 

XMALED(all) 

XMSUR(all) 

YFAW(all) 

YFEM(all) 

YFNY(all) 

YHEN(nst) 

YHFAC(all) 

YHSUR(all) 

YMALED(all) 

YMALE(all) 

YMSUR(all) 


DEFINITION 

total  number  of  sage  grouse  hens  March  15 

total  number  of  male  sage  grouse  March  15 

pounds/acre  of  each  plant  type  ingested  by 
deer  each  week 

overwinter  survival  rate  for  hens 

overwinter  survival  rate  for  males 

X-axis  of  Figure  6.3 

X-axis  of  Figure  6.2  (females  health  index) 

X-axis  of  Figure  6.2  (fawns  health  index) 

shrub  cover  (X-axis  of  Figure  6.9) 

forb  and  grass  biomass  (X-axis  of  Figure  6.8) 

X-axis  of  Figure  6.2  (males  health  index) 

shrub  cover  (X-axis  of  Figure  6.10) 

#  fawns/female  as  function  of  health  index 
(Y-axis  of  Figure  6.3) 

female  winter  survival  as  function  of  health 
index  (Y-axis  of  Figure  6.2) 

fawn  winter  survival  as  function  of  health 
index  (Y-axis  of  Figure  6.2) 

initial  conditions  yearling  hens  on  5 
identified  populations 

hen  survival  modification  factor  as  function 
of  shrub  cover  (Y-axis  of  Figure  6.9) 

hen  survival  over  March  15  to  September  30 
as  function  of  forbs  and  grasses  (Y-axis 
of  Figure  6.8) 

male  winter  survival  as  function  of  health 
index  (Y-axis  of  Figure  6.2) 

initial  conditions  yearling  males  on  5 
identified  populations 

male  survival  March  15  to  September  30  as 
function  of  %  shrub  cover  (Y-axis  of 
Figure  6. 10)    ,,  _ 

00020 


VARIABLE 
YNESMN 

YNX(all) 
YNY(all) 


DEFINITION 

minimum  yearling  nesting  success  if  forbs 
density  too  low 

forb  and  grass  biomass  (X-axis  of  Figure  6=6) 

yearling  nesting  success  as  function  of  forbs 
and  grasses  average  density  (Y-axis  of 
Figure  6.6) 


r\  n  n  o  ' 


G 


3.   PROGRAM  LISTINGS 

The  computer  program  listings  are  in  three  sections: 

3.1)  FORTRAN  source  code; 

3.2)  data  files;  and 

3.3)  common  block. 

The  model  is  structured  hierarchically  (Figure  1) ,  and  runs 
under  a  simulation  language  called  SIMCON  (Hilborn,  1973)  .   When  the 
model  is  run,  SIMCON  first  calls  an  initialization  subroutine  (UINIT) 
to  read  in  the  common  block  and  allocate  storage  (performed  by  the 
subroutine  CMREAD)  and  store  the  variables'  initial  values  (performed 
by  the  subroutine  DFAULT) .   The  main  line  subroutine  UMODEL  is  called 
during  each  time  step,  which  in  turn  calls  each  of  the  submodel 
subroutines.   Table  1  shows  the  relationship  between  subroutines,  data 
files  and  common  blocks. 


r.  (\f\00 


SIMCON 


t 

i 

SAVAI 

DAVE . C 

BEEF.C 

WAT.C 

WILD.C2 

CHOW.C 

GROW . C 

SAVAL.D 


GROW.D 


CHOW.D 


DAVE . DD 


BEEF . D 


DAVE.D 


DAVE . SET 


WILD.D2 


WAT.D 


-   file  names 


subroutine  names 


■*■     -      subroutine  calls 


data  input 


-   connected  files 


Figure  1:   Structure  of  the  Saval  Ranch  integrated  model.   Files 
and  subroutines  are  listed  in  the  order  they  appear  in 
the  program. 


r\  n  n  o  ?. 


0 


Table  1:   Data  files  and  common  block  files  associated  with  each 
submodel. 


SUBMODELS 

SUBROUTINE 
NAMES 

DATA  FILES 

COMMON  BLOCK 
FILES 

Mainline 

UMODEL 

SAVAL.D 

SAVAL . C 

Vegetation 

VEG 

DAVE.D 
DAVE.DD 

DAVE.C 

Hydrology 

WAT 

WAT.D 

WAT.C 

Cattle  and 
Economics 

VIANDE 

EATUM 

OBESE 

BEEF . D 
CHOW.D 
GROW.D 

BEEF . C 
CHOW . C 
GROW.C 

Wildlife 

WILD 

WILD. D 2 

WILD.C2 

D 
D 
1 


r)  c\  n  ?  4 

"  .'  ■>:,■'  \J    *—"■ 
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3.1  FORTRAN  SOURCE  CO  D 


INITIALIZATION  ROUTINE 

This    routine    reads    in    the    common    block    for    all    FORTRAN    subroutines, 
and    initializes    all    variables    to    the    values    listed    in    the    data    filei 

1  SUBROUTINE    UINIT 

2  CA  L  L  CM  R  E A  D  (  '  SAUA  L  .  C  '  ) 

3  CALL  DFAULT ( ' 1 =SAUAL . D  ' ) 

4  RETURN 

5  END 

MAINLINE    SUBROUTINE 


This  subroutine  is  called  every    time  step  (week).  It  then  call* 
each  of  the  other  subroutines  in  the  order  shown. 

6  SUBROUTINE  UMODELC ITIME) 


t 


♦  CONTINUE  WITH  SAv'AL.C  RETURN 
7.05         IFC ITIME. GT.l)  GO  TO  300 
7.07         I FCIMUNCH.LT. 2)  GO  TO  30  0 

7.1  TPREF=0. 
7.15         DO  900  IR=1,3 

7.2  DO  900  IF=1,7 
DO  90  0  IW=1 ,2 


7.3  TPREF=TPREF+FPREF(IF)*NPREF(IW)*RPREF(IR) 
7.35  90  0   CONTINUE 

7.4  DO  950  IR=1,3 
7.45  DO  950  I F=l ,7 

7.5  DO  950  IW=1 ,2 

7.55  PREFF(IF,IR, IW)=(FPREF<IF)*WPREF(IW)*RPREF(IR))/TPREF 

7.6  950   CONTINUE 
7.65  30  0   CONTINUE 

8  IF(ITIME.EQ.1)IWEEK=0 

8.2  I  F(  ITIME.  EQ.  LAND.  I  LUMP.  EQ.l)  CALL    LUMP 

9  INEEK=IWEEK+1 

10  IF(IWEEK.EQ.39)IWEEK=1 

11  IM0NTH=(IWEEK+l)/4+l 

11.05  C  NRITE(3,100)     IWEEK 

11.1  100    F0RMATCI3) 

11.2  IF(INAT.GT.0)CALL    WAT( ITIME) 

11.6  C  WRITE(3,100)     IWEEK 

12  IF(IUEG.GT.0)CALL    VEG(ITIME) 

12.5  C  WRITE(3,100)     IWEEK 

13  IF( IBEEF.GT.0)CALL    VIANDE 

13.1  C  WRITEC3,100)     IWEEK 

13.2  IFCIWILD.GT.0)CALL    WILD  00^25 

13.6  C  WRITE(3,100)     IWEEK                                  '"  " v 

14  RETURN 
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0 


0 

1 
1 

D 

G 
Q 


15 

END 

16 

* CONTINUE 

WITH 

VEG.S  RETURN 

17 

♦CONTINUE 

WITH 

BEEF.S  RETURN 

15 

♦CONTINUE 

WITH 

WAT.S  RETURN 

20 

* CONTINUE 

WITH 

WILD.S2  RETURN 

21 

♦CONTINUE 

WITH 

CHOW.S  RETURN 

22      ^CONTINUE  WITH  QROW.S  RETURN 
U  E  G  E  T  A  T  I  0  N    SUBMODEL 


This  submodel  is  described  in  section  3  of  the  report. 

1  SUBROUTINE  MEG (1 TIME) 

2  * CONTINUE  WITH  SAUAL . C  RETURN 

3  IF(ITIME.NE.l)GO  TO  IS 

4  IF(ILUMP.EQ.O}GO  TO  5 

5  CALL  LUMP 

6  GO  TO  18 

7  C 

8  C  CHANGE  RA  FROM  TOTAL  PRODN  TO  BIOMASS  Ai  START 

9  C  OF  SEASON.   SET  INITIAL  ROOT  BIOMASS. 

10  C 

11  5  DO  10  I=1,NPAS 

12  DO  10  J=1,NRT 

13  DO  10  IF-1 ,7 

14  IF(RA(I , J, 8) .EQ.0)GO  TU  10 

15  ROOT(I , J,IP)=0 

16  IF(IP.GT.b)GO  TO  10 

17  ROOTCI ,J,IP)*RPB*RA(I ,J,IP) 

18  10     CONTINUE 

19  C 

20  C  IMPLEMENT  ACTIONS 

21  C 

22  18  IF(IUM.EQ.0)GO  TO  45 

23  DO  42  I=1,NPAS 

24  IF(IOC(I) .EQ.0)GO  TO  42 

25  DO  40  J=l ,NRT 

26  DO  40  K=l,4 

27  IF(VMANIP(I , J,K) -EQ.0)GO  TO  40 

28  IF(K.EQ.l .AND.^MANIPCI , J,K) .LE.0)GO  TO  35 

29  PREM=UMAN I P ( I , J , K ) /RA ( I , J , 8 ) 

30  C 

31  C  SHRUB  KILL  FOR  BURNING  OR  SPRAYING 

32  C 

33  IF(K.EQ.4)G0  TO  25 

34  DO  20  IP  =  1  ,3 

35  IF(K.GT.1)RA(I , J , I P) =RA( I , J , I P)*(l . -SHKI LLC K)*PREM) 

36  IF(K.Eg.l)YSEEDCIP)=RACI , J , I P)*( 1 . -SHKILL< K)*PREM) 

37  20  CONTINUE 

38  C 

39  C  REMOUE  LAND  IF  PLOUGHING/SEEDING  OR  MINING. 

40  C  KEEP  TRACK  OF  AMT .  OF  LAND  IN  SEEEDING  RESERVE. 

41  C 

42  25  IF(K.EQ.2.0R.K.EQ.3)G0  TO  40                  00026 
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43  IF(K.EQ.1)RA(  I  ,  J,14)=RA(  I  ,J,S)*PREM 

44  RAC I , J , 8) =RA( I , Jf8)*<i . -PREM) 

45  GO  TO  40 

46  C 

47  C  ADD  LAND  BACK  INTO  RA  IF  SEEDING  READY  FOB.  GRAZING. 

48  C  USE  WEIGHTED  AVERAGE  OF  PLANT  BIOMASSES,  AND 

49  C  FIXED  MID-SUMMER   BIOMASS  IN  SEEDED  AREA. 

50  C 

51  35  DO  33  IP=1  ,7 

52  ARET  =  -1 . *VHAN I P ( I  .  J  ,  1) 

53  I F ( RA ( I , J , 8 ) . LE . 0 ) GOTO  38 

54  RAC I , J, IP)=(RA( I ,J,IP)*RA( I ,  J,8)  + 

55  1  USEED  C I P ) *ARET ) / ( RA ( I , J , 8 ) +ARET ) 
5h  38           CONTINUE 

57  RA( I ,J,14)=RA(I , J,14)-ARET 

53  RA(I , J,S)=RA(I , J,S)+ARET 

53  40     CONTINUE 

SO  42     CONTINUE 

61  C 

62  C  HAY  PRODUCTION.  K=PASTURE  (TREM .  , GANCE. , HAYSTACK) 

63  C 

64  45     IF(IWEEK.NE.20)GO  TO  100 

65  THAYP=0 . 

6S  DC  60  K  =  l  ,3 

67  I=3+K 

68  DO  50  M=l,5 

69  50        H  Y  ( M )  =  H  A  Y  Y  (  HAY  OP  (  K  )  ,  M ) 

70  I F ( K . EQ . 1 ) WPA ( 1 ) =FTOT ( 2 ) /AMEAD ( 1 ) 

71  I F ( K . EQ . 2 ) WPA ( 2 ) =FTOT ( 1 ) /AMEAD  <  2 ) 

72  I F ( K . EQ . 3 ) WPA ( 3 ) =FTOT ( 3 ) /AMEAD ( 3 ) 

73  HAYP C  K ) =SLP ( WPA ( K )  , HX , HY , 5 ) *AMEAD( K ) 

74  T  HAY  P =T  HAY  P+  HAY  F  (  K ) 

75  60     CONTINUE 

77  C  GROWTH  AND  MORTALITY  OF  MAJOR  PLANT  TYPES. 

73  C  TREAT  CHEATGRASS , PERENNIAL  GRASSES  AND  FORBS, SHRUBS 

79  C  AS  THREE  GROUPS. 

80  C 

81  100    DO  475  I=1,NPAS 

82  IF(IOC( I) .EQ.0)G0  TO  475 

83  C 

84  C  IF  SEED  SET,  TOTAL  UP  PASTURE  FORAGE  PRODN 

85  C  PERENNIAL  GRASSES  AND  FORBS  ONLY. 

86  C 

87  IF(IWEEK.NE.IFDATE)GO  TO  102 

88  GTPROD(I)=0 

89  DO  101  IP=4,6 

90  101    GTPRODC I)=GTPROD(I )+TPROD( I , IP) 

91  102       IF(IWEEK.NE.l)GO  TO  105 
32  DO  10  3  IP=1 ,7 

93  103       TPROD( I  ,  IP)=0 . 

94  105      TRED=SLP ( TMAX ( I ) , GTEMPX , GTEMPY , 3 ) 
35  DO  450  J=l ,NRT 

96  IF(RA(I ,J,8) .EQ.COGO  TO  450 

37  c  _       -\  n  n  o  "• 

98  C    COPY    MOISTURE    X-AXIS    APPROPRIATE    TO    RANGE    TYPE  •s\J\J£.t 


I 
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39  C  FROM  MASTER 

1100  DO  110  K  =  l,5 

101  110       G  S OX ( K ) =  G  S OXM ( J , K ) 

1 0  2  C 

110  3  C  COMPUTE  TOTAL  PLANT  EI0MAS3  ON  RANGE  SITE. 

104  C 

105  TOTPB( I , J)=0 

106  DO  115  IP=1 ,7 

107  TOTPBC I , J)=TOTPB(I , J)+RA(I , J, IP) 

108  115       CONTINUE 
103  DO  40  0  IP=1,7 

110  IF(RA(I , J, IP) .EQ.0)GO  TO  400 

ill  IF( IWEEK.EQ.38)G0  TO  30  0 

112  WMAX=SO( I , J,l) 

113  NMIN=AMIN1 ('  SO( I , J , 1 ) , SO( I  , J , 2) , SO( I  , J , 3) ) 

114  C  CHEATGRASS  -  USES'  MOISTURE  IS  UPPER  10"  ONLY 

115  IF(IP.EQ.7)Q0  TO  120 

US  I  F (  S 0 (  I  ,  J  ,  2 )  .  GT  .  WMAX )  WMAX = S  G  (  I  ,  J  ,  2 ) 

117  C  GRASSES  AND  FORES  -  USE  MAX.  OF  MOISTURE  IN  0-10"  AND  10-20" 
US  IF( IP.GE.4.AND. IP.LE.6)G0  TO  120 

113  C  SHRUBS  -  USE  MAX.  OF  MOISTURE  IN  ALL  3  SOIL  LAYERS 

120  I F ( SO ( I , J , 3 ) . GT . WMAX ) NMAX=SO (I , J , 3 ) 

121  C 

122  C  COMPUTE  POTENTIAL  PLANT  GROWTH  RATE  WITH  NO  COMPETITION 

123  C 

124  120  GMRED ( I , J , I P ) =SLP ( WMAX , GSOX , GSOY , 5 ) 

1 25  GPLANT ( I  , J , I P ) =GMAX (IP) *TRED*GMRED ( I , J , I P ) 

126  C 

127  C  RECOMPUTE  PLANT  GROWTH  RATE  TO  ACCOUNT  FOR  COMPETITION. 
12S  C  POTG=POTENTIAL  GROWTH  (LB/A) 

129  C  REDUCE  PLANT  GROWTH  RATE  IF  NOT  ENOUGH  WATER  TRANSPIRED, 

130  C  BUT  LEAVE  GROWTH  RATE  AS  IS  IF  WATER  SURPLUS. 

1 31  C 

132  RQEVIN=RA(I , J, IP)*RQEULB( I F)*0 .000276 

1 33  AVEY I N=PEYAF ( I , J )*RA C  I  ,  J  ,  I  P ) /TOTPB ( I , J ) 

134  GCOMP(I , J,IP)=0. 

135  I F  (  RQEY I N  .  NE  .  0  )  GCOMP  (  I  ,  J  ,  I  P )  =AUEU  I IM/RQEV I N 

136  IF(AUEYIN.GE.R0EVIN)GO  TO  140 

137  GPLANT ( I , J , I P ) =GPLANT ( I , J , I P ) *AUEU I N/RQEU I N 
133  140  CONTINUE 
133                      PMR=SLP ( GPLANT ( I , J , I P ) , PMX , PMY , 2 ) 

140  C 

141  C  KILL  SOME  SHRUBS  IF  SOIL  MOISTURE  IN  ALL  SOIL  LAYERS 

142  CIS  GREATER  THAN  HMMR  %  (30%????) 

143  I F( I P . LE . 3 . AND . WMIN . GT . HMMR) PMR=1 . -GMRED( I , J , I P) 

144  C  BIOMASS  CHANGE  WITH  NO  GRAZING...  PRODUCTION  ACCUMULATED 

145  C  ONLY  UNTIL  SEED  SET 
14b  DPB=(GPLANT(I , J , I P) -PMR)*RA( I , J, IP) 

147  IF(IWEEK.LE.IFDATE)TPROD(I ,IP)=TPROD(I ,IP)+DPB 

148  C  BIOMASS  CHANGE  WITH  GRAZING... 
143  DPB=DPB-CE( I ,J,IP)*(1 .+CWST)-WE( I » J, IP) 

150  CE( I , J, IP)=0 .0 

151  WE (I , J, IP)=0 . 

152  C 

153  C  ADJUST  DPB  WITH  TRANSLOCATION  TO  OR  FROM  ROOiS.  ,, .^  n  «  n 

154  C  BIOMASS  FROM  SHOOT  TO  ROOT  AFTER  FDATE,  FROM  :\j'JCO 
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155 
156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

1  by 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

20  4 

20  5 

20  6 

20  7 

20  8 

20  9 

210 


C    ROOT    TO    SHOOT    BEFORE .ADJUST    TOTAL    PRODN . 

C    CHEATGRASS    HAS    NO    ROOT    STORAGE 

C 

IF( IP.GT.6)G0   TO    250 

DROOT=0 . 

PNGR=DPB/RA( I , J „ IP) 

I F ( I WEEK- 1 FDATE ) 1 60 , 1 60 , 20  0 

L  «en»«BKian*a«DBiiKiiKocBB«osonrjftDO»coQ"aijiiGt;EiL' 
l-i  t  ■  .  t  t  «  i  bBlFOKL  bbLL'  OC.l>*caee>c*siDOBSBDos>G 
L -    BBB0BBB8BBaaaaBBBeeBB»BB»ae*a«JBBeanBBBQBBBB 

160  TPRTS=PRTS*C 1 . -FLOAT ( I WEEK) /FLOAT ( I FDATE) ) 

DROQT=-l  .-A-ROOTC  I  ,  J,  IP)*TPRTS 

C 

C    IF    STANDARD    RATE    OF    TRANSFER    (TPRTS)    SUFFICIENT 

C   TO   MAINTAIN    BIOMASS,    DON'T    DECREMENT    RESERVES    FURTHER. 

C 

IF(DPB-DROOT.GE.0)GO   TO    ISO 

I F( ROOTC I , J , I P)+DPB . GE . 0 ) DROOT=DPB 

I FC ROOTC  I  , J,IP)+DPB.LT.0)DROOT=-0.5*ROOi  (I  ,  J,  IP) 

C   MAKE    CHANGE    IN   ABOVE    GROUND    BIOMASS   MORE    POSITIVE. 

C    (DROOOT    IS   NEGATIVE) 

C 

180  DPB=DPE-DROOT 

TPRODC I f I P) =TPROD( I, IP) -DROOT 
60   TO    230 

Lo>Booo£jBBBfiraoaa*aaoiiioaooflooooBODnoQDOeoeoBEsa 
C  ■  b  c  a  b  b  a  a  *  A  FT  E  R  SEED  O  ET  eD«a.B»BB«08BB«,BDB«Be 
WaBBoaBBaBBBeBBSUoeiscaaBeBaa      BBeoaBHoBauo      oano 

200  DROOT=RA( I , J, IP)*PSTR 

C 

C   ACCOUNT    FOR    ROOT    TO    SHOOT    OR    SHOOT    TO    ROOT   TRANSFERS 

C 

230  ROOTC I , J,IP)=ROOT(I , J , I P)*( 1 . -PMR)+DROOT 

I F( ROOTC I , J, IP) .LT.0)ROOT(I , J , IP) =1.0 

I FC ROOTC I , J, IP) .GT.ROOTMXC IP)) 
1  ROOT C .  I , J , I P ) =ROOTMX (IP) 

C 

C    CHANGE    PLANT    BIOMASS.       RESET    IF    WENT    BELOW   ZERO    OR 
C    OUTRAGEOUSLY    HIGH. 
C 
250  RACI ,JSIP)=RACI s J, IP)+DPB 

I  F  C  RA  C  I  ,  J  ,  I  P  )  .  LT  .  0  . )  RA  (  I  ,  J  ,  I  P )  =1  .  0 

I FC  RAC I , J , I P) . GT . PLNTMXC IP)) 
1  RACI , J,IP)=PLNTMXCIP) 

GO   TO    400 
C 

C    PLANT    AND    ROOT    WINTER    SURVIVAL. 
C 
300  RACI ,J,IP)=RACI ,J,IP)*PWS(IP) 

IFC IF.GT.6)G0    TO    40  0 

ROOTC I , J, IP) "ROOTC I  *  J, IP)*RWS 
40  0  CONTINUE 

450  CONTINUE 

475  CONTINUE 

c  o  n  n  2  ° 

C    CALCULATE    %    COMPOSITION    BY    WT .  J u w       ' 
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211  C 

212  DO  530  I=1,NPAS 

213  IF(IOC( I ) .EQ.0)GO  TO  530 

214  GFTOT=0 . 

215  DO  473  IP=1,7 

216  478    GPCWT ( I  , I P ) =0  . 

217  DO  50  0  J=l ,NRT 

21 8  I F ( RA ( I , J , 8 ) . EQ . 0 ) GO  TO  50 0 

219  FTOT=0  . 

220  DO  480  IP=1 ,7 

221  GPCWT < I , IP}=GPCWT(I , IP)+RA(I , J, IP) 
480       PTOT=PTOT+RA ( I , J , I P ) 

GPTQT=GPTGT+FTOT 


G 


iiii 


<££w 


£24  DO  430  IP=1 ,7 


£c 


PCWTCI , J,IP)=0 

226  IF(PTGT.NE.O)PCWT( I ,J,IP)»RA(I , J , I P)/P I  0 l *100 

430  '   CONTINUE 


'd'cL'd< 


2  3  6 


266 


50  0    CON  I INUE 


223  DO  520  IP=1,7 

230  I  F (  GPTOT  .  NE  .  0  .  )  GPCWT (  I  ,  I  P )  =GPCNT £  I  ,  I  P) /GPTO' 


231  520    CONTINUE 

232  530    CONTINUE 


:33      C 

<;34      C  CALCULATE  %  COVER 


a 


C 


c 


CALL  COVER 


238  C  %  PROTEIN  CONTENT. 

233  C 

240  DO  60  0  IP=1,7 

241  IF(IP.LE.3)VP< IP)*SLP<IWEEK»SMPX,SyPY»2) 

242  I F ( I P . GE . 4 . AND . I P . LE . 6 ) VP ( I P ) =SLP ( I WEEK ,  PVPX , PVPY , 2 ) 

243  I F  < I P . EQ . 7 ) V P ( 7 ) =SLP ( I WEEK , CVPX , CMPY , 2 ) 

244  600    CONTINUE 

245  RETURN 

246  END 

247  SUBROUTINE  COVER 

248  ^CONTINUE  WITH  SAVAL . C  RETURN 
243  C 

250  C  THIS  SUBROUTINE  SCALES  %    COVER  OF  GRASSES , SHRUBS , 

251  C  TREES , CANOPY  AND  GROUND  ACCORDING  TO  %    WT ,  A 

252  C  TRULY  FALSE  ASSUMPTION. 

253  C 

254  DO  106  I=1,NPAS 

255  IF(IOCCI) .EQ.0)GO  TO  106 

256  DO  100  J=l ,NRT 

257  IFCRACI , J ,8) .EQ.0)QO  TO  100 

258  PC(I , J,1)=0. 
253  C  %    GRASS  COVER 

260  I F ( PWTMAX ( J , 2 ) . EQ . 0 ) GOTO  20 

261  PC (I  , J , 1)  =  ( PCWT ( I , J , 4 ) +PCWT ( I  ,  J , 5 ) +FCWT ( I  , J , 6 ) 

262  1  +PCWT ( I , J , 7 ) ) /PWTMAX ( J , 2 ) *PCMAX2 ( J , 3 ) 

263  C  %    SHRUB  COVER 

264  20     PC( I  , J,2)=0  . 


I F ( PWTMAX ( J , 1) . EQ . 0 ) GO  TO  40 

PC  (I  ,  J,2)=(PCWT(I  ,J,1)+PCWT(I  ,  J,2)+PCWT(I  ,  J,  3))/     U  U  0  J  C 


I 
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1  PWTMAX ( J , 1 ) * P CM AX  2 ( J , 1) 

C  %    TREE  COVER 

40       PC( I , J,3)=0  . 

C  ************  NEXT  STATEMENT  TEMPORARILY  REMOVED  ****** 

C         I F ( J . EQ . 3 ) PC ( I , 9 , 3 ) =PCMAX2 ( 5 5 1 ) 

C  %    CANOPY  COMER 

PCCI , J,4)=PC(I , J,2)+PC(I ,J,3) 

C  %  GROUND  COMER  ( GRASS , LITTER , ROCK  AND  CANOPY) 
275      C 
27b      C  RLF=RATIO  OF  ( GRASS+LIT I ER) : URASS  -  40%? 

277      C 

£73  PC  C I  , J . 5 )  =PC  (  I  ,  J  ,  1)  *RLP+PCMAX2  (  J ,  4)  -i-PC  (I  ,  J ,  4 ) 

275      C 


£  b  ,-■ 

£33 

233 

270 

271 

272 

273 

ii  7  4 

•i  c  U 


^Bb 


iCOO 


C  N  0  RM A  L 1 2  E  " P  E  R  C ENT A G E  S "  . 

231  C 

232  TOTPCT=PC(  I  ,  J,5)+PCt:  I  ,  J,4) 

233  I FCTOTPCT.LT .100=) GO  TO  100 

234  PC ( I , J , 4 ) =PC ( I , J , 4 )  /"T0TPCT*1 0  0 

235  PCi I , J,5)=PC(I , J,5)/TOTPCT*100 . 
100       CONTINUE 
10  3     CONTINUE 

RETURN 
233  END 

290  SUBROUTINE  LUMP 

231  ^CONTINUE  WITH  SAVAL . C  RETURN 

232  C 

233  C  THIS  SUBROUTINE  GROUPS  THE  TEN  SITES  IN  THE  SET  FILE  TO 

234  C  FORM  NRT  RANGE  SITES.  ROW  K  OF  NEWRS  SPECIFIES  THE  OLD 

235  C  RANGE  SITES  THAT  ARTE  GROUPED  INTO  THE  KTH  NEW  RANGE 

236  C  SITE  CATEGORY.  AT  THE  SAME  TIME  THIS  SUBROUTINE 
297  C  REDUCES  MID-SUMMER  BIOMASS  TO  EARLY  SPRING.  NEW 
293  C  RANGE  SITE  ARRAYS  ARE  WRITTEN  OUT  TO  UNIT  3  TO 
293  C  FORM  A  NEW  SET  FILE. 
30  0  C 

301  C  LOOP  OVER  PASTURES... 

302  C 

30  3  DO  100  I=1,NPAS 

30  4  C 

30  5  C  LOOP  OVER  NEW  RANGE  SITES 

30  6  C 
30  7  DO  SO  JNEW=1,NRT 

303  RSANEW=0  . 
303  C 
310  C  LOOP  OVER  OLD  RANGE  SITES  THAT  MAKE  UP  NEW  RANGE  SITE 

OJ.1  (_• 

312  DO  40  JO=l ,10 

313  J  0  L  D  =N  EW  R  3  (  JN  EW  ,  J  0  ) 

314  IF( JOLD.EQ.0)GO  TO  45 

315  RSAOLD=RA( I , JOLD , 3) 
31b  IFCRSAOLD.EQ.O . ) GO  TO  40 
317                 R  SAN  EW = R  SAN  EW+ R  S A 0 L  D 
313  C 
315  C  LOOP  OVER  PLANT  TYPES  CALCULATING  WEIGHTED  SUMS  OF  PLANT 

320  C  BIOMASS.  RANEW  IS  INITIALIZED  TO  ZERO  IN  SET  FILE. 

321  C 

DO  10 IP  =  1 


™ny, 
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OiO 


331 
332 


RANEW ( I , JNEW , I P) = RANEW ( I , JNEW , I P ) + 

324  1  RAC I , JOLD, IP)*RSAOLD 

325  10  CONTINUE 

326  C 

327  C  LOOP  OVER  SOIL  MOISTURE  LAYERS  CALCULATING  WIGHTED  SUMS 
32S  C  OF  AVAILABLE  WATER  CAPACITY. 

329  DO  20  IS=1,3 

3  0  AW  CN  EW  (  I  ,  JN  EW  ,  I  S )  =A  W  CN  E  W  (I  ,  JN  E  W  ,  I  S )  + 

1  AWC <  JOLD , I S)*RSAOLD 

20  CONTINUE 

333  40        CONTINUE 

334  C 

335  C  CONVERT  WEIGHTED  BIOMASS  SUM  TO  WEIGHTED  AVERAGE  SPRING  ETOMASS 
33b  C  CONVERT  WEIGHTED  AWC  SUM  TO  WEIGHTED  AVERAGE  AWC. 

337  45         DO  50  IP=1,7 

338  IF(RSANEW.EQ.0)GO  TO  50 

3  3  9  RAN  EW  (  I  ,  JN  EW ,  I P  )  =SP  B*  <  RAN  EW  (  I  ,  JN  EW  ,  I  P )  /  R  SAN  EW )  ** 

340  1  SPBEXP 

341  50         CONTINUE 

342  DO  55  I S=l , 3 

343  IFCRSANEW.EQ.6)GO  TO  55 

344  AW  CN  EW  (  I  ,  JN  EW  ,  I  S  )  ■  A  W  CN  EW  (  I  ,  JN  EW  ,  I  S  )  /  R  SAN  EW 

345  55         CONTINUE 

346  C 

347  C  WRITE  OUT  NEW  VALUES  IN  SET  FILE  FORM. 

348  C 

343  WR I TEC  3 , SO )  I  , JNEW , C  RANEW ( I , JNEW , I P ) , I P=l , 7 ) , RSANEW 

350  60        FQRMATC'S  RAC  ,12,"  ,'  ,U,'  .ALL)*'  ,8F6.0) 

1  WR I TE ( 3 , 70 ) I , JNEW , ( AWCNEW ( I , JNEW , I S) , I S«l , 3 ) 

70        FORMAT (  ;  S  AWCNEW (  '  ,  I  2  ,  '  ,  '  ,  1 1  ,  "  ,  ALL )=''  f  3F5  .  2 ) 

353  80        CONTINUE 

354  100   .  CONTINUE 

355  RETURN 

356  END 


35^ 


HYDROLOGY  SUBMODEL 

This  submodel    is    described    in    Section    4    of    the    report. 

1  SUBROUTINE  WAT (ITIME) 

2  ^CONTINUE  WITH  SAVAL . C  RETURN 

3  C 

4  C  SETUP  INITIAL  CONDITIONS  IN  FIRST  RUN 

6  IFCITIME.NE.l)  GOTO  10 

7  DO  6   1=1,  NPAS 

8  DO  5   J=l ,  NRT 

9  DO  4   K=l ,  3 

I  o  SOWAT  U , J , K )  =  0 . 9*AWCNEW ( I  , J , K ) *WATRES ( K ) 

II  4  CONTINUE 

12  IF  ( (I .NE.S) .AND. ( I .NE.9))  GOTO  5 

13  SNOPACY  I  ,  J)=l  . 

1 4  FREZDP  (I  ,  J )  =6 .  r  r\  n  7 ") 

15  5  CONTINUE  UUUJu 
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16 

6 

CONTINUE 

17 

10 

CONTINUE 

IS 

NSTEP  =  NSTEP  +  1 

13 

NDAYS  =  7 

20 

IWFLG  =  0 

21 

IF  (NSTEP. NE. 2)  GOTO  30 

C 

ZERO  ACCUMULATORS  AT  START  OF  YEAR 

22.2 

TSEDG  =0.0 

22.4 

TSEDM  =0.0 

23 

DO  20   1=1,  3 

24 

FTOTC I )  =  0 . 

25 

20 

CONTINUE 

26 

c 

27 

c 

ZERO  WEEKLY  ACCUMULATORS 

28 

30 

DO  60   I =1 ,  NPAS 

23 

TMAX(I)  =  0. 

30 

60 

CONTINUE 

31 

DO  65   1=1,  3 

32 

FWTOT(I)  =  0. 

33 

65 

CONTINUE 

34 

NQUAL  =  0 . 

35 

DO  80   1=1,  NPAS 

36 

DO  75   J=l ,  NRT 

36.2 

REMAP (I , J)=0 . 

37 

DO  70   K=l ,  3 

38 

SO (I , J,K)  =  0  . 

33 

70 

CONTINUE 

40 

75 

CONTINUE 

41 

80 

CONTINUE 

42 

DO  30   1=1,  NPAS 

43 

PRECIP(I)  =  0, 

44 

50 

CONTINUE 

45 

IF  < NSTEP. LE.NWEEKS)  GOTO  10  0 

46 

c 

47 

c 

SET 

UP  WINTER  STEP 

48 

NDAYS  =  WINDAY 

43 

IWFLG  =  1 

50 

NSTEP  =  0 

51 

c 

52 

c 

START  MAIN  LOOP  OVER  DAYS  IN 

53 

c 

54 

c 

55 

100 

DO  5000   IDAY=1,  NDAYS 

56 

J DAY  =  J DAY  +  1 

57 

IF  ( JDAY.GE.366)  JDAY  =  1 

58 

c 

53 

c 

INI 

TIALIZE  RAIN  PER  PASTURE 

60 

DO  110   1=1,  NPAS 

61 

DFREC(I)  =  PREC(JDAY) 

62 

110 

CONTINUE 

63 

DPREC(S)  =  PREC(JDAY)  *  PMULT 

64 

' 

DPRECC3)  =  PREC(JDAY)  *    PMULT 

65 

c 

66 

c 

67 

c 

START  MAIN  LOOPS  OVER  PASTURE 

63 

L 

THE  CURRENT  TIME  STEP 


S  AND  RANGE  TYPES         "  0 G 3  J 


D 
I 
I 

0 
I 

0 
G 
I 
I 

D 
I] 

n 
o 

Q 
D 
G 
I 
I 
I 
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69  C 

70  DO  100  0  IPAS-1,  NPAS 

71  DO  90  0   IRTYPE=1,  NRT 

72  I F  ( RA ( I  PAS , I RTYPE ,  S )  . EQ . 0 )  GOTO  90 0 

73  1 F ( ( I OC ( I  PAS )  . EQ . 0 ) . AND . ( I  PAS . NE . 9 ) )  GOTO  90  0 

74  C 

75  C  EU A POT  RAN  S  P I  RAT I  ON 

76  C 

77  C  POTENTIAL  EVAPOTRANSFI RATI  UN 

78  RADLAN  =  20  0.  +  60  0. *8 I N ( ( JDAY/365 . ) *3 . I 41 53) 

79  I F  ( LUPTR ( I  PAS )  . EQ . 1 )  DMTEMP  = 

50  1   CCTMAXU JDAY)*TMULT)  + 

51  2    CTMINUJDAY)*TMULT))  /  2. 

82  I F  ( LUPTR ( I  PAS ) . EQ . 0 )  DMTEMP  =  ( TMAXL ( JDAY ) +TM I NL ( JDAY ) ) /2 . 

83  ET  P  OT  =  <  (  0  .  0 1 4*  DMT  EM P  -  0.37)* R A  D  LAN )  /  5  8  0 

84  I F  ( ETPOT . LT . 0 . )  ETPOT  =  0 . 
b!3  C 

86  C  PERCENT  ACTIVE  CANOPY 

87  PACAN  =  S I N ( ( ( C FLOAT ( JDAY ) ) -66 . ) /3£5 . )*3 . 1 41 53 ) 

88  I F  ( ( JDAY . LT . 66 ) . OR . ( JDAY . GT . ( NNEEKS*7+66 ) ) )  PACAN  =  0 . 

89  C 

90  C  ACTUAL  ET  BY  PLANTS 

31  I F  ( ( I T I ME . NE . 1 )  . AND . ( I OC ( I  PAS )  . EQ . 1 ) )  PCDUM=PC ( I  PAS ,  I  R  I YPE 

) 

92  IF  ( ( I T I ME . EQ . 1 )  . OR . ( I OC ( I  PAS ) . EQ . 0 ) )  PCDUM=PCMAX ( I R  i YPE , 5 ) 

33  ETPLNT  =  PCDUM  *  TCOF  *  ETPOT  *  PACAN  /  100. 

33.1  ETSUM  =  0. 

33.2  DO  210   1=1 ,  3 
ETSUM  =  ETSUM  +  ( PACAN*SOWAT( I  PAS , I RTYPE , I ) )  / 

,4  1           <  AWCNEW ( I  PAS , I RTYPE , I ) *WATRES ( I ) ) 

93.5  210  CONTINUE 

33.6  ETPLNT  =  ETPLNT  *  ETSUM  /  3. 

94  C 

95  C  ET  FROM  BARE  SOIL 
ETSO I L  =  (10  0. -PCDUM ) * ( SOWAT ( I  PAS , I RTYPE  , 1 ) /WATRtS ( 1 ) ) 

1  *  ETPOT  /  (AWCNENi:iPAS?IRTYPE.,l)  *  10  0.) 

98  C 

99  C  REMOVE  ET  FROM  SUlL  LAYERS 

1 0  0  SOWAT ( I  PAS , I RTYPE , 1 ) =SONAT ( I  PAS , I RTYPE , 1 )  -  ETSO I L 

101  I F  ( SOWAT ( I  FAS , I RTYPE , 1 )  . LT . 0 . )  SOWAT ( I  PAS , I RTYPE , 1 ) =0  . 

101.2  ETPREM  =  ETPLNT 

10  2  DO   40  0   ILEAK=1,  3 

103  IF  ( ETPLNT. LT. SOWAT (I  PAS, I RTYPE, I  LEAK))  GOTO  300 

104  ETPLNT=ETPLNT  -  SOWAT ( I  PAS , I RTYPE , I  LEAK) 

105  SOWAT (I  PAS, I RTYPE, I  LEAK )=0  . 

106  GOTO  400 

107  300  CONTINUE 

108  SOWAT ( I  PAS , I RTYPE , I  LEAK) =SOWAT( I  PAS , I RTYPE , I  LEAK) -ETPLNT 
108.2  PEyAP(IPAS,IRTYPE)  =  ETPREM  +  PEVAPC I  PAS, I RTYPE) 

109  GOTO  50  0 

110  400  CONTINUE 

110.2  P  EVA  P ( I P A  S , I RT Y  P  E )  =  ET  P  R  EM  -  ET  F  LNT 

110.4  1    +  PEVAP ( I  PAS , I RTYPE ) 

111  50  0  CONTINUE  r                A 

±  cl  !_■                                                                ■ ;  -  v  «-»  » 

113  C                      SNOW  PACK  MANIPULATION 


93 

93, 


3  b 
97 


1 27 
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114  C 

115  IF  (DMTEMP.6T.32)  GOTO  60  0 

116  0 

117  C  ACCUMULATION 
US  SNOPAC ( I  PAS , I RTYPE ) =SNOPAC ( I  PAS , I RTYPE ) +DPREC ( I  PAS ) 

119  DPRECT <I PAS)  =  0. 

120  GOTO  70  0 

121  C 

122  C  MELTING 

123  60  0  CONTINUE 

124  SMELT  =  (DMT EM P  -  32)  /  10. 

125  I F  ( SMELT . LT . SNOPAC ( I  PAS , I RTYPE ) )  GOTO  61 0 
DPRECT ( I  FAS ) =DPREC ( I  PAS )  +  SNOPAC ( I  PAS , I RTYPE) 
SNOPAC < I  PAS , I RTYPE )  =  0 . 

128  GOTO  620 

129  610  CONTINUE 

1 30  SNOPAC ( I  PAS , I RTYPE ) -SNOPAC ( I  PAS , I RTYPE )  -  SMELT 

131  DPRECT(IPAS)  =  DPREC(IPAS)  +  SMELT 

132  620  CONTINUE 

133  70  0  CONTINUE 

134  C 

135  C  RUNOFF  /  INFILTRATION 

136  C 

137  C  RUNOFF 

1 38  APAR=CNA ( I PTR ( I RTYPE ) ) 
133  8PAR=CNB( I PTR( I RTYPE) ) 

140  CN  =  APAR  -  BPAR*PCDUM 

141  SMAX  =  (100  0./CN)  -10. 

141.1  WSOW  =  0 
141.12  AWSOW  =  0 

141.2  DO  720   1=1,  3 

141.3  HSOW "W S  OW+ S ON AT ( I  PAS , I RTY PE , I ) /WAT R E S ( I ) 
1 41 . 32  AWSOW=AWSOW+AWCNEW ( I  PAS , I RTYPE , I ) 

141.4  720  CONTINUE 

141.5  WSOW  =  WSOW  /  3. 
141.7  AWSOW  =  AWSOW  /  3. 

142  S=SMAX*( AWSOW  -  WSOW)  /  AWSOW 

144  RUNOFFC I  PAS , I RTYPE)  =  ( ( DPRECT < I  PAS) -0 . 2*S)**2)  / 

145  1     ( DPRECT ( I  PAS)  +  0.8*S) 

1 46  RUNOFF ( I  PAS , I RTYPE )  =  RUNOFF  < I  PAS , I RTYPE )  *  RFUDG 

147  IF  (DPRECT(IPAS) .LT.( 0.2*3))  RUNOFF CI PAS , I RTYPE) =0 . 

148  C 

143  C  INFILTRATION 

150  WAT  IN  =  DPRECT (I  PAS)  -  RUNOFF (I  PAS „ I RTYPE) 

151  IF  (WATIN.LT.O . )  WATIN=0. 

152  DO  800   IFIL-1 ,  3 

153  I F  < ( WAT I N+SOWAT ( I  PAS , I RTYPE , I F I L ) )  . LT . ( WATRES ( I F I L ) 

154  1  *AWCNEW ( I  FAS , I RTYPE , I F I L ) ) )   GOTO  750 

1 55  NAT  I N=WAT I N- ( ( WATRES ( I F I L ) *AWCNEW ( I  PAS , I RTYPE ,  I  F I  L )  )  - 

156  1  SOWAT < I  PAS, I RTYPE , I FI L) ) 

1 57  SOWAT ( I  PAS , I RTYPE , I F I L ) -WATRES ( I F I L ) *AWCNEW ( I  PAS , I R i YPE , I F I L ) 

158  GOTO  80  0 
153  750  CONTINUE 

160  SOWAT  (  I  PAS , I RTYPE , I F I L ) =SOWAT ( I  PAS , I RTYPE , I F I L )  +  WAT I N 

161  GOTO  850 

162  800  CONTINUE  (\  (]  f>  7  ^   


0 


[j 


0 

I 

0 
0 

0 

D 
0 

0 
0 


let       4    11:26    1384       sour ce . code    Page    12 


163  S5G  CONTINUE 

164  C 

165  C  DEPTH  OF  FROZEN  SOIL 

166  C 

167  FREZDPC I  FAS , I RTYPE) =FREZDF( I  PAS , I RTYPE)  + 

168  1      0.1  *  (32-DMTEMP) 

169  IF  C FREZDF < I  PAS , I RTYPE )  . LT . 0  .  )  FREZDP ( I  PAS , I RTYPE ) =0 . 

170  I F  ( FREZDP ( I  PAS , I RTYPE)  . GT . 30 . )  FREZDPC I  PAS , I RTYPE) =30 

171  C 

172  C  EROSION 

173  C 

174  RFAC  =  10.  *  CDPRECTC IPAS)**0 .0  374) 

175  CFAC  =  10**  (  -0  .  552-C  0  .  016*PCDUf-1)  ) 

1 76  SO I LOF ( I  PAS , I RTYPE ) =ER ( I RTYPE ) *RFAC*XLS( I RTYPE ) * 

177  1       CFAC*PMECH 

178  30  0  CONTINUE 
173  C 

ISO  C  END  OF  SOIL  TYPE  LOOP 

181  C 

182  10  00  CONTINUE 

183  C 

184  C  END  OF  PASTURE  LOOP 

185  C 

186  C 

187  C  FLOWS  IN  GANCE  MAHALA  AND  NORTH  FORK 

188  C 

183  FLGNC  =  0. 

130  FLMAH  =  0 . 

191  DO  1100   1=1 ,  NRT 

132  FLGNC  =  FLGNC  +  RUNOFFC 3 , I )*ARGNuC I ) 

133  FLMAH  =  FLMAH  +  RUNOFFC 3 , I )*ARMAHC I ) 

134  1100  CONTINUE 

135  FLMAH  =  FLMAH  *  FLFUDG  *  CUFS 

136  OS  =  C100  0./OCN)  -  10. 

137  ORUN  =  ( CDPRECTC 3) -0 .2*0S)**2)  /  CDPRECTC3)  +  0.8*OS) 

138  IF  CDPRECTC 3) . LT . C 0 -2*0S) )  ORUN  =  0. 

139  FLGNC=C FLGNC  +  ORUN*OAR)  *  FLFUDG  *  CUFS 
20  0  FTOTC1)  =  FTOTCD  +  FLGNC*AFCON 

201  FT0TC2)  =  FT0TC2)  +  FLMAH*AFCON 

20  2  FT0TC3)  =  FT0TC3)  +  NFFLONCNSTEP+1 )  *  AFCON 

203  FWTOTC 1 ) =FWTOTC 1 )+FLGNC*AFCON 

204  FWTOT  C  2 ) =FWTOT  C  2 ) +FLMAH*AFCON 

20  5  FWT0T(3)=FWT0TC3)+NFFL0NCNSTEP+1)*AFC0N 

206  C 

207  C  WATER  QUALITY 

208  C 

203  C  UPLAND  EROSION  REACHING  STREAMS 

210  UEGNC  =  0. 

211  UEMAH  =  0 . 

212  DO  120  0   1=1,  NRT 

213  UEGNC  =  UEGNC  4-  SOI  LOFC  3  ,  I  )*ARGNCC  I  ) 

214  UEMAH  =  UEMAH  +  SOI LQF<9, I )*ARMAH< I) 

215  120  0  CONTINUE 

216  UEGNC  =  UEGNC*SEDDEL 

217  UEMAH  =  UEMAH*SEDDEL 

217.2  TSEDG    =    TSEDG    +    UEGNC  H 0 H  ^  i j 
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217. 

218 

213 

220 

221 

222 

i.'2o 

224 

■—:  ■- :  ET 

226 
ciii!  7 
22b 
223 
230 
231 
232 
d  3  3 
234 

•-i  •— :  ET 

d'  O  r"' 

233 
233 
240 
241 
242 
243 
244 
245 
246 
247 
248 
243 
251 


254 
255 
256 


253 
260 
261 
262 
263 
264 
265 
266 
267 
268 
263 
270 
271 

C.7  c. 


TSEDM  =  TSEDM  +  UEMAH 
C 
C  BANK  EROSION  INTO  STREAMS 

BETOT  =  CAT FA C  *  CU(3)  *  RM(3) 
C 
C  WATER  QUALITY 

WAT  0  A  L  =  (  B  ET  0T+  U  E  GN  C+  U  EM  A  H )  / 
1  (FLMAH  +  FL6NC  +  0.0  01) 


C 

c 
c 


BANK  STABILITY 

WTWOOD  =  0 . 
DO  1400   1=1,  3 

IF  ( ( I T I ME . NE . 1 ) )  PCDUM1 =PC (3,1,1) 
IF  ( ( I T I ME . NE . 1) )  PCDUM2=PC (3,2,1) 
IF  ( ( I T I ME . EQ . 1 ) )  PCDUM1 =PCMAX (1,1) 
IF  ( ( I T I ME . EQ . 1 ) )  PCDUM2=PCMAX (2,1) 
J=l  +  I 

WTW  O  O  D  =WT  W  0  0  D  +  J*  (  (  P  C  D  UM 1 + P  C  D  UM  2 )  /  2  .  ) 
140  0  CONTINUE 

WTWOOD  =  WTWOOD  /  4 

BSTAB  =  BSTB1 * WTWOOD  -  BSTB2*CU<9) 

IF  ( BSTAB. GT. 100 . )  BSTAB  =  100. 


MAKE  UP  WEEKLY/BY  TIME  STEP  TOTALS 


C 
C 
C 
C  TEMPERATURE 

DO  1500   1=1 ,  NPAS 

IF  ( ( I .EQ.8) .OR. (I .EQ.3) )  GOTO  1500 

TMAX(I)  =  TMAX(I)  +  TMAXL(JDAY) 
150  0  CONTINUE 

TMAX(S)  =  TMAX(S)  +  TMAXL(JDAY)  *  TMULT 

TMAX(S)  =  TMAX(S)  +  TMAXL(JDAY)  *  TMULT 
C 
C  SOIL  MOISTURE 

DO  170  0   1=1,  NPAS 

DO  1650   J=l,-  NRT 

DO  160  0   K=l  ,  3 

SO( I , J,K)  =  SOFDG*100 .*(SOWAT(I ,J,K)/( (AWCNEW( I , J,K)+0 .0001) 
1  *WATRES(K)))  +  SO(I,J,K) 

160  0  CONTINUE 
1650  CONTINUE 
1700  CONTINUE 
C 
C  WATER  QUALITY 

WQUAL  =  WQUAL  +  WATQAL 
C 

C             END  OF  DAY  LOOP 
c  . 

C 

50  0  0  CONTINUE 
C 

C        CALCULATE  AVERAGE  VALUES 

C  TEMPERATURE  „,  r,  ,,  -,  -, 

DO  600  0   1=1,  NPAS  «J!JuJ» 

-t.i.w,-   T  ',   —   TMAVf  T  'i    /   Kin, 
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0 
0 

0 

0 

D 
0 

G 
I 
I 


274 

60  0  0  CONTINUE 

£l7D 

C 

27b 

C  SOIL  MOISTURE 

c  i    '' 

DO  70  0  0   I =1 ,  NPAS 

278 

DO  630  0   J=l ,  NRT 

27£1 

DO  6S0  0   K=l  ,  3 

230 

SO ( I , J , K )  =  SO ( I , J , K )  / 

NDAYS 

2S1 

680  0  CONTINUE 

2S2 

630  0  CONTINUE 

'd  'd  3 

700  0  CONTINUE 

283 

C 

230 

C  WATER  QUALITY 

231 

WQUAL  =  NQUAL  /  NDAYS 

232 

C 

233 

C 

234 

RETURN 

255 

END 

CATTLE         AND         ECONOMICS         SUBMODEL 

This  submodel    is    described    in    Section    5    of    the    report. 

1  SUBROUTINE  VIANDE 

2  %  CONTINUE  WITH  SAMAL.C  RETURN 

3  C 

4  I F ( I WEEK . LT . I STR ) RETURN 

5  IFdWEEK.GT.  ISTR.AND.IWEEK.LT.  I  ENR)GU  TO  100 
£  IF(IWEEK.EQ.IENR)  GO  TO  1000 

7  I F (I WEEK . GT . I ENR ) RETURN 

8  C 

3      C .SHIFT  PASTURES 

10  C 

11  DO  30  1=1,3 

12  IF(IPCNTU)  .EQ.IPYN(D)  GO  TO  35 

1 3  I PCNT ( I ) = I PCNT ( I ) +1 

14  GO  TO  30 

15  35    CONTINUE 

16  IGR(I )=IPTO( IGR(I) ) 

17  IPCNT(I)=1 
13  30    CONTINUE 
13  IUSE=0 

20  NPAST=0 

21  TGWA=0 . 

22  TGWAO=TGWA 

23  IGRT=0 

24  C 

25  C SPRING  CALLING,  ETC. 

26  C 
C  OW  ( 1 )  =  C  OW  ( 1)  +  C  OW  B+  C  OW  (  2  ) 
C:  OWN  ( 1  )  =  (  C  OWN  ( 1 )  *  C  OW  ( 1  )  +  C  OWW  (  2 )  *  C  0 W  (  2  )  )  /  (  C  OW  ( 1  )  +  C  OW  (  2 )  ) 


i 


cc 


23  CuW(3)=YRB 

30  COWW ('  3 )  =YEARW  „  n  r.  -y  {) 

UUU  J  O 
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31  CGW(2)=CALF+CLFB 

32  C0WW(2)=CALFW 

33  TCOW=CGW ( 1 ) +COW ( 2 ) +CON ( 3 ) 

34  100   CONTINUE 

35  C 

36  C... SELECT  THE  PASTURE,  IF  NECESSARY 

37  C 

38  IF( IGRT.LT.l)  GO  TO  105 
32  IFCNPAST.EQ.7)  GO  TO  60  3 

40  I F ( NPAST . EQ . 2 )  RETURN 

41  GO  TO  10  7 

42  105   CONTINUE 

43  NPAST=NPAST+1 

44  IFCNPAST.EQ.7)  GO  TO  603 

45  IFCNPAST.EQ.9)  RETURN 

46  DO  103  1=1 ,NFMAX 

47  IF(IGR( I ) .NE. NPAST)  GO  TO  103 

48  IUSE=I 

49  GO  TO  10  4 

50  103   CONTINUE 

51  104   CONTINUE 

52  107   CONTINUE 

53  IF(IOC( I  USE) .NE.l)  GO  TO  603 

54  I  F  (  I M  UN  C  H  .  GT  .  1 )  CALL  EAT  UM 

55  I F ( I MUNCH . GT . 1 )  GO  TO  1 923 

56  C 

57  C. ... .CALCULATE  TOTAL  HERD  REQUIREMENTS"-  FOR  THE  WEEK 

58  C 

53  TTREQ=0 . 

60  DO  40  1=1,3 

SI  TREQ( I)=( (YII+SLI*COWW( I )* . 454)/ . 454)*DPN*C0W( I ) 

62  TTREQ=TTREQ+TREQ(I) 

63  40    CONTINUE 

64  DO  150  K=l,3 

65  TDN(K)=0. 

66  BEAT(K)=0. 

67  150   CONTINUE 

68  C 

63  C..... CATTLE  DAYS  ON  RIPARIAN  HABITAT 

70  C 

71  DO  151  K=1,NPAS 

72  CU(K)=0.0 

73  151   CONTINUE 

74  CU( I  USE)  = 

75  >(C0W(1)4-C0W(2)+CQH(3))*DPW/(RA(  I  USE  ,  3  ,  8)4- .  0001 ) 

76  C 

77  C...  START  MUNCHING  ON  FORAGE,  FIRST  OVER  RANGE  TYPES 
7S  C 

79  DO  50  0  IRR=1,NRT 

80  IR=NRT-IRR+1 

81  IF(RA(IUSE,IR,S)  .LE.1.0E-5)  GO  TO  500 
32  C 

83      C LOOP  OVER  SPECIES 

84  c                                   an*  to 

85  DO    110    K=l  ,7  iJUUJ/ 

86  CE( I  USE, IR,K)=0 . 


0 

I 

Q 
0 
G 
0 
I 
i 
II 
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87  110   CONTINUE 

88  DO  40  0  IP=1 ,7 
83  IPP=IPREF(IF) 
30  C 

|    91      C LOOP  OVER  WATER  NEARNESS  TYPES 

32  C 

33  DO  30  0  IW=9,IWM 

94  I F ( I W . EQ . 3 . AND . I  CRT . EQ . 0 ) TGWA=TGWA+RA ( I  USE , I R , 8 ) *RA ( I  USE , I R , ! 

W) 

35  IF( IN.EQ.3) 

96  >X=RA ( I  USE , I R , 8 ) *RA ( I  USE , I R , I PF ) *RA ( I  USE , I R , I N ) 

37  IF(IW.EO.IO) 

93  >      x=RA ( I  USE , 1 R , I PP)*RA( 1  USE , I R , 8 ) * ( 1 . -RA ( I  USE , I R , I N-l ) ) 

99  XTEST=-1*TTREQ/(X+0.Q0001) 

10  0  IF  (XTEST.GT.174.)  XTEST=174. 

101  IF  (XTEST.LT.-180 .)  XTEST=-180 . 

102  XX =X* ( 1 . 0 - EX P ( XT  E ST ) ) 

103  CE(IUSE,IR,IPP)= 

104  >  CE C I  USE , I R , I PP ) +XX 

105  C 

106  C CALCULATE  INGESTION  OF  NUTRIENTS  BY  CATTLE 

107  C 

108  DO  200  K=l,3 

103  XXX=XX*TREQ(K)/(TTREQ*DPW*COW< K)+ . 00001 ) 

110  TREQ ( K ) =AMAX1 ( TREQ ( K ) -XXX ,0.0) 

111  C  WRITE (6,610)  X, XX, XXX 

112  C610   FORMAT ( IX, 3F1 5.6) 

113  T  DN C  K ) =T  DN ( K )  +  ( S  LT  DN ( IFF) * ( U  P ( I P  P ) / 100.) +YT  DN ( I P  P ) ) *XXX 

114  BEAT(K)=BEAT(K)+XXX 

115  20  0   CONTINUE 

116  TTREQ=TTREQ-BEAT(1)-BEAT(2)-6EAT(3) 

117  IF(TTREQ.LT. 1 .0E-5)GO  TO  401 

118  300   CONTINUE 

119  40  0   CONTINUE 

120  401   CONTINUE 

121  DO  450  K=l,7 

122  CE( I  USE , I R , K) =CE( I  USE , I R , K)/( RA( I  USE , I R , 8)  +  . 0001 ) 

123  450   CONTINUE 

124  IF(TTREQ.LT.1.0E-5)GO  TO  501 

125  50  0   CONTINUE 

126  501   CONTINUE 

127  C 

128  C CALCULATE  GROWTH 

123  C 

130  1333  CONTINUE 

130.2  IF(IGROW.GT.O)  CALL  OBESE 

130.4  IF(IGROW.GT.O)  GO  TO  603 

131  DO  560  K=l ,3 

132  DO  550  1=1,7 

133  I F ( CONN ( K ) . GE . WGAR ( I , 1 ) *2 . 2 . AND . COWW ( K ) . LT . WGAR ( I +1 , 1 ) *2 . 2 ) 

134  >  GO  TO  555 

135  550   CONTINUE 

136  555   CONTINUE 

137  DO  558  J  =  2,6  .                                   ')  (j  U  4  U 

138  XWG(J-1)= 

139  >  (WGARC 1  +  1 ,J)-WGAR(I  , J) )/( WGARC 1+1 , 1 ) -WGAR ( I  ,  1 )  )* 
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140  >  ( ( COWW ( K ) /2 . 2 ) -WGAR (I , 1 ) ) +W6AR (I , J ) 

141  558   CONTINUE 

1 42  WGA I N (  K )  =SLP ( TDN (  K ) /2  .  2  , XNG  ,  YNG  , 5)*2 .  2-A-DPN 

143  C  OWN ( K ) =  C  DWW ( K ) +W GA I N ( K ) 

144  560   CONTINUE 

145  C 

146  C. ... .DECISION  RULES  FOR  MOVING  TO  NEW  PASTURES 

147  C 

148  603   CONTINUE 

149  IGRT=IGRT+1 

150  I F( I GRT . GE . I PW(NPAST) ) I GRT=0 

151  RETURN 

152  C 

153  C.....THIS  IS  THE  FALL  PERIOD 

154  C 

155  1000  CONTINUE 

156  CALFS=C0WW(2)*PCS*C0W(2) 

157  C0W(2)=C0W(2)*(1 .-PCS) 

158  COWS* 

153  >  C OWN ( 1 ) *AMAX 1 ( 0 . 0 , ( C  ON ( 1 ) + C  0 W ( 2 ) ) - C OWMAX ) 

160  COW ( 1 ) =COW ( 1 ) -COWS/ ( COWW ( 1 ) + . 0  0  0  0 1 ) 

1 61  YEARS=COW ( 3 ) *COWW C 3 ) 

162  COW(3)=0. 

163  P  B  S  L  D = C  A  L  F  S+  C  OW  S+ Y  EA  R  S 

164  BEEF=PPCOW*CQWS+FPCLF*CALFS+PPYR*YEARS 

165  C 

166  C. ... .CALCULATE  CALMING  RATE 

167  C 

168  CALF" 

169  >  SLPCCOWW(l) ,XCLF,YCLFS6)*C0W(1) 

170  PCALF=(CALF/(COWd)+.00001))*100. 

171  v'  IF(ISTOP.GT.O)  RETURN 

172  C 

173  C. . . . .HAY  STUFF 

174  C 

175  THAYG=0.0 

176  DO  1010  1=1,3 

177  "  THAYG=THAYG+HAYP( I)*2000. 

178  1010  CONTINUE 

179  THAYR=0.0 

150  DO  1020  1=1,2 

151  HAYRC I )=DPWN*( (YI I+SLI*COWW< I )* . 454)/ . 454)*C0W< I ) 

1 82  HAYR ( I ) =HAYR ( I ) +HPP ( I ) *AMAX1 ( 0 . 0 , CM I N ( I ) -COWW ( I ) ) *COW ( I ) 
1 S  3  T  HAY  R  =T  HAY  R+  HAY  R  (  I  ) 

184  1020  CONTINUE 

185  HAY  B  =AM I N 1 ( T  HAY  R -T HAY  G , BUYHAY* 20  0  0 . ) * I AHB 

186  HAYS=AMAX1 ( THAYG-THAYR , 0 . 0 ) * I AHS 

187  T  HAY  =T  HAY  G+  HAY  B 

188  HAY  B = H A Y  B/ 200  0. 
188.2  I F ( HAYB . LT . 0 . ) HAYB=0 . 

183  H A Y  S = HAY  S/ 2000. 

130  C 

131  C.....GROW  COWS  OVER  WINTER 


132      C 

133 

194 


133  DO  10  30  I  =1  ,  2  ..'•  J  C  4  ' 

CGWWC I )=COWW( I )+ 


Bet 
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0 
Q 

0 
D 
D 

0 

E 

0 

1 
I 
I 
I 
I 


135 

196 
137 
198 
199 

20  0 
201 
202 
203 
20  4 
20  5 
20  6 
207 


3 

4 

5 

b 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

23 
29 
30 
31 
32 

34 
35 
36 
37 


40 


>SLP( HAYR( I )*THAYV  ( 1 . QE-5+(DPWN*C0W< I )*THAYR) )  ,XWNG , YNNG , 3)*DP 
WN 
10  30  CONTINUE 

IFCIECON.LT.DRETURN 
COSTS= 

>  T  HAY  G*  D  G  H+  HAYB*  D  B  H+ F  C+M  C*TC  QW+ 

>  •    DAN*  (  T  GN  A  -T  GN  A  0 ) 
R  EU = B  E  E  F+  D  S  H*  HAY  S 
RNP=C0STS-(FC/2. ) 
GVO=Xf-1l*RNP 
CVA=XM2*GV0 
RNA=XM3*RNP 

RETURN 
END 


*CONT 
C 


40 


SUBROUTINE  EATUM 
INUE  WITH  SAUAL.C  RETURN 

TTREQ=0 . 

DO  40  1=1,3 

TREQC I )=( (YI I+SLI*COWW< I )*.454)/.454>*DFW*C0W< I ) 

HCI )=1 .0/( (TREQC I )+. 0000001 )/(CON( I )+. 0000001) ) 

TTREQ=TTREO+TREG( I ) 

CONTINUE 


C 
C 

c 


90 


OT 


95 
100 


.BOTTOM  OF  RESPONSE 

DO  100  1=1,3 

TDNCI )=0 . 

FBOT=l . 

DO  90  IF=1,7 

DO  90  I R=l , 3 

DO  90  IW=1 ,2 

FBOT  =  FBOT+ALPH( I )*H< I )*PREFF( I F , I R , IW)*RA( I  USE , I R , I F) 

TDIE(IR,IR,IH)=0. 

TTAT(I ,IF,iR,IN)=0. 

CONTINUE 

DO  95  IF=1 ,7 

DO  95  IR=1,3 

DO  95  IW=1 ,2 
X=ALPHCI)*FREFF(IF,IR,IW)*(CQWCI)/<RA(IUSE,IR,8)+.000001))/FB 

IF(X.GT.100 . )X=100 . 

TDIE(IF,IR,IW)=TDIE(IF,IR,IW)+X 

TTATCI ,IF,IR,IW)=X 

CONTINUE 

CONTINUE 


C 

C. 
C 


.CALCULATE  CONSUMPTION,  LOSS 

DO  20  0  IF=1,7 

XX=SLTDN ( I F ) * ( UP ( I F ) /1 0  0 . ) +YTDN (IF) 

DO  20  0  IR=1 ,3 

DO  20  0  IW=1 ,2 

X=l  .  0-EXP(-TDIE(IF, I R , IN) ) 

DO  150  1=1 ,3 


004 


4b 

46  . 

I 

46. 

o 

4b  . 

3 

47 

48 
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41  Y=(TTAT(I  , IF, IR, IW)/TDIE< IF,IR, IW)+ .000001 )*X*RA( I  USE, IR,8)*X 
X 

42  >*RA( I  USE, IR,8+IW)*RA(IUSE, IR, IF)/(COW( I )+. 000001) 

43  T  DN ( I ) =T DN ( I ) +Y/ ( D PW+ .0000001) 

44  150   CONTINUE 

45  CE< I  USE , I R , I F) =CE( I  USE , I R , I F)*( 1 . -( (X*RA< I  USE , I R ,8)*RA< I  USE , I 
R , S+IW) )/RA( I  USE , IR , 8)  +  . 00001) ) 

45 . 2  CE< I  USE , I R , I F) =CE( I  USE , I R , I F)*( 1 . -< (X*RA( I  USE , I R , 8)*RA< I  USE , I 

R,8+IW))/RA( I  USE, IR, 8)+. 00001)) 
20  0   CONTINUE 

DO  399  IR=1 ,NRT 

DO  399  IF=1  ,7 

CE( I  USE, IR, IF)=RAC I  USE, IR, IF)*(1 .-CE( I  USE, IR, IF) ) 

RETURN 

END 

1  SUBROUTINE  OBESE 

2  ^CONTINUE  NITH  SAMAL . C  RETURN 

3  DO  50  1=1,3 

4  WGAIN(I)=0. 

5  50    CONTINUE 

6  DO  100  1=1,3 
6.5  DO  85  IR=1 ,NRT 

7  DO  80  I F=l  ,  7 

8  TT  = ( 3 LT  DN ( I F ) ) *  < U P ( I F ) / 10  0. ) + YT DN (IF) *  100. 
3  EM=TT*EMP 

1 0  F=EXP C FP1 -FP2*EM ) 

11  EMN=EMNP/F 

12  GEN=GENP1-6ENP2*F 

1 3  REMN=PM1*C0WW ( I ) **PM2 

14  RIN=REMN/EMN 

15  SURP=(CE( IUSE, IR, IF)*RA( IUSE, IR,8)/(C0W( I )*DPW+. 0001) )-RIN 

16  DN  E  G = S  U R P* G EN 

17  GO  TO  (1 ,1  ,2) , I 

18  1     CONTINUE 

19  X=Gia)*SQRT(G2(l)+G3(l)*<DNEG/C0WW( I )**PM2) )-G4(l) 

20  GO  TO  10 

21  2     CONTINUE 

22  X=G1  (  2 )  *SQRT  (  G2  (  2  )  +G3  (  2 )  *  (  DNEG/  (  CONN  (  I  )  **F'M2 )  )  )  -G4  (  2  ) 

23  10    CONTINUE 

24  WGAIN(I)=WGAIN( I )+X 

25  SO    CONTINUE 
25.5     85    CONTINUE 

26  100   CONTINUE 

27  DO  20  0  1=1,3 

28  CONN ( I ) =COWW ( I ) +NGA I N ( I ) *DPW 

29  20  0   CONTINUE 

30  RETURN 

31  END 

WILDLIFE    SUBMODEL 

."  0043 

This  submodel  is    described  in  Section  6  of  the  report. 
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1  SUBROUTINE  WILD 

1.5  * C ONT I N  U E  W I T H  S  A  V  A L . C  RET U RN 

2  C 

3  C   SECOND  VERSION  OF  SAVAL  WILDLIFE  MODEL  CREATED  AFTER  WORKSHOP 

4  C   JAN  3 ,1982 

5  C 

6  M- I MONTH 

6.1  I  Ft!  IWEEK.NE.DGO  TO  3 

6 .2  DO  2  1=1,3 

6.3  IC(I)=0 

6.4  2     CONTINUE 

6.5  3     CONTINUE 

7  1     I F ( I WEEK . NE . 1) GOTO  10 

8  C 

3  C   *  *  *  *  *  *  *  *    DEER     MODEL   *  *  *  *  *  *  *  *  *  * 

10  C 

11  C   OVER  WINTER  MORTALITY 

12  C 

13  H E A L A = D  EN  E R A/ ( AN E  E  D* 3  8 ) 

14  HEALF=DENERF/(FNEED-A-3S.  ) 

1 5  DMA  L  E = DMA  L  E* S  L  P ( H  E A  L A , XM A  LED, YM  A L  E  D  ,  2 ) 
IS  DFEM=DFEM*SLP<  HEALA ,XFEM , YFEM ,2) 

1 7  FAWN =FAWN* S L  P ( H E A L  F , X FN Y , Y FN Y , 2 ) 

18  C 

IS  C   ADD  FAWNS  INTO  ADULT  POPULATION 

20  C 

21  DM  A  L  E = DM  A  L  E+  FAWN*  F  P  E  R  C 

22  DFEM-DFEM+FAWN* ( 1 . -FPERC ) 

23  TDEER=DMALE+DFEM 

24  FAWN=0 

25  DENERA=0 

26  DENERF=0 

27  DSEAS=0 

28  DO  5  1=1,3 
28.2  ITRIG(I)=0 
23  IC(I;=0 

30  5     CONTINUE 

31  10    CONTINUE 

32  C 

33  C   SEE  IF  FORES  ARE  GROWING  IN  RESPECTIVE  PASTURES 

34  C   LOWER  SHEEP 

36  DO  20  J=l ,3 

37  IF(GPLANT(6, J, 4) . GT . 0 ) ITRI G( 1) =1 

38  20    CONTINUE 

39  IFCITRIG(l) .GT.0)IC(1)=IC(1)+1 

40  C 

41  C   IF  MORE  THAN  8  WEEKS  HAVE  PASSED  AND  NO  GROWTH  EVIDENT  THEN  MOVE 

42  C   DEER  ONTO  LOWLANS 

43  C 

44  IF( IWEEK.GT.S.AND.ITRIG(l) . EQ . 0 ) I C( 1 ) =2 

45  IF< IC(1) . LT.2)G0T0  105 

46  I F ( DSEAS . GE . 2 ) GOTO  50 

47  C 

48  C   DEER  ARE  MOVING  INTO  LOWER  PASTURE  . 

49  DSEAS=1  >  j  U  U  4  4 


w'JJ 


w'O 
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50  D  E  E  R = F  R  A  C*  <  D  F  EM+  DM  A  L  E  > 

51  FAWN=U 

52  IF< IC(1) .EQ.2)G0T0  10  0 
50    CONTINUE 

4  C 

5  C   SEE  IF  FORBS  ARE  GROWING  IN  UPPER  SHEEP 

7  DO  60  J=l,3 

I F ( GPLANT ( 7 , J , 4 ) . GT . 0 . . OR . GPL ANT ( 5  ,  J ,  4 ) . GT . 0 . ) I TR I G ( 2 ) -1 

53  60    CONTINUE 

60  IFC ITRIGC2) . GT . 0  ) I C( 2) =1 C( 2)4-1 

61  IF( IC(2) .LT.2)G0T0  100 

62  IF(DSEAS.GE.3)G0T0  SO 

63  C 

64  C   DEER  ARE  MOVING  INTO  MIDDLE  ELEU .  PASTURES 

65  DbErtb=2 

SS  DEER=DFEM+DMALE 

67  65    CONTINUE 

68  IF(IC(2) .EQ.2JG0T0  100 

69  C 

70  C   SEE  IF  FORBS  ARE  GROWING  IN  INDEPENDANCE  PASTURES 

71  SO    CONTINUE 

73  DO  85  J=l , 3 

74  IF(GPLANT(7? J,4) . GT . 0 . ) ITRI G( 3)=1 
85    CONTINUE 

IF(ITRIG(3)  .GT.O)  I  C(  3)  =1  C(  3)-fl 

IFCICC3) .LT.2)G0T0  100 

DSEAS=3 

79  DEER=DFEM+DMALE 

SO  100   CONTINUE 

81  C 

82  C   ALLOW  DEER  TO  EAT  FORAGE 

83  GOTO  499 

84  105   CONTINUE 

85  C 

86  C   FAWNS  ARE  BORN  IN  WEEK  10  (JUNE  1) 

S7  I F( I WEEK . EQ .10) FAWN=SLP( HEALA ,XFAW , YFAW , 2)*DFEM 

88  C 

89  C   HUNT  MALE  DEER 

30  HDAYS= HF  RA  C*  P  E RM I T* DAY  H/3 7 

91  XKI LL=SLP ( DMALE , HHX , HHY ,3) 

32  TK I LL=HDAYS*XK I LL 

93  DMALE=DMALE-TKILL 

94  IF(DMALE.LT.l . )DMALE=1 . 
98  GO   TO    525 
93  C 


76 


o 


100 


C  DEER  INGESTION  MODEL 

101  C 

102  439   CONTINUE 

103  IFLAG=DSEAS 

104  TOTIN-0. 

10  5  TOTUEG=0 . 

106  C   AMOUNT  DEER  WANT  TO  EAT 

10  7  T  D  EM  D A = D  E  E  R* D I N GW ( 2 ) 

108  TDEMDF  =  FAWN*DINGW(1)  00045 

109  T  GT I N =T  D EM  D A+T  D  EM D F 

110  IX=NPERP( I  FLAG) 


let   4  11:26  1984   source. code  Page  22 

C   LOOP  OVER  PASTURES 

DO  50  0  1=1  ,  IX 

I I-NPA8DC I  FLAG, I ) 

JX=NPERT(IFLAG) 
C   LOOP  OVER  RANGE  TYPES 

DO  50  0  J=1,JX 

JJ=NTYPE( I  FLAG, J) 
C   LOOP  OVER  7  PLANT  GROUPS 

DO  50  0  K.=l  ,7 
C 
C   SUBTRACT  VEG  EATEN  BY  CATTLE  THIS  WEEK  IF  ANY 

RRA( I  I  , JJ , K) =RA( I  I  ,  JJ , K) -CE( I  I  , JJ , K) 

IFCRRAC I  I , JJ.K) .LT.0.)RRA( I  I , JJ,K)=0 . 
C   OVREALL  WEIGHTED  FOOD  AMOUNT  AND  PERCENT  OF  DIET 

AMT=RRA ( I  I , J J , K ) *RA <  I  I  ,  J  J  ,  8 ) 
126      C   TOTAL  BIGMASS  WEIGHTED  BY  PREF  IN  AREA  USED  DURING  CURRENi  SEA30 


111 

112 

113 

114 

115 

116 

117 

118 

113 

120 

121 

1 22 

123 

124 

X  w  o 

N 

127      C 


128 

129 

130 

131 

132 

133 

134 

135 

136 

1  o7 

1 38 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

T  OTV  E  G =T  OT V  E  G+  P  R  E  F ( K , M ) *AMT 
50  0   CONTINUE 

DO  505  1=1, IX 

II=NPASD(IFLAG,I ) 

DO  505  J=1,JX 

JJ=NTYPE(IFLAG, J) 

DO  50  5  K=l,7 
C 
C   DETERMINE  FROF'-'N  OF  EACH  VEG  TYPE 

AMT=RRA ( I  I  , J J , K ) *RA ( I  I , J J , 8 ) 

FERINCI I , JJ,K)=0. 

I F ( TOTVEG . ST . 0 . ) PER I N ( I  I , J J , K ) =PREF ( K , M ) *AMT/ I OTVEG 
C 
C   AMOUNT  OF  VEGETATION  EATEN 

WECII ,JJ,K)=0 
C 
C   TOTAL  BIOMASS  OF  EACH  VEG.  TYPE  EATEN  BY  DEER 

IF(RA(II  ,JJ,8)  .GT.0)WE(II  , J J , K) =PERIN( I  I  , J J  ,  K)*TOTIN/RA(I I  ,JJ 
,8) 
u 
C   CHECK  THAT  DEER  NOT  EATING  MORE  THAN  AVAILABLE 

IFCWEC I  I , JJ,K) .GT.RRAC I  I  , JJ,  K) )WE< I  I  , JJ,K)=RRA( I  I  , JJ,K) 
510   CONTINUE 
C 
C   ADD  WHAT  DEER  GET  INTO  ENERGY  POOL(  PER  INDIVIDUAL) 

I FCTOTIN . LE . 0 . . OR . DEER . LE . 0 . ) GOTO  505 

DENERA=DENERA+WE ( I  I , J J , K ) *RA ( I  I  , J J , 8 ) * ( TDEMDA/TOT I N ) *CAL ( K , M ) 

/DEER 

I FCTOTIN. LE.O . . OR . FAWN . LE . 0 . )GOTU  505 

DENERF=DENERF+WE(I I  , J J , K)*RA( I  I  , J J , S)*(TDEMDF/TOTIN)*CAL( K ,M) 
/FAWN 
50  5   CONTINUE 

GOTO  105 
C 
C  *  +    *       SAGE    G  R  0  U  8  E   M  0  D  E  L   *  *  *  *  *  *  * 

C 

C   CALCULATE  VARIOUS  AVERAGES  ";  0  0  4  6 

C   AVERAGE  FORES  AND  GRASSES   APRIL  1   JUNE  15  iJ  "-  J 
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1(53  525   IF(  IWEEK.GT.DGOTO  550 

164  DO  540  1=1 ,5 

165  FA1J15(I)=0 

166  GA1J15(I)=0 

167  FJ1S1CI)=0 
16S  GJ131(I)=0 

169  CM1M3K  I)=0 

170  FM15S3(I)=0 

171  CM15S3C I)=0 

172  GM15S3(I)=0 

173  540    CONTINUE 

174  550   CONTINUE 

175  I  FCINEEK.GT.il) GOTO  630 

176  I  FdWEEK.LT.  3)  GOTO  630 

178  C   LOOP  OVER  STRUTTING  GROUNDS 

179  DO  620  I =1,1  SAGE 

180  C   SET  UP  FOR  CORRESSPONDING  PASTURES 
1S1  II=ISTRUTCI) 

182  GA=0 

183  C   CALCULATE  AVERAGE  FORBS  OVER  2  RANGE  SITES  -  LOAMY  AND  CLAY 
1S4  FA«(RACI I ,1 ,4)+RA(II ,2,4) )/2. 

185  C   CALCULATE  AVERAGE  DENSITY  OF  GRASSES  OVER  2  RANGE  SITES  AND 

186  C   3  PLANT  GROUPS 

187  DO  619  K«5,7 

188  GA=GA+  C  RA  C  I  I  ,  1  ,  K )  +RA  C  I  I  t  2  ,  K )  )  /6  . 

189  619   CONTINUE 

190  C   NOW  AVERAGE  OVER  TIME   -  9  WEEKS 

191  FA1J15C I )=FA1J15C I )+FA/9. 

192  GA1J15C  I  )=GA1J15C  I  )+GA/9. 

193  620   CONTINUE 

194  630   CONTINUE 

195  C 

196  C   AVERAGE  FORBS  AND  GRASSES   JUNE1   SEPT1 

197  I FCIWEEK.LT. 12) GOTO  640 

198  IFCIWEEK.GT.24)G0T0  640 

199  DO  635  1=1 „ I  SAGE 

200  II=ISTRUTCI) 

201  C   AVERAGE  FORBS  OVER  2  RANGE  SITES  IN  EACH  PASTURE 

202  FA=CRAC 11,1 , 4)+RAC II ,2,4) )/2. 

203  GA=0 

204  C  AVERAGE  GRASSES  OVER  2  RANGE  SITES  AND  3  PLANT  GROUPS 

205  DO  634  K=5,7 

206  GA=GA+C'RAC  11,1  ,K)+RA(II  ,2,K))/6. 
20  7  634   CONTINUE 

208  C   NOW  AVERAGE  OVER  THE  TIME  PERIOD 

209  FJ1S1CI)=FJ1S1CI )+FA/13. 

210  GJ1S1C  I  )=GJ1S1C  D+GA/13. 

211  635   CONTINUE 

212  640   CONTINUE 

213  IFCIWEEK.GT.il) GOTO  650 

214  C   AVERAGE  %  SHRUB  COVER   MAR  1   MAY  31 

215  DO  645  I =1,1  SAGE 

216  II=ISTRUTCI) 

217  C   AVERAGE  %    COVER  OVER  2  RANGE  SITES 

218  C   MARCH  15  TO  MAY  31  .,  ^  a  y» -> 

21 9  CO=  C  PC  C  I  I  ,  1  ,  2 )  +FC C  I  I  ,2,2)) /2  .  ,  ■•'  J  A  i 
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220  C   AVERAGE  OVER  TIME  PERIOD  -  11  WEEKS 

221  CM  1 M  3 1  (  I  )  =  CM  1 M  3 1  (  I  )  +  C  0/ 1 1  . 

222  645   CONTINUE 

223  650   CONTINUE 

224  C 

225  C   AVERAGE  PORBS , GRASSES  AND  %  COVER   MAR  15   SEPT  30 

226  I P ( I WEEK . GT . 28 ) GOTO  660 
DO  655  I =1,1  SAGE 


a 
a 

o 

0 
0 
0 


£-£.  • 


228  1 1=1  STRUT (I) 

223  C   AVERAGE  PORE  DENSITY  OVER  2  RANGE  SITES 

230  PO= ( RA ( I  I  ,1,4 ) +RA (II  ,2,4)) /2 . 

231  C   AVERAGE  COVER  OVER  2  RANGE  SITES 

232  CO=(PC( 11,1 , 2)+PC( II ,2,2) )/2 . 

233  GA=0 . 

234  C   AVERAGE  GRASSES  OVER  RANGE  SITES  AND  3  TYPES  OP  GRASS 

235  DO  654  K=5 , 7 

236  GA=GA+ ( RA ( I  I  , 1 , K ) +RA (II , 2 , K ) ) /6 . 

237  654   CONTINUE 

238  C   NOW  AVERAGE  OVER  TIME  PERIODS 

239  FM15S3( I )=PM15S3( I )+P0/28 . 

240  CM15S3( I )=CM15S3( I )+C0/28. 

241  GM1 5S3 ( I ) =GM1 533 ( I ) +GA/28  . 
241.2  655   CONTINUE 

242  660   CONTINUE 

243  IP( I WEEK. NE. 37) GOTO  1000 

244  C 

245  C   NOW  DETERMINE  POPULATION  DYNAMICS  OVER  LAbT  YEAR 

246  C   CHICK  NUMBERS   SEPT  31 

247  C   PIRST  NESTING  SUCCESS  BY  STRUTTING  GROUND 

248  DO  70  0  1=1,5 

249  T0T=FA1J15( I)+GA1J15(I ) 

250  TT=0 

251  IP(TOT.GT.O)TT=PA1J15(I)/TOT 

252  ANEST ( I ) =SLP ( TOT , ANX , ANY , 3 ) 

253  I F ( TT . LT . . 1) ANEST ( I ) =ANESMN 

254  YNEST ( I ) =SLP ( TOT , YNX , YNY , 3 ) 

255  I F ( TT . LT . . 1 ) YNEST ( I ) =YNESMN 

256  C 

257  C   NOW  CHICKS  PER  HEN 

258  T0T=FJ1S1(I)+GJ1S1(I ) 

259  CHEN( I ) =SLP(TOT , CHHX , CHY ,3) 

260  TT=0 

261  IF(TQT.GT.0)TT=FJ1S1( I )/TOT 

262  I F ( TT . LT . . 1 ) CHEN ( I ) =CHENMN 

263  C   CALCULATE  #  OF  CHICKS 

264  ■  CHICK? I )=( ANEST ( I )*AHEN( I )+YNEST( I )*YHEN( I ) )*CHEN( I ) 

265  700   CONTINUE 

266  C 

267  C   DETERMINE  HEN  SURVIVAL   MAR  15   SEPT  30 

268  C   OPTIMAL  ACCORDING  TO  PORBS  AND  GRASSES  <S    25%  COVER 

269  DO  705  1=1 ,5 

270  T0T=FM1 5S3 ( I ) +GM1 5S3 ( I ) 

271  SURVH(I  )=SLP(T0T,XHS'UR,YHSUR,3) 

"272  C   MODIFY  ACCORDING  TO  COVER                               nf!A/!n 

,!  U  H  O 


273  FAC=SLP(CM15S3( I ) ,XHFAC , YHFAC , 4) 

274  C   MALE   SURVIVAL 


UUUHo 
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£.(    w' 


0 
D 
I 


SURVM< I ) =SLP( CM1M31 ( I ) ,XMSUR , YMSUR , 5 ) 

276  70  5    CONTINUE 

277  C 

278  C   HUNT  THE  GROUSE  -  FALL  POPULATION 

279  DO  70  6  1*1,5 

280  AHEN ( I ) =AHEN ( I ) * ( 1 . -SHUNT ) 

281  AMALE  (  I  ■)  =AMALE  C  I  )*<  1 .  -SHUNT ) 

282  YHEN ( I ) =YHEN ( I ) * ( 1 . -SHUNT ) 

283  YMALE ( I ) =YMALE  C I ) * ( 1 . -SHUNT )  ■ 

284  706   CONTINUE 

2S5  C   CALCULATE  NUMBER  FOR  EACH  POPN  SEPT  30  AND  ADD  RECRUITS 

286  XTT=TDTHEN+TOTMAL                                                 Li 

288  TOTHEN=0 

2 8 3  TOTMAL=G                                                     n 

250  DO  710  1=1  ,5 

231  AHENFC I ) =( AHEN( I )+YHEN( I ) )*SURVH( I )*FAC 

232  AMALEF ( I ) = ( AMALE ( I ) +YMALE ( I ) ) *SURVM CD 

233  C  I 
294  C  CALCULATE  POP'N  AFTER  WINTER (MAR  15)  ■ 
235  AHENC I ) =AHENF( I )*WHENS 

296  AMALE ( I ) =AMALEF ( I ) *WMALES 

297  YHEN ( I ) =CH I CK ( I )*  C 1 . -CMALE ) *WHENS  ■ 
2  9  8  YMA  LE(I)=CHICK(I)* CM A  L  E* WM A L  E  S 

299  T  OT  H  EN  =T  OT  H  EN+ A  H  EN  <  I  )  +  Y  H  EN  ( I  ) 

300  TOTMAL=TOTMAL+AMALE ( I ) + YMALE ( I ) 

301  710   CONTINUE 
30  2  SAGEQ=TOTHEN+TOTMAL 
30  3  100  0  CONTINUE 
30  4  RETURN 
30  5  END 


049 


n  f\  C\  A  0 


1 

G 


i 


1 

D 
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3.2  DATA  FILES 

3   A    M   A    L     .     D 


This  is  the  main  data  file  called  by  the  initialization  subroutine 
UINIT.  Data  files  associated  with  each  subroutine  are  attached  to 
this  file  by  the  use  of  ^CONTINUE  statements.   The  data  files  are 
listed  accordinq  to  the  order  they  appear  in  SAMAL . D . 


1  S  GPLANT(ALL,ALL,ALL)=0 . 
1.02  ^CONTINUE  NITH  GROW.D  RETURN 

1.05  ^CONTINUE  WITH  CHOW . D  RETURN 

1.1  S  IBEEF=1 

1.2  S  IWAT=1 

1.3  3  IWILD=1 

1.4  3  IUEG=1 

2  ^CONTINUE  WITH  DAVE.DD  RETURN 

3  3  DPRECT(ALL)=0  . 

4  3  CU(ALL)=0. 

5  ^CONTINUE  WITH  EEEF.D  RETURN 

6  * CONTINUE  WITH  DAYE.D  RETURN 

7  ^CONTINUE  WITH  WILD.D2  RETURN 

8  * CONTINUE  WITH  WAT . D  RETURN 

3.01  SAME  FILE=-S 

3.02  MACRO  GO 
8.0  4  RUN  1  13 

8.06  SET  #NYSKIP=38 
8.08  RUN  20  476 

8.1 

8.12  MACRO  IRR 

8.14  3  HAY0P(ALL)=2 

8.16  3  ISTR=7 

8.18  S  IENR=37 

8.2 

8.22  MACRO  STOCKER 

8.24  S  HAYOP(ALL)=3 
8.26  S  IAHB=0 

5.25  S  COW(ALL)=500  300  1 
AT  33  'S  YRB=150/ 
AT  77  'S  YRB=30  0' 

8.34  AT  115  '3  YRB=450 ' 
8.36  AT  153  '3  YRB=60  0' 

5.35  AT  49  -S  SHKILL(3)«0.-6' 


8.32 


8.385  AT  43  'S  IMM*1' 
8.39  AT  50  "3  IVM=0' 
8.4    AT  50  'S  SHKILL<3)*O.0' 


•s 


42   AT  43  'S  <v'MANIP(5,2,3)=2203' 


8.44   AT  50  'S  VMANI P( 5 , 2 , 3) =0 ' 


e 


.46   AT  49  '  S  UMAN I  P  (  7  ,  2  ,  3 )  =3568 '  ■■.  c\  n  r,  f) 


8.48   AT  50  "S  UMAN I P( 7 , 2 , 3) =0 ' 
MACRO  FIRE 


JO 

JO 
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8.63 

8  IAHB=1 

3 .  68 

AT  30  •"  S 

UMANIP(3,1 ,2) =2000 

8  .  69 

AT  30  "S 

IVM«1' 

8 . 7 

AT  31  '8 

IUM=0'" 

8.73 

AT  30  'S 

SHKILL(3)=i .0" 

8 .  78 

AT  31  "S 

UMANIP(3,1 ,2)=0 .0" 

8.83      AT    31     ' S    SHK I LL ( 3 ) =0 . 0 


G    R    0   W 


1  S  IGROW=0 

2  S  EMP*. 036155 

3  8  EMNP=77. 

4  S  FP1=2.2577 

5  S  FP2=.2213 

6  3  GENP1=2.54 

7  S  GENP2-.0  314 

8  8  PM1=.0  77 
3  8  PM2=.75 

10  8  61 (ALL) =73. 099    33.526 

11  8  G2(ALL)=. 002773    .003133 

12  S  G3< ALL) =.02736    .0506 

13  S  G4(ALL)=3.8538    2.2146 

CHOW    .    D 


1  8   NPREF(ALL)=2    1 

2  S    RPREF(ALL)=1    2    3 

3  8    FPREF(ALL)=0 .1     .05    .05    .2    .3    =1 

4  8   ALPH(ALL)=126 

5  8    IMUNCH«2 

D   A   U  E    ,    D    D 


1 

8  IOC(ALL)=l  14*0 

2 

8  NRT*3 

3 

S  NPA3=3 

4 

8  RA(ALL,ALL,ALL)=0 

5 

8  CE(ALL,ALL,ALL)=0 . 

6 

8  WE (ALL, ALL, ALL )*0. 

7 

8  FTOT(ALL)=4000  5000 

8 

SAME  FILE=-S 

BEE 

F  .  D 

1 

8  COWW(ALL)=700  350  7 

1.1 

8  BUY HAY* 20 00. 

X  o  r— 

8  IUSE=0 

1.5 

S  IECON=l 

C 

S  ISTOP=0 

3 

S  CQWWMN=700 . 

800  0 


700 


005 


0 

o 


1U 

a 

11 

c 

12 

3 

13 

s 

14 

Q 

15 

s 

16 

b 

17 

s 

IS 

s 

13 

o 

20 

s 
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CE CALL, ALL, ALL )=0 

YI 1=4.1 

SLI=.0136 

CALF=350 . 

COW (ALL) =50  0  0  1 

IFYN(ALL)=7*1  2  2 

IGR(ALL)=2  2  S  4  3  1  3  6  5 

IPCNT(ALL)=7*1  2  2 

IPT0(ALL)=3  14965223 

IPN(ALL)=3  3  2  2  6  6  4  4  0 

IPREF(ALL)=5  6  7  2  14  3 

DPN  =  7 

!GRT=0 

IPAST(ALL)=1  2  3  6  4  5  7  3  8 

IWM=10 

YTDN(ALL)=.4 

SLTDN(ALL)=1.3 

21  S  YWG ( ALL ) =-1 .0  0.0  0.5  0.7  0.3 

22  3  FTHR=50 

23  S  XCLF(ALL)=0  60  0  SCO  350  10 DO  20  0  0 

24  3  YCLF(ALL)=0  0  .85  .85  .6  .6 

25  8  PCS=.8 

26  8  PAS=. 08 

27  S  CALFW=SG . 

28  8  HPC=.025 
23      8  DPWN=156 

30  3  IAHS=1 

31  3  IAHB=1 

32  8  XWNG(ALL)=G  22  10  0 

33  3  YWNG(ALL)=-1 .5  0  5 

34  8  WGAR(1,ALL)=0 

35  8  WGARC2,ALL)=100  0  .5  1.  1.5  1.8 

36  3  WGAR(3,ALL)=150  0  1.  1.5  2.5  2.7 

37  3  WGAR(4,ALL)=20  0  0  1.3  3.5  3.8  4.0 

38  8  NGAR(5,ALL)=250  0  2.3  3.3  4.  4.6 
33      8  NGAR(6,ALL)=300  0  2.6  4.5  4.7  5.2 

40  S  MGAR(7,ALL)=350  0  2.3  5.1  5.4  6.0 

41  S  WGAR(8,ALL)=333  0  4*1000 

42  8  TGWA=0 . 

43  3  TGWAO=0 . 

44  8  CONB=0 . 

45  S  YRB=1 . 

46  3  CLFB=0 . 

47  8  PPC0W=.4 

48  3  PPCLF=.6 
43      8  PPYR=.63 


1 


50  3  DGH=0.0i: 

51  S  DBH=46. 

52  8  FC=2.5E5 

53  8  UCC=60 . 

54  3  DAN =5 

55  8  DSH=60  . 

56  3  XM1=2.5 

57  3  Xf-12=.4 

58  8  XM3=.i5 


CMIN(ALL)=700.  70  0.  ()  Q  Q  5  2 


Oct 
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60 
61 
62 
63 
64 
65 
.66 
67 
68 
63 
70 
71 
72 
73 
74 
75 


DAVE 


S    RA(1  .  .  .  5,ALL,9)=1 . 
S    RAC1 . . .5,ALL,10)=0 
S    RAC6,ALL,9)=.5 
S    R  A  i'6i  ALL,  10)  =  .  5 


^LL,_ 
R A (6, ALL, 10)=. 5 
S    RA(7,ALL,3)=.65 
S    RA(7,ALL,10)=.35 
S    RAC8,ALL,9)=.S 
S    RA ( 8 5  ALL , 1 0 )  =  . 2 
S    RAO,  ALL  ,9)  =  .65 
S    RAO. ALL,  10)  =  .35 


S 

s 
s 


,hLL  , 
I STR=5 

t  r-K  ;rj  — -~t  c 

i  lisK-jj 

YEARW=7Q0 . 
NPMAX=9 
COWMAX=1QOG 
HPP(ALL)=5. 


1 

2 

3 

4 

5 

6 

7 

8 

3 

10 

11 

12 

13 

14 

15 

16 

17 

IS 

19 

20 

21 

22 

23 

24 

25 

26 

27 

29 
30 
31 
olc 
33 
34 
35 
36 


S 

s 
s 

c 

s 
s 


s 

s 
s 
s 
s 
s 
s 
s 

3 

s 
s 
s 

s 

s 
s 

o 

s 

C; 

s 


RQEVLB(ALL)*=3* 

ILUMF-G 

1^-1=0 SNITCH  FOR 

RANEW CALL, ALL, ALL )=0 . 
NEWRS CALL, ALL )=0 
GSOXM ( ALL , ALL ) =0 . 
FCMAX2 CALL, ALL )=0 . 
PCM AX (ALL, ALL ) =0 . 
PWTMAX CALL, ALL )=0. 
AWCNEW CALL, ALL, ALL )=0 . 
GS0XM(1,ALL)=16.  17.  21. 
GSOXM (2, ALL) =13.  IS.  24. 
GSOXM (3, ALL) =9.  12.  16. 
PCMAX2(1,ALL)=30.  45.  15 
PCMAX2(2,ALL)=35.  20.  30 
PCMAX2(3,ALL)=25.  5.  60. 
PCMAXC1 ,  ALL) =15.  30.  0. 
PCMAXC2,ALL)=30 .  35.  0. 
PCMAX ( 3 , ALL ) =60 .  25 .  0 
PWTMAX (1, ALL) =45. 
PWTMAX C 2, ALL) =65. 
PWTMAX ( 3 , ALL ) =70 . 


5  4*10  LB.  OF  WATER  REQ./LB.  OF  NEW  GROWTH 


MEG.  MANIP.  NO  IF  =0 


90 

90 

90. 

cr 

.  0 
0. 
30  . 
35 . 
25. 


10  0.  —  CLAY  PAN  EA  R  S  T  EN  S I  ON 

100.  • LOAMY 

100.  ---  RIPARIAN 

5.  CLAY PAN <% 

15.  LOAMY* OF 

10.  RIP<GRASS 


25. 
50. 
50  . 
70. 
3) =3  5  10 
5) =4  6  7  S  9 
,2)=1  2 

PLNTMX (ALL ) =3*1 0  0  0  0.  4*20 0 0  . 
SPB=1.0   CONSTANT  FOR  SETTING  INITIAL 


NEWRS CI ,1 
NEWRS (2,1 
NEWRS < 3,1 


jj.  — 

SO  .  - 

95.  - 
CLAY PAN 
LOAMY.  AND 
RIPARIAN 
GROUPING  OF 


COVER 

SHRUB, BARE, 
ROCK, LITTER 
CLAY PAN  <%    COVER 
LOAMY < OF  GRASS , SHRUB , 


R I P< TREE , CAN  0 P Y , G  RN  D 
WT.  %    OF  SHRUBS 
GRASSES- 


OLD  RANGE  SITES 


PLANT  BIOMASS 


SPBEXP=1.0  EXPONENT  FOR  SETTING  INITIAL  PLANT  BIOMASS 

AMEAD CALL) =621 .  829.  1225. 

HMMR=90 . 

HAYPCALL)=0 . 

THAYF-0 . 

ROOTMX  CALL) =3*1 00  0.  400.  3*1 0  0  0 

GM AX ( ALL ) - 3* 1 . 5  4*  3 

GTEMPX ('  ALL )  =40  .  80  .  10  0. 

GTEMPY(ALL)=0.  0.3  0.25  TEMPERATURE  LIMITATION        ,  n 


MAX .  ROOT  STORAGE 


w  J  •*> 


)ct  4  1; 


0 
1 

0 
1 

0 


1'  o 


40 

41 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

5S 

53 

60 

61 

62 

63 

64 

65 

££ 

67 

68 

63 

70 

71 
-n  ^-. 

{  c 

DAME 


uu 


s 

s 
s 

b 

s 
s 

s 
s 
s 

s 


s 
s 

b 


.384   da  1 3.  files-  Pag©  5 


O  r, 


n  n  f- 

u  u  5 


GSOY ( ALL) =0 .  0.2  1 .  1 .  0.5 MO  I STURE  L I M I  TAT  I  ON 

PM Y ( A  L  L ) =  0 . 0 1  0 . %  MORTAL I TY  WITH  L  OW  G  R OWT H 

PMX(ALL)=0.  0.0  5 

SPB=0 . 1  SPR I NG  PLANT  B I OMASS  (INITIAL) /PRODN . 

CWST=0.1  CON  WASTAGE 

REQGA=0.01  %GROWTH  REQUIRED  AFTER  FLOWERING 

IFDATE=13  WEEK  OF  FLOWERING 

PSTR=0.1  %  OF  SHOOT  B I OMASS  SENT  TO  ROOTS 

PRT3=0.1  %  OF  ROOT  BI OMASS  TO  SHOOTS/WEEK 

RWS=0.7  SURMIMAL  OF  ROOT  POTENTIAL  B I OMASS  OVER  WINTER 


PWS(ALL)=6*0.2 
HAYY'd ,ALL)=0 . 


0.05 P L ANT  W I  NT  E  R  S  U RM I MA L 

1 .  1  .25  1 .  0.5 


HAYY(2,ALL)=0 


1.5  1.75 
HAYY ( 3 , ALL ) =0 .  2 .  2.5 
HX(ALL)=0.  2 
USEED(ALL)=0 
HAYOP(ALL)=l 
£R(ALL)=0.34 
Siv!PX(ALL)=l  . 
SMPY(ALL)=15 
Piv'PX(ALL)=l  .  33. 
Plv'PY(ALL)=15.  3. 
CUPX(ALL)=1.  10. 
CMPY(ALL)=15.  0. 
MMANIF(ALL,ALL,ALL)=Q  . 
SHK ILL (ALL) =3*0. 8  1.0 
PAST CI . • .6) =1911 ■  1564.  3139.  3464. 
.13) =5803.  3302.  20  53.  2663 


O  a    w'  a    Ob 

30.  30.  30.  800. 

1.1.  l=POOR 

0 .36  0 . 38 
38. 

10.  SHRUB  %  PROTEIN 

33.  PERRENNIAL  GRASSES 


10.  2. 

I RR . , 2=G00D .  I RR , 3=FERT 


FORES 


YOUR  CHEATIN-  GRASS.. 


PASTC7. . 
RPB=10. 
RLF-1 . 4 


4372  . 


4638.  6448. 
1803. 


START 


ROOTS/ABO'v'E  GROUND  AT 

(GRASS+LITTER)/GRASS 

S  PC (ALL, ALL, ALL ) =0 . 

MACRO  W DARLING 

PRINT  RA( 1,2,1)  RA(1,2,2)  RA( 1,2,4)  RA(1,2,5)  RA(1,2,6)  RA(1 

* CONTINUE  WITH  DAME. SET  RETURN 

SET: 


,2,7) 


1 

s 

RA  (  1  , 1  , 1  .  .  .  8 )  = 

0  . 

13. 

16. 

3. 

73. 

1  . 

0. 

533. 

2 
3 

S 
S 

AWCNEW(1 ,1 ,ALL)= 
RA(1 ,2,1 . . .8)= 

0.13 

0  . 

0  .13 
13. 

0  .08 
16. 

3. 

73. 

1  . 

0  . 

333. 

4 
5 

S 

C 

AWCNEW(1,2,ALL)= 
RA (1,3,1. . .8)= 

0.14 
2 . 

0  .11 
42. 

0  .08 
43. 

134. 

5 . 

150. 

5. 

45. 

6 
7 

S 

s 

AWCNEW(1,3,ALL)= 
RA  (  2  , 1  , 1  .  .  .  8  )  = 

0.15 
0. 

0.15 
13. 

0.14 

16. 

3. 

73. 

1. 

0. 

782 . 

8 

9 

s 

C; 

AWCNEW(2,1,ALL)= 
RA  (  2  ,  2  , 1  .  .  .  8  )  = 

0.13 
0. 

0.13 
13. 

0.08 
16. 

3. 

73. 

1. 

0. 

782. 

10 

11 

b 

s 

AWCNEW(2,2,ALL)= 
RA  (  2  ,  3  , 1  .  .  .  8 )  = 

0.14 
0  . 

0.11 
0  . 

0.08 
0  . 

0. 

0  . 

0. 

0. 

0  . 

12 
13 

s 

b 

AWCNEW(2,3,ALL)= 
RA(3,1,1. . .8)= 

0  .0 
0. 

0  .0 
8  . 

0.0 

2  . 

11  . 

1. 

8. 

0. 

470. 

14 
15 

AWCNEW(3,1,ALL)= 
RA (3,2,1. . . 8 ) = 

0  .13 
0  . 

0  .13 

21  . 

0.08 

11  . 

3. 

U  ■ 

4. 

265b  ■ 

16 
17 

s 

AWCNEW(3,2,ALL)= 
RA  (  3  ,  3  , 1  .  .  .  8  )  = 

0  .14 

0  . 

0  .11 
0  . 

0.08 
468. 

86  . 

0  . 

132. 

0  . 

16. 

IS 
19 

c 

AWCNEW(3,3,ALL)= 
RA  (  4  , 1  , 1  .  .  .  8 )  = 

0.16 
0. 

0.16 

0.16 

11  . 

1  . 

8. 

0  . 

1113. 

Oct   4  12:00  1984   data. files  Page  6 


20  S  AWCNEW ( 4 , 1 , ALL) =  0.13  0.13  0.08 

21  S  RAC 4 , 2 , 1 . . . 8) ■     0 .    33 .  11. 

22  S  AWCNEW C  4 , 2 , ALL )  =  0.14  0.12  0.09 

23  S  RA ( 4 , 3 , 1 . . . 8 ) =     6  .   146^  2 . 

24  S  AWCNEW (4,3, ALL ) =  0.13  0.13  0.10 

25  S  RA(5,1,1...S)=     0.     9.  3. 

26  S  AWCNEW (5,1, ALL )=  0.13  0.13  0.08 

27  S  RAC 5, 2,1 . . .8)=     2.    51.  17. 

28  S  AWCNEW (5,2, ALL ) =  0.14  0.12  0.03 
25  S  RA(5,3,l..-8)=     0.   198.  119. 

30  8  AWCNEW ( 5 , 3 , ALL ) =  0.13  0.13  0.11 

31  8  RA  (  6  , 1  , 1  .  .  .  8 )  =     0  .     8  .  2  . 

32  8  AWCNEW(6,1 ,ALL)=  0.13  0.13  0.08 

33  8  RAC 6,2,1 .. .8)=     0.     0.  0. 

34  3  AWCNEW(6,2,ALL)=  0.14  0.11  0.08 

35  8  RAC 6, 3,1 . . .8)=     6.   139.  25. 

36  S  AWCNEW ( 6, 3 , ALL )=  0.13  0.13  0.10 

37  8  RAC 7 s 1 s 1 ■ • ■ 8) =     0.     9.  2. 

38  8  AWCNEW (7,1, ALL) ■  0.13  0.13  0.08 

39  8  RA ( 7 , 2 , 1 . . . 8 ) =     1 .    43 .  16. 

40  8  AWCNEWC7,2,ALL)=  0.14  0.12  0.09 

41  8  RAC 7,3,1 .. .8)=     6.   145.  2. 

42  8  AWCNEW(7,3,ALL)=  0.13  0.13  0.10 

43  SRA(S,1,1...8)=     2.    28 .  14. 

44  8  AWCNEW ( 8 , 1 , ALL) =  0.08  0.08  0.01 

45  8  RA ( 8 , 2 , 1 . . . 8 ) =     0 .    15.  48 . 

46  S  AWCNEW (8,2, ALL ) =  0.15  0.14  0.11 

47  8  RA ( 8 , 3 , 1 . . . 8 ) =     1 .    29 .  0. 
43  8  AWCNEW (8,3, ALL ) =  0.15  0.15  0.15 

49  S  RAC 9, 1,1.. .8)*     3.    29.  14. 

50  8  AWCNEW (9,1 ,ALL)=  0.08  0.08  0.02 

51  8  RA(9,2,1. . .8)=     6.    31=  30. 

52  8  AWCNEW(3,2SALL)=  0.14  0.11  0.09 

53  8  RAC 9,3,1. . .8)=     2.    36.  153. 

54  8  AWCNEWC9,3,ALL)=  0.15  0.15  0.15 


14. 

6  . 

& . 

o  ■ 

2*30 

51. 

19. 

14  . 

18. 

121 

1 2  . 

-Hi 

8. 

0  a 

1961 

17. 

10  . 

6. 

4. 

2203 

56. 

12. 

45. 

9  . 

208 

11. 

1. 

8. 

0  . 

693 

0  . 

80  . 

0  . 

0  a 

370  3 

52. 

IS. 

"jo 

£-■3  a 

17. 

302 

11. 

2 . 

8  a 

0  . 

d.  ij/ 

15. 

6. 

6 . 

5. 

3568 

51  . 

18. 

15. 

17. 

523 

18. 

11  , 

4. 

1  . 

3713 

57 , 

7  a 

5. 

0  a 

1984 

46 . 

13. 

141. 

4. 

106 

19. 

11. 

4. 

1  a 

6250 

31  . 

7. 

15. 

2a 

3495 

LIS. 

6 . 

1   f   O  « 

7. 

157 

W  I  L  D  .  D  2 

1  8  NPERPCALL)=5  3  5  —  NUMBER  OF  PASTURES  UTILIZED  BY  DEER  IN  3  PER 

I  CDS 

2  8  NPASDC1,ALL)=3  4  5  6  7  —  PASTURES  USED  BY  DEER  IN  PERIOD  1 

3  3  NPASDC2,ALL)=5  7  13  DITTO  PERIOD  2 

4  S  NPASD(3,ALL)=5  7  8  9  13 DITTO  PERIOD  3 

5  S  NPERTCALL)=2  2  3  NUMBER  OF  RANGE  SITES  USED  BY  DEER  IN  3  PERI 

ODS 

6  S  NTYPEC1,ALL)=2  3—  RANGES  SITES  USED  PERIOD  1 

7  S  NTYPEC2,ALL)=2  3 RANGE  SITES  IN  PERIOD  2 

S  NTYFE(3,ALL)=1  2  3  RANGE  SITES  USED  IN  PERIOD  3 

S  PREF (1, ALL) =3*. 1  6*5.—  PREFERENCE  PARAMETERS  FOR  PLANT  GROUPI  B 
Y  MONTH 

10  S  PREF ( 2 , ALL ) =3* .2  3*1 0  3*2 

11  3  PREF ( 3 , ALL ) =3* . 2  4* .01  2* . 2 

12  S  PREF  ( 4 ,  ALL )  =9*20  .  ,r,    r 

13  S  PREFC5,ALL)  =  .01  10  10  .01  OJOuJ 


a 
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14  S  PREF (  6  ,  ALL)  =  .  0 1  2*1 0  .01  ■'}  P\  n  C  / 

15  S  PREF (7, ALL)  =  .01  2*20  .01  '-'""  'J 
lb      S  DINGW<ALL)=12.4  20.2   IB/WEEK  ING  FAWN  AND  ADULT 

17  S  FLEAVE(ALL)=1 .   AMOUNT  OF  FORAGE  NOT  EATEN 

IS  S  DFEM=150   INITIAL  NUMBER  OF  FEMALE  ADULT  DEER 

13  S  DMALE=10  0    INITIAL  NUMBER  OF  MALE  ADUL i  DEER 

20  S  FAWN =50    INITIAL  NUMBER  OF  FAWNS 

21  S  CAL(1,ALL)=.55  .64  .73  .8b NUTRITIONAL  VALUE  FOR  EACH  PLANT  TY- 
PE BY  MONTH 

22  S  CAL( 2 , ALL) = . 8b  1.0 

29  S  FNEED=18      ;'  "  "  "  "  "  "  "  "  "  FAWN 

30  S  XMALED(ALL)=0  1 

31  S  YMALEDCALL)=0  .95   MALE  WINTER  SURVIVAL  AS  FUNCTION  OF  HEALTH 

INDEX 

32  S  XFEM(ALL)=0  1 

33  S  YFEM(ALL)=0  .3    FEMALE  WINTER  SURVIVAL  AS  FUNCTION  OF  HEALT 

H  INDEX. 

34  S  XFNY(ALL)=0  1 

35  S  YFNY(ALL)=0.  .8    FAWN  WINTER  SURVIVAL  AS  FUNCTION  OF  HEALT 

H  INDEX. 

36  S  FPERC=.4  PROP'N  OF  FAWNS   I  HAT  ARE  MALE 
S  KFRAO.25  PROPORTION  OF  GMU  COVERED  BY  SAVAL  PROJECT  AREA 
S  SHUNT=.l  PROPORTION  OF  SAGE  GROUSE  POPULATION  REMOVED  IN  FALL  B 

Y  HUNTING 

g  FRAC=.25    FRACTION  OF  DEER  POPULATION  OVERWINTERING  EAST  OF  SAVA 

L  PROJECT 

40  S  PERMIT=50G    #  DEER  HUNTING  PERMITS  IN  GMU 

41  S  DAYH=4    #  DAYS/ HUNTER/ YEAR 

42  S  HHX(ALL)=C  46  10  00  DENSITY  OF  BUCKS  ON  RANGE  FOR  SLP  FUNCTION 

43  S  HHY(ALL)=  0  .4  1    DEER  HUNTING  SUCCESS  PER  UNIT  OF  EFFORT  - 

HUNTER  DAYS 

44  S  XFAWCALL)=0  1 

S  YFAW(ALL)=0  1.2 #  FAWNS/FEMALE  AS  FUNCTION  OF  HEALTH  INDEX 

S  ISTRUT(ALL)=6  7  7  5  4   PASTURES  CONTAINING  SAGE  GROUSE  STUTTI 

NG  GROUNDS 

47  S  ISAGE=5    NUMBER  IDENTIFIED  SAGE  GROUSE  POPULATIONS  IN  MODEL 

48  S  ANX(ALL)=0  150  50  0 

49  S  ANY < ALL) =  0  .35  .8   ADULT  NESTING  SUCCESS  AS  FUNCTION  OF  FORB 

S+GRASSES  AVERAGE  DENSITY 

50  S  YNX(ALL)=  0  150  50  0 

51  S  YNY(ALL)=0  .2  .65   YEARLING  NESTING  SUCCESS  AS  FUNCTION  OF  F 

ORBS+GRASSES  AVERAGE  DENSITY; 

52  S  ANE3MN=0   MINIMUM  ADULT  NESTING  SUCCESS  IF  FORBS  DENSITY  TOO  LOW 

53  S  YNESMN=0   MINIMUM  YEARLING  NESTING  SUCCESS  IF  FORBS  DENSITY  TOO  L 
OW; 

54  S  CHX(ALL)=  0  100  50  0   CHICKS/HEN 

55  S  CHY(ALL)=0  1  5 NUMBER  OF  CHICKS  PER  HEN  AS  FUNCTION  OF  AVERAGE 

FORBS+GRASSES  DENSITY 

56  S  CHENMN=0     MIN  CHICKS/NEST 

57  S  AHEN(ALL)=2  1  8  6  38   INITIAL  CONDITIONS  ADULT  HENS  ON  STUTTIN 

G  GROUNDS  FOR  5  IDENTIFIED  POPULATIONS 

58  S  AMALE(ALL)=8  14  53  19  170   —  INITIAL  CONDI TIONSS  ADULT  MALES  ON 

STRUTTING  GROUNDS  FOR  5  IDENTIFIED  POPULATIONS 

59  S  YHEN(ALL)=1  1  4  3  15   INITIAL  CONDITIONS  YEARLING  HENS  OH    5    1 

DENTIFIED  POPULATIONS 

60  S  YMALE(ALL)=4  7  25  9  75   INITIAL  CONDITIONS   YEARLING  MALES  ON 


38 
39 


45 
46 


D 
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5  IDENTIFIED  POPULATIONS 

€1      S  XHSUR(ALL)=G  100  500   HEN  SEASON  SURVIVAL 

62      S  YHSUR(ALL)=.6  .6  .9   HEN  SURVIVAL  OVER  MARCH  15  TO  SEPT  30  AS 

FUNCTION  OF  FORBS+GRASSES 
£3      S  XHFAC(ALL)=0  5  25  50 

64  S  YHFAC(ALL)«.25  .25  1  =5  HEN  SURVIVAL  MODIFICATION  FACTOR  AS  F 

UNCTION  OF  SHRUB  COVER 

65  S  XMSURCALL)=0  5  15  30  50 

Sb      S  YMSUR(ALL)=.35  .35  .65  .65  .5  —  MALE  SURVIVAL  MARCH  15  TO  SEPT 
30  AS  FUNCTION  OF  %  SHRUB  COVER 

67  S  WHENS=.9   OVERWINTER  SURVIVAL  RATE  FOR  HENS 

68  S  WMALES=.8   OVERWINTER  SURVIVAL  RATE  FOR  MALES 

S3  S  DENERF=1E3   INITIAL  CONDITION  -  KCAL  10**3  INGESTED  BY  FAWN  DEER 

IN  PREVIOUS  YEAR 

70  S  DENERA=1E3   INITAL  CONDITION  -  KCAL  10**3  INGESTED  BY  ADULT  DEER 
IN  PREVIOUS  YEAR 

71  S  PERI N< ALL, ALL, ALL )=0 INITIAL  COND'N  -  PRCP-'N  OF  EACH  PLANT  TYP 

E  IN  DEER  DIET  BY  PASTURE  AND  RANGE  TYPE 

72  S  WE (ALL, ALL, ALL )=0  —  POUNDS/ACRE  OF  EACH  PLANT  TYPE  INGESTED  BY 
DEER  EACH  WEEK.  BY  PASTURE  AND  RANGE  SITE 

73  S  RRA (ALL, ALL, ALL) =0  DUMMY  VARIABLE  TO  REFLECT  DENSITY  OF  PLAN! 

S  AFTER  REMOVING  CATTLE  GRAZING  -  OTHER  WISE  SAME  AS  RA 

74  s  FK-4QO0   HALF  SATURATION  PARAMETER  IN  EXTRA  MORTALITY  FUNCTION  AP 
PLIED  TO  FAWN  DEER 

75  S  CMALE=.4   PROPORTION  OF  SAGE  GROUSE  THAT  ARE  MALE  IN  FALL 

76  S  SK=50  0  0   HALF  SATURATION  PARAMETER  IN  EXTRA  MORTALITY  FUNCTION  AP 
PLIED  TO  SAGE  GROUSE 

77  S  TOTHEN=0   TOTAL  NUMBER  OF  SAGE  GROUSE  HENS  MARCH  15 

78  S  TOTMAL=0   TOTAL  NUMBER  OF  MALE  SAGE  GROUSE  MARCH  15 

79  S  SAGEG  =0   TOTAL  NUMBER  OF  SAGE  GROUSE  ON  SAVAL  RANGE  MARCH  15; 

W  A  T  .  D 


1.2 

SET 

XLS(ALL)=4„5  2.8  0.16  0. 

2 

SET 

LUPTRC1. . .12)  ■ 

7*0  2*1  0 

3 

SET 

OAR  =  2.0  7 

4 

SET 

OCN  =  80 

5 

SET 

CNAC1 . . .3)  =  83 

87  91 

6 

SET 

CNBC1...3)  =  0.28  0.18  0.13 

7 

SET 

iptrci . . .io)  =  ; 

2332322112 

8 

SET 

WATRESCl. . .3)  ■ 

10  .  10  .  20 . 

9 

SET 

PMECH  ■  1 

10 

SET 

FLFUDG  =  1  . 

10  .2 

SET 

SOFDG  ■  1  . 

10.4 

SET 

RFUDG  =  1  . 

11 

SET 

ARMAH(1 . . .10) «G 

0.1  0  0  1.34  0  0  0.63  0.85  1 . 

ii7 

12 

SET 

ARGNCC1 . . .10)=0 

0.16  0  0  1.97  0  0  1.66  0.64  2 

.19 

12.05 

SET 

ARMAH(ALL)=2.62 

1.48  0.1  0 . 

12.1 

SET 

ARGNC(ALL)=4.16 

2.30  0.16  0 . 

12.2 

SET 

RUNOFF (ALL, ALL) 

=  0  . 

12.4 

SET 

SOILOF(ALL,ALL) 

=  0  . 

13 

SET 

SEDDEL  =0.2 

14 

SET 

CATFAC  =  0  .00  2 

15 

SET 

ESTB1  =  1 

let   4  12:00  1284   data. files  Page  9 


1 


lb 

SET  BSTB2  =0.2 

17 

SET  CUFS  =  27 

18 

SET  AFCON  =  2 

19 

SET  NSTEF  =  0 

20 

SET  J DAY  =65 

21 

SET  NINDAY  =10  6 

cd  ki 

SET  NWEEKS  =37 

d 

SET  SNOPAC(ALL,ALL)  =  0. 

22 » 

4 

SET  FREZDP(ALL,ALL)  =  0. 

C1.C-    m 

6 

SET  SOWAT(ALL,ALL,ALL)  =  0. 

J~>  <~i 

SET  NFFLOWC1 . . .12) =15  18  18  22  28  32  32  38  43  36  35  30 

24 

SET  NFFL0WC13. . .38) =28  23  18  10  8  5  3  2  2  2  2  2  0 

SET  PMULT  =2.3 

25. 

05 

SET  RM(ALL)  =  0 . 

25. 

1 

SET  RM(3)  =12. 

-~tC- 

2 

SET  TMULT  =  0.85 

25. 

4 

SET  TSEDG  =0. 

25. 

6 

SET  TSEDM  =  0 . 

25. 

-? 

SET  FTOT(ALL)  =  0. 

25. 

b 

S  TCOF  =  0.42 

26 

S  PRECC1 . .  .22)=0.  0.  0.2  0.12  12*0  .  0.1  0.12  0.01  3*0. 

27 

S  PRECC2S. . .64) =0.15  0.3  0.5  26*0.  0.4  0.7  11*0. 

28 

S  PREC ( 65  ...  82 ) =0 . 0  5  0.04  4*0 .  0.01  0.0  3  0 . 1 6  6*0 .  0.1  0 .  0 . 

23 

S  PREC (83. . .104)=0 .22  0.03  0.04  16*0.  0.04  0.03  0.22 

30 

S  PRECC105. . .120)=0 .  0.  0.1  7*0.  0.16  0.03  0.01  3*0. 

31 

S  FRECC121 . . .142)=0 .04  0.05  0.  11*0.  0.4  0.1  0.22  0.  0.  0.14  0.16 

Oil 

S  PREC ( 1 43  ... 1 81 ) =7*0  .  0.08  8*0  .  2*0  .1  6*0 .  0.02  0.04  0.14  1 2*0 . 

w'C> 

3  PRECC182. . .210)=2*0 .1  3*0.  0.02  13*0.  0.14  0.04  0.1  0.  0.1 

34 

S  PREC ( 211 .. .254 )=31*0 .  0.4  12*0. 

35 

S  PREC ( 255 . . . 236 ) =0 . 1  8*0 .  0.1  0.18  1 0*0 .  0.02  4*0 .  2*0 . 1  1 4*0 . 

36 

S  PREC (237. . .313)=0 .  .14  0.04  0.02  6*0.  2*0.1  11*0. 

37 

S  PREC( 320. . .333) =0.04  0.03  0.22  0.  0.1  6*0.  0.16  0.03  0.01 

38 

S  PREC ( 334 . . . 365 ) =4*0 .  0.04  0.05  22*0 .  0.7  0.4  2*0 . 

33 

3  TMAXLC1. . .10)=36  40  40  42  41  34  31  38  33  27 

40 

S  TMAXLC11 . . .20) =17  30  42  37  40  33  32  20  24  33 

41 

8  TMAXU21.  .  .30)=35  30  40  45  35  31  23  24  24  21 

42 

S  TMAXU'31  .  .  .40)  =27  40  42  44  44  42  42  31  38  38 

43 

S  TMAXU41  .  .  .50)=41  44  42  41  42  33  38  40  43  37 

44 

S  TMAXL(51. . .60)=35  37  35  33  44  43  52  55  44  43 

45 

S  TMAXL(61 . . .70)=42  38  38  33  33  36  36  33  42  46 

46 

S  TMAXU71.  .  .80)  =37  30  41  48  34  23  46  35  37  47 

47 

S  TMAXL(81 . . .30)=32  33  38  33  37  38  28  33  49  26 

48 

8  TMAXLOl.  .  .100)=31  36  38  43  50  37  37  38  41  43 

43 

8  TMAXLdOl  .  .  .110)=40  45  43  54  62  56  65  70  63  72 

50 

S  TMAXU'lll . . .120)=74  67  53  54  45  53  64  67  67  60 

51 

3  TMAXL(121 . . .130) =53  58  65  64  64  63  60  53  55  45 

52 

S  TMAXL(131 . . .140)=43  48  51  56  55  55  43  53  65  63 

53 

S  TMAXL(141. . .150)=75  73  71  40  35  33  47  57  57  54 

54 

S  TMAXLC151 . . .160) =55  51  51  40  50  50  52  55  66    74 

55 

S  TMAXU161.  .  .170)  =75  73  68  60  62  56  57  71  81  73 

S  TMAXU171  .  .  .180)  =78  80  78  76  64  76  76  70  70  80             Q00» 

56 

57 

S  TMAXU181  .  .  .130)  =81  67  67  72  72  77  7/    80  ,'  3  73 

53 

S  TMAXLC131 . . .233) =83  73  84  73  76  77  82  8b  82  86 

53 

S  TMAXU201.  .  .210)  =32  85  32  34  33  30  30  93  32  S3 

60 

S  TMAXL(211 . . .220) =85  86  84  85  83  76  82  85  84  81 

61 

6  THAXL(221 . . .230 )=86  82  86  83  87  88  85  67  77  81 

4    12:00    1384       dat. 


files    Pa': 


10 


62 

S  TMAXL(231 . . 

.240) =64 

64 

76 

82 

80 

80 

■•■'  ':■• 

/  3  b^  i-  3 

63 

S  THAXLC241 • ■ 

.250) =75 

71 

63 

68 

76 

73 

i  © 

o  •_'  o  ■-'  b  1 

64 

S  TMAXLC251 . . 

.  260 ) =62 

f   £l 

69 

52 

64 

68 

67 

70  74  75 

65 

S  TMAXLC261  .  . 

.270) =80 

7d 

63 

62 

56 

67 

69 

71  76  77 

66 

S  TMAXU271  .  . 

. 2S0 ) =78 

/  £. 

"7  "7 

bk' 

73 

"7  "7 

81 

73  76  73 

67 

S  TMAXLC  281 .  . 

.230) =78 

7o 

7  7 

63 

55 

40  36  38 

68 

S  TMAXU231  .  . 

.300) =56 

59 

61 

60 

62 

53 

63  51  40 

63 

S  TMAXU301  ■  . 

. 31 0 ) =38 

50 

57 

62 

64 

ETC1 

DO 

59 

65  67  69 

70 

S  TMAXL(311 . . 

.320) =62 

65 

49 

47 

55 

36 

37 

38  37  35 

71 

3  TMAXU'32'1  .  . 

. 330 ) =38 

47 

40 

42 

42 

41 

33 

30  26  41 

72 

S  TMAXLC331 . . 

.340) =41 

42 

48 

43 

41 

45 

46  35  26 

73 

S  TMAXL( 341 . . 

.350 )=25 

17 

■i!  f 

oJ 

34 

44 

46 

46  46  48 

74 

S.TMAXLC351 . . 

. 360 ) =52 

51 

41 

4o 

42 

bb 

41 

40  43  43 

75 

s'TMAXL( 361 . . 

.365) =43 

59 

43 

48 

51 

76 

8  TM I NL « 1 . • . 1 0 ) =1 5  8  16  £ 

24 

21 

16 

26 

27 

1 2 

77 

S  TMINLC11 . . . 

20) =11  11 

3C 

31 

•ne 
c~  — 

IE 

11 

—3  -7 

73 

3  TMINL( 21 . . . 

30 ) =4  6  2 

7 

8  2 

-1 

-1 

3  - 

16 

73 

S  TMINLC31 . . . 

40) =-8  IS 

:  10 

2fc 

i  1£ 

;  is 

2= 

'  17 

'  14  5 

SO 

S  TMINLC41  .  .  . 

50) =11  13 

-tC 

i   dw 

2A 

2£ 

S.  .' 

24 

2E 

30  30 

81 

S  TMINLC51  .  .  . 

60) =20  21 

IE 

!  15  12 

'  3C 

27 

'd  8  2  'd 

82 

8  TMINLC61  .  .  . 

70 ) =20  28  25 

27 

£lw 

.  2C 

2C 

2C 

22  20 

S3 

8  TMINU71  .  .  . 

80) =13  13  23  27 

'12   £ 

1  12 

!  2i 

I    18  17 

85 

8  TMINLC81. . . 

90  )  =21  23  22 

:  H 

i    12 

1   E 

.  IE 

(  13  16  12 

86 

8  TMINLC31  .  .  . 

100)=12  20  1 

4  1 

7  25  24  20  1 

2  16  17 

87 

8  TfilNLClOl  .  . 

.110) =22 

25 

'"JO 

31 

26 

19 

27  23  28 

88 

S  TMINLdll  .  . 

.120) =31 

C1  -w1 

•-ir 

23 

35  34  30 

89 

S  TMINLC121. . 

.130) =23 

38 

32 

3b 

35 

34 

38 

37  33  30 

90 

S  7MINLC131 . . 

.140)=  30  23  28  32  3fc 

>  32  25  32  33 

31 

8  TMINLC141  .  . 

.150) =36 

43 

33 

30 

27 

iC  c-- 

30 

27  23  34 

32 

3  TMINL(151. 

.160) =36 

i'b 

>•.'  £- 

30 

30 

■d'd1 

29  b  k'  b  7 

33 

S  TMINL(161. , 

.170) =40 

42 

42 

34 

32 

31 

34 

37  42  45 

94 

8  TMINL(171 . 

.180) =40 

46 

40 

c!b 

33 

b2 

39 

33  27  37 

35 

8  TMINLC1S1. , 

.190)=44 

50 

51 

46 

44 

38 

37 

36  45  43 

36 

S  TMINLC191. 

.200) =41 

42 

38 

42 

43 

3  / 

41 

2  48  41 

37 

S  TMINLC2Q1 . 

.210)=41 

45 

48 

53 

51 

49 

46 

47  48  43 

98 

8  TMINLC211 . 

. 220 ) =53 

48 

51 

43 

41 

51 

bb 

40  47  44 

33 

8  TMINLC221 . 

.230) =45 

42 

38 

40 

45 

43 

45 

43  33  40 

100 

8  TMINLC231. 

,  .240) =35 

35 

31 

37 

43 

3b 

33 

33  40  45 

101 

8  TMINU241. 

.250) =34 

36 

41 

35 

31 

36 

38 

38  41  50 

102 

8  TMINL(251 . 

■ . 260 ) =41 

33 

41 

44 

40 

3'tL 

42 

36  36  36 

103 

S  TMINU261. 

,  . 270 ) =35 

42 

34 

23 

27 

21 

26 

27  23  33 

104 

S  TMINL(271 . 

. .280) =34 

42 

31 

35 

41 

31 

3b 

3b  34  35 

105 

8  TMINU281. 

. . 230 ) =32 

31 

40 

33 

33 

29 

13 

22  26  26 

106 

8  TMINL(231 . 

.  .300) =24 

23 

25 

22 

2  b 

20 

7 

14  25  26 

107 

8  TMINL(301 . 

.  .310) =23 

16 

11 

16 

24 

23 

30 

31  26  25 

108 

8  TMINU311  . 

. .320)=31 

33 

25 

23 

34 

31 

25 

21   2  11 

103 

S  TMINLC321. 

. .330)=  3 

3 

8 

12 

7 

6 

26 

13   1  -3 

111 

S  TMINLC331 . 

.  .  3  4  U  )  —    5 

21 

30 

20 

3 

31 

30 

11   6  8 

112 

8  TMINL(341 . 

. .350)=  8 

1  • 

-8  ■ 

-2  ■ 

-4 

2 

4  0  -1  13 

113 

8  TMINU351  . 

. . 360 ) =22 

20 

20 

20 

18 

19 

dO 

2b  20  23 

114 

S  TMINLC361 . 

. .365) =23 

26 

25 

19 

20 

n  n  ^ n 
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1 

C 


M    M    0    N      B    L    0 


A  V  A  L 


This  is  the  main  common  block,  which  is  called  by  the  initialization 
routine  UINIT.  Common  block  files  associated  with  subroutines  are 
attached  to  SAVAL.C  by  means  of  ^CONTINUE  statements;  these 
files  are    listed  according  to  their" order  of  occurrence  in  SAUAL.C. 

1.2  ^CONTINUE  NITH  DAVE.C  RETURN 

2  ^CONTINUE  NITH  BEEF . C  RETURN 

4  * CONTINUE  NITH  NAT . C  RETURN 
4.2  ^CONTINUE  NITH  NILD.C2  RETURN 

5  C  OMM  ON  I N EEK , I M ONT  H 

S  C  OMM  ON  I M  E  G , I N I L  D , I N AT , I B  E E  F 

7      ^CONTINUE  NITH  CHON.C  RETURN 


o 


* CONTINUE  NITH  GRON.C  RETURN 


DAVE 


wl 


C  OMM  ON  G  C  OM  P  (  9  ,  3  ,  7  ) 

C  OMM  ON  RAN  EN ( 9 , 3  , 1 4 ) 

C  OMM  ON  T  OT  F  B ( 9 , 3 )  , P OT  G , R  Q E V I N , A M E V I N , RQEVLB ( 7 ) 

COMMON  I , IP, J,K 

COMMON  ANCNEN ( 9 , 3 , 3 )  , I  MM 
£  COMMON  I  LUMP 

7  COMMON  NRT , NPAS , I OC ( 1 5 ) , AMEAD ( 3 ) , HMMR , PLNTMX ( 7 ; 

3  '  C  OMM  ON  N  EN  R  S ( 1 0  , 1 0  )  , S F  B EX  P 

9  COMMON  GMRED (9,3,7), GPCNT (13,7), ROOTMX ( 7 ) 

10  COMMON  RPB,RLP 

11  COMMON  ANC(10,3)  , RAO, 3, 14)  ,  SPB  ,  ROOT(  9  ,  3  ,  7) 

1 2  COMMON  VMAN I P ( 9 , 3 , 4 ) , PAST (13), SHK I LL ( 4 ) , PREM 

1 3  COMMON  ARET , Y SEED ( 7 )  , THAYP , HY ( 5 ) , HAYY ( 3 , 5 )  , HX ( ! 

1 4  -  COMMON  HAYOP ( 3 )  , NPA ( 3 ) , HAYP ( 3 ) , PSTR , PRTS , RNS 

15  COMMON  TPR0D(13,7) ,TRED , REQGB , REQGA , I FDATE 

IS  C  OMM  ON  GT  EM  PX ( 3 ) , GT EM P Y ( 3 ) , NM AX , G  P  L ANT (9,3,7) 

1 7  COMMON  GMAX ( 7 ) , GSOX ( 5 )  , GSOXM ( 1 0 , 5 ) , GSOY ( 5 ) 

18  C  OMM  ON  PM  R  ,  PMX  (  2  ),  PMY  (  2  ),  CN  ST  ,  D  R  0  OT  ,  D  P  B 

1 9  COMMON  PNGR , PNS ( 7 )  , PTOT , PCNT ( 9 , 3 , 7 )  , UP ( 7 ) 

20  COMMON  SUPX  (  2 )  ,  SVPY  (  2 )  ,  PUPX  (  2 )  ,  FVPY  (  2 )  ,  CU PX  ( 2 ) 

21  COMMON  CVPY ( 2 ) , PC ( 9 , 3 , 5 ) , FNTMAX (10,2), PCMAX ( 1 0 

22  COMMON  PCMAX2(10,5) ,TOTPCT , ER( 10 ) ,GTPROD(13) 


B  E  E  F  .  C 


pi  ' 


\R9 


i;  V  w  "  ' 


1  C  OMM  ON  C  ON  (  3  )  ,  C  ONN  (  3 )  ,  C  A  L  F\-\ ,  HAY  B  ,  HAY  S  ,  H  A  YT 

2  COMMON  I  YPN (9 )  , I PTO ( 9 )  , I PCNT ( 9 ), I GR ( 9 ), I PN ( 9 ), I PREF ( 7 ) 

3  COMMON  DPN , Y I  I  , SL I  , TTREQ , TREQ ( 3 )  , I GRT , I  PAST ( 9 )  , I  USE 

4  C  OMM  ON  I  NM  , C  E ( 9 , 3 , 7 )  , B EAT ( 3 )  , YT DN ( 7 )  , S  LT  DN ( 7 ) 

5  C  OMM  ON  YN  G ( 5 )  , N GA I N ( 3 )  ,  F  L  E  FT  , F  0  L  D , F  STA  RT , FT  H  R 


12 : . 


1 9 8 4   co  mm  o  n , bl o  c  k  Page 


6 
7 
8 

Q 

10 

11 

12 

13 
14 
15 
16 
17 
18 


,DPWN, IAHS 


C GMM  ON  X  C L  F ( b )  , Y  C L  F ( 6 )  , C A  L  F , F  C S , P A  S , H  P  C 

C  OMM  ON  I A  H B , HAY  R ( 3 )  , T  HAY R , XWN G ( 3 ) , Y WN G ( 3 ) 

COMMON  TON ( 3 ) , XWG ( 5 ) , I PYN ( 3 )  , NGAR (8,6) 

C  OMM  ON  CA  L  F  S  ,  C  ON  S  ,  Y  E  A  R  S  ,  PP  C  ON  ,  F  P  C  L  F  ,  P  P  Y  R 

C  OMM  ON  C  0  ST  S  ,  D  G  H  ,  F  C  ,  U  C  C  ,  DAN  ,  T  G N  A  ,  T  6WA  0 

COMMON  REU  ,  DSH  ,  RNP  ,  Gk-'O  ,  XM1  ,  CUA  ,  XM2  ,  RNA  ,  XM3 

COMMON  BEEF, P CALF 

C  OMM  ON  Y  R  B  ,  C  L  F  B  ,  C  ON  B  ,  D  B  H 

COMMON  NPAST 

COMMON  YEARN , NPMAX , CU (3 ) , CONMAX , HPF ( 3 ) , CM I N ( 3 ) , THAYG 

COMMON  I  EN R, I  ST R 

C  OMM  ON  C ONNMN , I  ST  0 P , I E  C ON , P  B  S  L  D 

COMMON  BUY HAY 


00060 


W  A  T  .  C 


4 

5 

7 

8 

3 
10 
11 
12 
13 
14 
14, 
14, 
14, 
16 
16 
16, 
17 
18 


COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 
COMMON 


PREC(365) ,  XLS (10) 

TM AX  L ( 3 6  5 )  ,  TM I N L ( 3 6 5 )  ,  CN A ( 3 )  ,  UN B ( 3 ) 
SNOPAC (3,3), DPREC (13),  DPRECT (13),  SONAT (9,3,3) 
NATRES(3),  S0IL0F(3S3),  RUNOFF (3, 3) 
ARMAH (10),  ARGNC ( 10),  FTOT ( 3 ) ,  FM I N ( 3 ) ,  NFFLOW ( 3£ 
NTN  O  O  D  ,  B  STA  B ,  B  ST  B 1 ,  B  ST  B  2 ,  TM AX (13) 
OCN,  OAR,  PMULT,  TMULT ,  FLFUDG ,  ORUN ' 
U E  GN  C ,  U EMA  H ,  C AT  F A  C ,  BETOT ,  W A l Q A L ,  N  Q UA  L 
CUFS,  AFCON,  I  PAS,  I RTYPE ,  RADLAN,  ETPOT 
PACAN,  ETPLNT,ETSOIL,  I  LEAK,  IFIL,  AFAR,  BPAR 
S,  SMELT,  DMT EM P,  NAT IN,  RFAC ,  CFAC 


OS,  RM(13) 

PCDUM2,  FREZDP(9,3) 
NINDAY,  NNEEKS 


FLGNC,  FLMAH, 

PCDUM,  PCDUM1 
SEDDEL ,  PMECH , 
SO (3, 3 ,3) 

NSTEP,  J DAY,  NDAYS, 

LUPTR(13),  IFTR(lU) 

SOFDG,  FNT0T(3) 
PRECIP(13),  RFUDG,  TSEDG, 

FEVAFO,3),  ETPREM,  TCOF , 


IWFLG,   I  DAY 


TSEDM 
ETSUM 


W  I  L  D  .  C  2 


1 
2 
3 
4 
5 

6 

7 
8 

9 
10 
11 


■3) 


COMMON  HEALA , HEALF , DENERA , DENERF , ANEED , DEER , FNEED , DMALE 

COMMON  XMALED ( 2 ) , YMALED ( 2 ) , DFEM , XFEM ( 2 ) , YFEM ( 2 ) , FAWN , XFNY ( 2 ) 

COMMON  I TR I G ( 3 ) , FFERC , YFNY ( 2 ) 

COMMON  I C( 3) , DSEAS , FRAC , HDAYS , HFRAC , PERM I T , DAYH , KI LL 

COMMON  HHX ( 3 )  , HHY ( 3 )  , TK I LL , XFAN ( 2 )  , YFAN ( 2 )  , I  FLAG , TOT I N , NPERF ( 

COMMON  NPASD(3,13) 

COMMON  NPERT ( 3 ) , NTYPE ( 3 , 1 0 ) , RRA ( 1 3 , 3 , 7 ) , TOTVEG , PREF (7,12) 

COMMON  PER I N ( 1 3 , 3 , 7 ) , TDEMDA , D I NGN ( 2 ) , TDEMDF 

C  OMM  ON  F  L  E AU  E ( 7 )  , CA  L ( 7 , 1 2 ) 

COMMON  I  STRUT ( 5 )  , FA1 Jl 5 ( 5 )  , F Jl SI ( 5 )  , G Jl SI ( 5 )  , CM1M31 ( 5 ) 

C OMM  ON  FM 1583(5)  , CM 1 5  S 3 ( 5 )  , GM 1533(5), ANX ( 3 )  , ANY ( 3 )  , AN  E  ST ( 5 )  , A 
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1 2  COMMON  YNEST ( 5 ) , YNX  C  3 )  , YNY ( 3 )  ,  CHEN ( 5 ) , CHX ( 3 )  , CHY ( 3 ) , CHENMN 

1 3  COMMON  CH I CK ( 5) , AHEN ( 5 ) , YHEN ( 5 ) , SURUH ( 5 ) , XHSUR ( 3 ) , YHSUR ( 3 ) , FA 

C 

1 4  COMMON  XHFAC ( 4  )  , YHFAC ( 4 )  , SURU-1 ( 5 )  , XMSUR ( 5 )  , YMSUR ( 5)  , AHENF ( 5 ) 

15  C  OMM  ON  AM A  L  E ( 5 ) , YM A  L  E  <  5 )  ,  W  H  EN  S  , WM A L  E  S , CM A  L  E , YN  E  3MN ,  S  A  G  E  G 

16  COMMON  TOTHEN,TOTMAL 

17  C  OMM  ON  GA 1 J 1 5 ( 5 )  , AM A  L  E  F ( 5 )  , S K , XM  FAN , F  K 
1 3  COMMON  WE ( 1 3 , 3 , 7 )  , SHUNT , I  SAGE , TDEER 

C  H  0  W  .  C 


G  R  0  W 


ON 


COMMON  H(3) ,ALPH(3) , PREFF( 7 , 3 , 2) ,TTAT( 3, 7,3,2) ,TDIE(7,3,2) 
COMMON  TPREF , FPREF ( 7 ) , WPREF ( 2 )  , RPREF  C  3 ) 
COMMON  I MUNCH 


COMMON  EM „ EMP , EMN , EMNP , GEN , GEN PI , GENP2 , FP1 , FP2 

COMMON  REMN , PM1 , PM2 , R I N , SURP , DNEG , Gl( 2 ) , G2 ( 2 ) , G3 ( 2 ) , G4 ( 2 ) , I GR 
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II.   Integrated  Sage  Grouse/Walker  Vegetation  Model 
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INTRODUCTION 


Following  the  first  Saval  workshop  in  November,  1981,  it 
was  recognized  that  there  was  a  need  for  a  less  complex  model  of 
key  system  interactions.   Therefore,  a  simple  model  was  constructed 
to  simulated  interactions  between  grazing,  vegetation,  and  sage 
grouse  on  an  APPLE  microcomputer.   The  rules  for  grazing  and 
vegetation  dynamics  are  drawn  from  Walker  et  al.  (1981),  which  is 
included  here  to  facilitate  understanding  of  the  program  listing. 
The  program  is  included  on  the  diskette  at  the  back  of  this  report. 

The  model  runs  on  an  annual  time  step,  tracking  changes  in 
the  numbers  of  sage  grouse  (adult  hens  and  males) ,  shrub  biomass 
and  growth  rate,  grass  biomass  and  growth  rate,  grass  consumption 
by  herbivores,  and  the  infiltration  rate  of  water  into  the  soil. 


n  o  n  <  7 
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SUMMARY 

(1)  Semi-arid  savannas,  wherever  they  occur,  have  generally  been  overgrazed  and 
encroached  on  by  bush.  A  model  is  developed  which  accounts  for  the  growth  of  woody 
vegetation  and  of  grasses,  and  analyses  the  competition  between  them  for  available  soil 
water. 

(2)  The  model  is  based  on  Walter's  two  layer  hypothesis.  Woody  vegetation  and 
grasses  compete  for  water  in  the  surface  layers  of  the  soil,  but  woody  vegetation  has 
exclusive  access  to  a  source  of  water  relatively  deep  underground.  Where  there  is  only  a 
small  biomass  of  grass  the  soil  surface  tends  to  become  impermeable  and,  in  these 
conditions,  the  model  shows  that  two  different  steady  states  may  develop:  with  a  lot  of 
woody  vegetation  alone,  or  with  a  relatively  large  biomass  of  grass  and  rather  little  woody 
vegetation. 

(3)  The  results  are  discussed  in  terms  of  the  concept  of  resilience.  The  continued 
existence  of  both  stable  states  under  ranching  conditions  seems  to  depend  on  periodic 
heavy,  or  over-,  grazing  which  allows  for  the  maintenance  of  unpalatable  or  unstable 
grass  species,  which  thus  set  a  minimum  to  grass  biomass — a  minimum  which  cannot  be 
reduced  by  herbivores. 

(4)  Comparison  of  the  dynamics  of  various  savanna  and  other  natural  systems  leads 
to  the  conclusion  that  the  resilience  of  the  systems  decreases  as  their  stability  (usually 
induced)  increases. 

INTRODUCTION 

In  this  paper  the  term  'semi-arid  savanna'  refers  to  those  regions  of  the  world,  which,  in 
their  natural  state,  have  a  predominant  continuous  grass  cover,  with  scattered  to  numerous 
trees  and  shrubs.  Plant  growth  is  water-limited.  The  term  does  not  include  those  areas  that 
would  develop  into  woodland  in  the  absence  of  fire  or  other  disturbance. 

A  series  of  radical  changes  occurred  in  these  semi-arid  savannas,  beginning  in  the  late 
1 800s  when  settlement  of  the  drier  parts  of  the  tropics  and  sub-tropics  by  Europeans  was 
at  its  peak.  These  changes  have  continued  and  show  a  remarkable  similarity  wherever 
semi-arid  savanna  is  found.  The  general  pattern  of  events  was  initiated  by  the  introduction 
of  domestic  cattle  or  sheep,  or  both,  and  resulted  in  a  drastic  reduction  in  numbers,  or  even 
elimination,  of  indigenous  large  mammalian  herbivores  and  their  predators. 

The  populations  of  cattle  and  sheep  increased  rapidly,  exceeded  the  carrying  capacity  of 
the  vegetation,  showed  increased  mortality  (usually  in  conjunction  with  a  period  of 
drought),   and   crashed.   They   then   recovered   somewhat,   and   fluctuated   around   a 
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considerably  lower  level  than  was  attained  at  their  peak.  This  process  was  accompanied  by 
marked  changes  in  the  vegetation,  usually  beginning  with  a  change  in  grass  species 
composition,  followed  by  a  decline  in  the  perennial  grass  cover,  and  an  increase  in  the 
biomass  of  annual  grasses  and  woody  vegetation.  In  the  most  extreme  cases,  woody 
vegetation  in  the  form  of  scrub  or  thicket  now  dominates  the  final  state,  virtually  excluding 
grasses. 

Even  if  all  the  large  herbivores  are  then  removed  the  vegetation  generally  remains  in  its 
new  state,  at  least  for  20  yr,  and  does  not  revert  to  its  original  state.  This  general  story  has 
been  described  for  many  areas  in  the  Sahel  zone  of  Africa  (Wickens  &  Whyte  1979)., 
southern  and  eastern  Africa  (Barnes  1979),  the  south-western  U.S.A.  (Buffington  &  Herbel 
1965;  Clark-Martin  1975),  the  drier  parts  of  India  (Singh  &  Joshi  1979)  and  Australia 
(Wilson  &  Graetz  1979).  In  most  cases  the  collapse  is  triggered  by  drought,  and  there  is 
generally  increased  soil  erosion  and  reduced  infiltration  of  water  into  the  soil.  Fire  is  not  an 
important  factor  and  seems  not  to  contribute  to  the  change  (in  contrast  with  the  wetter 
savannas). 

Our  objectives  are  to  examine  the  dynamics  of  these  systems  and  to  understand  how 
stable  they  are  likely  to  be,  and  how  they  will  behave  under  different  forms  of  management. 
We  begin  with  a  general  case,  based  on  southern  African  experience,  and  use  this  to  show 
how  different  situations  can  arise  from  inherent  differences  in  the  environment  or  structure 
of  the  system,  or  from  differences  in  management. 


DESCRIPTION  OF  THE  SYSTEM  AND  NECESSARY  ASSUMPTIONS 

The  many  kinds  of  woody  vegetation  can  be  broadly  grouped  into  the  microphyllous, 
xerophytic  Acacia  forms  and  deciduous,  longer-lived  broadleaved  forms  (perhaps  best 
typified  by  the  family  Combretaceae).  The  major  components  of  the  herbaceous  layer  are 
forbs,  annual  grasses  and  perennial  grasses.  These  details  can  be  simplified  for  our 
purposes,  however,  and  the  vegetation  regarded  as  a  mixture  of  woody  vegetation  and  of 
perennial  grasses,  maintained  in  a  stable  equilibrium  by  the  fact  that  they  tap  different 
sources  of  water  in  the  soil  (Walter  1971).  Grasses  are  more  efficient  than  trees  in 
extracting  water  from  the  upper  layers  of  the  soil  (topsoil),  but  below  the  grass  root  zone 
(in  what  is  here  called  the  subsoil)  the  woody  vegetation  has  nearly  exclusive  use  of 
whatever  water  gets  through.  Walker  &  Noy-Meir  (1981)  have  explored  this  hypothesis 
further  and  have  shown  that  in  some  conditions  a  single,  stable  equilibrium  of  woody 
vegetation  and  grasses  exists.  These  conditions  are  that  the  proportion  of  a  tree's  water 
requirements  obtained  from  the  subsoil  is  sufficiently  high,  and  that  rainfall  is  enough  for 
some  water  to  infiltrate  to  the  subsoil,  but  not  so  much  that  woodland  can  develop. 

Because  semi-arid  savanna  is  a  water-limited  system,  any  reduction  in  the  rate  of  water 
infiltration  into  the  soil  is  critical.  On  medium  to  heavy-textured  soils,  infiltration  decreases 
as  the  surface  pores  become  sealed  ('capping')  and  soil  erodes,  and  both  of  these  processes 
increase  in  rate  as  grass-cover  and  litter-cover  decline.  Kelly  &  Walker  (1976)  have 
demonstrated  that  the  rate  and  amount  of  infiltration  into  a  loamy  savanna  soil  is  about  ten 
times  greater  under  a  grass-litter  cover  than  it  is  through  a  bare  soil  surface.  During  heavy 
rainfall,  the  difference  is  crucial  and  leads  to  run-off  from,  and  erosion  of,  bare  soils. 
Similar  reports  have  been  made  for  the  semi-arid  ranges  in  the  south  western  U.S.A.  (e.g. 
Glendenning  1952;  Buffington  &  Herbel  1965).  The  biomass  of  grass  is  therefore  a  critical 
factor  in  determining  the  rate  and  amount  of  infiltration. 
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On  the  other  hand,  woody  plants  can  influence  the  water  content  of  the  subsoil.  The 
canopy  intercepts  rain  and  funnels  it  down  the  stem  and  into  the  soil  at  the  tree  base. 
Pressland  (1973),  working  in  an  Australian  savanna,  has  recorded  a  six-fold  increase  in  the 
amount  of  water  percolating  into  the  soil  below  trees  as  compared  to  open  areas,  and 
Reynolds  (1966)  has  noted  a  similar  phenomenon  in  a  British  woodland.  Further  means  by 
which  water  penetrates  preferentially  to  the  subsoil  are  the  quick-flow  channels  or  tongues 
of  preferred  flow  associated  with  wetting  front  instability,  (P.  de  Vries,  personal 
communication).  As  a  result  of  these  channels,  water  moves  rapidly  into  the  subsoil  before 
the  topsoil  is  at  field  capacity.  The  net  result  is  that,  although  less  water  will  infiltrate  the 
topsoil  if  grass  cover  declines,  there  is,  rather  surprisingly,  a  relatively  smaller  reduction  in 
the  amount  of  water  reaching  the  subsoil.  If  the  reduction  in  infiltration  is  not  too  drastic, 
then,  as  the  amount  of  grass  declines,  less  water  will  be  taken  up  from  the  topsoil  by  plants, 
and  so  more  will  penetrate  to  the  subsoil. 

Rainfall  is  here  assumed  to  be  generally  constant  but  interrupted  by  periodic  droughts, 
which  result  in  the  death  of  many  plants,  especially  those  grasses  and  woody  plants  in  the 
process  of  establishment. 

Cattle  and  sheep  are  the  major  herbivores  and  eat  negligible  amounts  of  woody 
vegetation.  In  Africa,  at  least,  cattle  do  not  inhibit  woody  plant  regeneration  (Taylor  & 
Walker  1978)  and,  if  anything,  increase  the  rate  of  dissemination  of  seeds  and  pods.  Under 
high  stocking  levels  the  common  pattern  of  range  deterioration  is:  reduced  biomass  of 
grass,  leading  to  reduced  infiltration,  which  further  reduces  the  rate  of  growth  of 
vegetation,  particularly  of  grasses.  The  decrease  in  the  biomass  of  grass  is  accompanied 
by  an  increase  in  that  of  woody  vegetation. 

This  shrub-dominated  state  is  remarkably  stable  and,  if  the  cattle  are  removed,  does  not 
revert  at  once  to  open  grassland.  Over  a  long  period  of  time  (how  long  is  not  known)  the 
system  would  probably  slowly  change,  with  a  gradual  improvement  in  the  soil  surface 
conditions,  which  would  re-establish  the  predominance  of  topsoil  water  and  allow  an 
increase  in  the  biomass  of  grasses. 

PROCEDURE 

Our  analysis  consisted  of  two  phases.  First,  we  identified  the  important  variables,  and  used 
graphs  to  understand  the  factors  which  govern  their  dynamics.  Secondly,  we  expressed  the 
behaviour  of  the  system  algebraically,  and  manipulated  it  mathematically  to  allow  a  better 
understanding  of  the  quantitative  response  of  the  system  to  changes  in  the  values  of  the 
variables  and  parameters. 

The  mathematical  procedure  was  similar  to  that  described  by  Ludwig,  Jones  &  Holling 
(1978).  The  state  variables  were  divided  into  'fast'  and  'slow'  categories,  according  to  the 
time-scale  over  which  they  change.  At  first,  the  values  of  the  slow  variables  were  fixed,  and 
the  fast  ones  allowed  to  reach  steady  state  values.  Thus  the  potential  equilibria  of  the  fast 
variables  could  be  expressed  in  terms  of  the  slow  ones.  This  procedure  greatly  simplified  the 
subsequent  analysis  of  the  dynamics  of  the  combined  system. 

SELECTION  OF  STATE  VARIABLES 

In  southern  Africa  there  are  five  essential  variables:  topsoil  and  subsoil  water  content  (cf. 
Walker  &  Noy-Meir  1981),  the  biomass  of  grass,  of  woody  vegetation,  and  of  herbivores. 
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It  would  be  possible  to  subdivide  the  vegetation  variables  since  it  is  well  established  that  (at 
least)  there  are  differences  in  growth  and  survival  rates  of  perennial  and  annual  grasses 
(Kelly  &  Walker  1976),  in  the  payability  to,  and  consumption  by,  herbivores  of  various 
grasses  (Gammon  &  Roberts  1978;  Kennan  1969)  and  in  the  susceptibility  to  drought  and 
competition  of  establishing  and  established  plants  (Glendenning  1952).  There  is  no  clear 
evidence  as  to  how  these  differences  affect  the  final  outcome  of  the  woody  vegetation-grass 
combination,  however,  and  we  therefore  proceed  with  just  the  two  broad  categories  and 
consider  the  differences  as  specific  effects  at  a  later  stage. 

A  simple,  but  still  effective  representation,  does  not  need  all  the  five  variables  explicitly. 
For  example,  the  two  most  rapidly  changing  variables,  topsoil  and  subsoil  water,  are 
totally  dependent  on  other  variables:  they  depend  on  infiltration  rate,  on  the  biomass  of 
grass  and  of  woody  vegetation,  and,  to  a  much  lesser  extent,  on  leaching  and  evaporation 
losses.  The  effects  of  changes  in  the  amounts  of  topsoil  and  subsoil  water  can,  therefore,  be 
expressed  by  changes  in  the  biomass  of  grasses  and  of  woody  vegetation,  and  of  infiltration 

rate. 

Infiltration  rate  itself  depends  largely  on  the  amount  of  grass  cover  at  any  one  site,  but  it 
is  convenient  to  retain  infiltration  rate  in  the  analysis.  There  is  some  (empirical)  evidence 
that  under  conditions  of  low  grass  cover,  a  moderate  amount  of  animal  hoof  action  breaks 
up  soil  surface  crusts  and  can  increase  infiltration.  This  effect  is  difficult  to  quantify  and  is 
not  generally  agreed  upon.  In  any  case  it  is  of  minor  significance  when  compared  to  the 
effect  of  grass  cover,  so  we  omit  it  in  this  analysis. 

Finally,  in  a  natural  ecosystem  the  population  of  herbivores  fluctuates  in  a  complex  way 
in  response  to  changes  in  the  other  system  variables.  Under  ranching  conditions,  however, 
the  herbivore  population  can  be  considered  as  a  constant.  Stocking  rates  are  set  and 
maintained  by  the  manager  and,  usually,  are  not  altered  until  drastic  changes  (such  as 
death  of  animals)  have  occurred.  The  dynamics  of  the  herbivore  population  can  therefore 
be  omitted,  and  the  analysis  reduces  to  an  examination  of  the  remaining  variables  under 
various  fixed  values  of  herbivore  biomass. 

Hence  the  effects  of  all  five  essential  variables  can  initially  be  reduced  to  the  interactions 
between  only  two  of  them,  i.e.  the  biomass  of  grass  (G)  and  of  woody  vegetation  (W),  and 
their  effects  on  infiltration  rate  (/).  These  variables  will  be  referred  to  henceforward  by 
their  symbols. 

GRAPHICAL  ANALYSIS 

Changes  in  the  relationships  between  G  and  W  are  linked  to  the  soil  water  variables  by 
changes  in  J  (assuming  an  average  annual  rainfall)  and  it  is  therefore  necessary  to  establish 
the  form  of  the  relationship  between  G  and  /.  In  the  absence  of  any  grass  there  will  still  be 
some  infiltration  (i0).  This  minimum  is  likely  to  be  about  10%  of  the  maximum  (Kelly  & 
Walker  1976).  The  rate  of  infiltration  will  initially  increase  rapidly  as  G  increases  and  then 
approach  the  maximum  asymptotically  with  increasing  G,  as  depicted  in  Fig.  1.  For 
convenience,  the  rate  is  expressed  as  a  proportion  of  the  maximum  rate,  and  is  thus 
dimensionless.  It  is  referred  to  as  a  'rate'  for  convenience.  The  value  of  i„  and  rate  of 
increase  will  vary  from  site  to  site,  largely  as  a  function  of  soil  texture.  On  loose,  sandy 
soils  both  the  value  of  i0  and  the  rate  of  increase  in  /  will  be  high,  and  both  will  be  tow  on 
clay  soils.  Unless  otherwise  specified,  the  effects  of  /  in  the  following  graphs  assume  a 
medium  or  heavy  textured  (clay)  soil. 
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Fig.  1.  Hypothesized  relationship  in  semi-arid  savanna  between  rate  of  water  infiltration  into  the 
soil  (/),  expressed  as  a  proportion  of  the  maximum  rate,  and  the  biomass  of  grass  (G). 

Grass  dynamics 

Grass  biomass  (G)  changes  more  quickly  than  does  that  of  woody  vegetation  (W)  and  it 
is  consequendy  analysed  first.  In  the  absence  of  woody  vegetation  and  herbivores  the 
dynamics  of  G  is  controlled  by  available  soil  water.  For  the  moment  we  ignore  the  effects 
of  I  (I0  in  Fig.  2).  When  G  is  small  little  water  is  used  and  the  annual  rate  of  increase  in 
grass  from  one  year  (t)  to  the  next  (t  +  1)  is  large.  As  G  increases,  so  does  water  uptake 
and  the  annual  increase  in  grass  becomes  smaller  until  it  finally  becomes  zero  and,  if  for 
some  reason  G  exceeds  the  biomass  which  can  be  sustained  by  the  available  water, 
negative.  It  is  convenient  to  call  this  the  recruitment  curve  for  G.  The  point  E,  where  R  = 
Gt+l/Gt  =  1,  represents  a  stable  equilibrium  if  all  other  variables  are  fixed. 

Inclusion  of  the  effects  of  infiltration  leads  to  a  depression  in  the  recruitment  curve  at  low 
values  of  G,  as  shown  by  I1  and  I2  in  Fig.  2.  At  higher  values  of  G  infiltration  is  reduced 
very  little  as  G  becomes  smaller,  because  the  litter  cover  remains  high.  Only  when  the 
biomass  of  G  is  small  is  the  effect  significant.  The  reduced  infiltration  leads  then  to  less 
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Fig.  2.  Hypothesized  recruitment  curve  (see  text)  in  semi-arid  savanna  for  grass  biomass  (G). 
The  point  E  represents  a  stable  equilibrium  if  all  other  variables  are  fixed.  /<,  is  the  relationship 
when  /  has  no  effect  on  recruitment.  /,  and  I2  represent  cases  in  which  reduced  infiltration  at  low 

values  of  G  affect  recruitment. 
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erowth  Whether  recruitment  continues  to  rise  as  G  becomes  smaller  (/,  in  Fig.  2),  or 
curves  downwards  (72)  depends  on  whether  or  not  the  reduction  in  the  amount  of  water 
infiltrating  is  more  or  less  than  the  reduction  in  water  uptake  W^  reduced  biomas^o 
Erass  Our  own  observations  of  recovery  rates  on  areas  where  G  has  been  kept  very  small 
throueh  overgrazing  suggest  that  both  conditions  are  possible  and  that  72  is  not  uncommon. 
In otter  wolTe  proportional  recruitment  of  grass  may  be  greater  at  intermediate  values 
of  gfa  WorSs  than  at  low  (or  high)  values.  Because  some  infiltration  will  always  occur 
t  Fig  iTthe  curve  will  never  be  depressed  below  1  at  low  values  of  G,  and  there  remams 

the  single  equilibrium  value  of  G,  at  E.  -— -*• 

wSdy  vegetation  competes  with  grass  for  water  in  the  top-soil,  and  there  is  ample 
evident  from  emi-arid  regions  all  over  the  world  (e.g.  Barnes  1979;  Walker,  Moore  & 
SSSTStT  MTdduAeal  &  Gomez-Gonzalez  1979)  to '^J^^t 
infiltration,  the  effect  of  increasing  Won  the  recruitment  curve  of  G £ ttiat  sh°™  m  F  ^ 
With  small  to  medium  values  of  W  the  effect  is  greater  at  high  values  of  G  than  at  small 
ones  because  when  there  is  little  grass  and  few  trees,  water  uptake  is  less. 

CombmingThe  effects  of  infiltration  and  woody  vegetation  results  in  the  recruitment 
cu^TS  in  Fig.  4.  These  are  qualitatively  different  from  the  individual  effects  because 
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Fig  3  Hypothesized  effect  on  the  recruitment  curve  for  G,  in  semi-arid  savanna,  of  increasing 
amounts  of  woody  vegetation  from  none  {WJ  to  woodland  or  thicket  (Wt). 
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F.G.4.  Hypothesized  combined  effects,  in  semi-arid  savanna,  of  woody  vegetauon  and  of 
Educed  infiltration  at  low  G  (h  of  Fig.  2)  on  the  recru.tment  curve  for  G. 
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Fio.  5.  Equilibrium  conditions  for  grass  biomass  (AG  =  0)  in  relation  to  grass  biomass  (G)  and 
woody  biomass  (W).  The  grass  zero-isocline  is  the  locus  of  points  for  which  G  is  in  equilibrium 
(R  =  Gt+l/G,  —  1-0;  AG  =  0)  at  the  given  value  of  W.  Arrows  show  the  direction  of  movement 
for  small  departures  from  the  isocline.  There  are  two  equilibria  for  W  between  Wa  and  #: 
the  upper  is  stable  and  the  lower  unstable. 


we  assume  that  at  some  large  values  of  W  there  are  now  two  intersections  of  the  recruit- 
ment curve  with  the  G(t+1)/Gt  =  1  line,  giving  the  usual  stable  equilibrium  point  at  high 
values  of  G,  and  also  an  unstable  one  at  low  values  of  G.  This  lower  equilibrium  is  unstable 
because  at  even  smaller  values  of  Gt,  the  recruitment  rate  R  (=  Gt+l/Gt)  is  less  than  1 
and  the  grass  biomass  declines.  At  Gt  values  just  greater  than  the  equilibrium,  however, 
R  >  1  and  G  increases.  The  locus  of  these  equilibria  as  a  function  of  W  can  be  shown  on 
a  graph  of  G  v.  W  (Peterman,  Clark  &  Holling  1979),  as  in  Fig.  5.  We  call  the  locus  of 
equilibrium  points  the  grass  zero-isocline.  It  is  central  in  our  argument. 

Herbivores  affect  grass,  directly,  by  eating  it  and  by  standing  on  it.  Trampling  becomes 
significant  only  at  very  high  levels  of  stocking,  and  is  trivial  in  comparison  to  the  effects  of 
defoliation.  It  is  therefore  not  considered  further.  Light  grazing  has  a  stimulating  effect  on 
grass  growth  (McNaughton  1979),  and  leads  to  a  higher  grass  biomass  than  occurs  when 
the  grass  sward  is  ungrazed  for  any  length  of  time,  as  the  sward  then  becomes  moribund 
and  the  proportion  of  old,  dead  grass  in  the  standing  crop  increases.  Kelly  &  Walker 
(1976)  found  that  under  complete  protection  from  grazing  for  more  than  5  yr  the  above- 
ground  seasonal  grass  production  in  two  successive  years  was  about  66%  and  80%  of  that 
under  very  light  grazing.  As  grazing  becomes  heavier,  however,  grass  production  decreases 
until  under  very  heavy  grazing  the  grass  plants  are  killed.  These  effects  of  increasing 
herbivore  biomass  (//)  on  the  grass  recruitment  curve  are  summarized  in  Fig.  6.  The 
stimulating  effect  of  Hl  does  not  occur  at  low  values  of  G  because,  although  light  grazing 
results  in  only  light  defoliation  if  there  is  a  lot  of  grass  (high  values  of  G),  it  causes  pro- 
portionately increasingly  severe  defoliation  at  lower  values  of  G,  and,  also,  because  the 
build-up  of  large  amounts  of  dead  grass  occurs  only  at  high  values  of  G.  At  low  values 
of  G,  dead  grass  falls  to  the  ground  as  litter,  and  has  a  beneficial  effect. 

It  may  at  first  seem  strange  to  conclude,  as  we  do  in  Fig.  6,  that  larger  amounts  of  grass 
are  apparently  less  able  to  sustain  high  levels  of  stocking  than  are  the  intermediate 
amounts.  However,  the  curves  represent  net  grass  recruitment,  which  is  approximately  the 
product  of  R  for  H0  with  G,  less  the  amount  eaten  (which  is  constant  until  it  approaches 
G).  This  product  is  maximal  at  some  intermediate  value  of  G.  The  depression  of  the  curve 
to  below  the  R  =  1  line  at  large  values  of  H,  once  again  implies  the  existence  of  an  unstable 
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Fig.  6.  Hypothesized  effects  in  semi-arid  savanna  of  increasing  herbivore  biomass  from  none 
(H0)  to  heavily  overstocked  (H«),  on  the  grass  recruitment  curve.  Note  the  order  of  tf0  and  //,  at 

high  values  of  G. 


equilibrium  at  small  values  of  G,  and  a  stable  one  at  larger  values.  Theoretical  and 
empirical  evidence  for  this  unstable  equilibrium  at  small  values  of  G  are  presented  in 
Noy-Meir's  (1975)  lucid  account  of  grazing  systems.  The  zero-isocline  of  G  on  a  graph  of 
G  v.  H  is  shown  in  Fig.  7. 

Finally,  the  combined  effects  of  I,  W  and  H  on  the  equilibrium  positions  (zero-isocline) 
of  G  are  shown  in  Fig.  8  on  the  graph  of  G  v.  W  for  various  values  of  H.  The  effects  of 
infiltration  (J)  are  implicit  in  these  curves. 

Woody  vegetation  dynamics 
Woody  vegetation  dynamics  are  far  simpler  than  those  of  grass  because  herbivores  do 
not  affect  them  significantly,  and  because  the  depressive  effect  of  reduced  infiltration  when 
there  is  little  grass  is  relatively  small.  The  woody  vegetation  has  exclusive  access  to  subsoil 
water,  so  high  values  of  G  cannot  depress  the  recruitment  curve  below  1  (Fig.  9),  and  the 
resulting  zero-isocline  for  W  appears  as  in  Fig.  10,  with  the  dotted  line  indicating  the  effects 
of  reduced  infiltration  at  low  values  of  G.  Woody  vegetation  cannot,  therefore,  be 
eliminated  by  competition  from  grasses  (though  the  reverse  is  possible).  This  conclusion  is 
based  on  the  more  detailed  analysis  by  Walker  &  Noy-Meir  (1981). 


Fio.  7.  The  grass  zero-isocline  in  relation  to  G  and  H,  resulting  from  the  effects  in  Fig.  6.  See 
text  and  Fig.  5  for  further  explanation. 
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Fig.  8.  The  grass  zero-isocline  (zero  recruitment;  R  =  Gt+S/Gt  =  1-0;  AG  =  0)  in  relation  to  G, 
W  and  H.  The  effects  of  infiltration  (I)  are  implicit  in  the  shapes  of  the  curves,  which  summarize 

those  of  Figs  5  and  7. 
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W(  t) 
Fig.  9.  Hypothesized  effect  on  the  recruitment  curve  for  woody  vegetation  (W),  in  semi-arid 
savanna,  of  different  biomass  of  grass  from  woodland  (G0)  to  nearly  open  grassland  (G3). 
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Fig.  10.  Equilibrium  conditions  (zero-isocline)  for  woody  biomass  (AW=  0)  in  relation  to  grass 

(G)  and  woody  {W)  biomass.  The  recurved  dotted  line  indicates  the  effects  of  reduced  infiltration 

at  low  values  of  G,  and  the  dashed  line  {WJ)  is  the  asymptote  (see  text).  This  graph  is  the 

counterpart  of  that  for  grass  in  Fig.  5. 


'"C\( 


070 


482  Stability  of  savanna 

Combined  dynamics 
The  dynamics  of  the  system  as  a  whole  can  now  be  summarized  as  follows.  In  average 
conditions  the  bulk  of  the  rainfall  enters  the  topsoil  and  is  preferentially  taken  up  by  grass 
roots.  Some,  however,  is  taken  up  by  tree  roots  and  when  there  is  heavy  or  continuous  rain, 
or  both,  some  penetrates  to  the  subsoil,  where  it  too  is  taken  up  by  tree  roots.  A  single 
reduction  in  the  amount  of  grass  to  some  very  low  value  does  not  cause  an  immediate 
change  in  the  equilibrium  conditions,  because  of  the  slower  response  of  infiltration 
characteristics  and  (even  slower)  response  of  the  amount  of  woody  vegetation.  If  the  grass 
biomass  is  maintained  at  a  low  value  for  a  sufficiently  long  time  (e.g.  by  several  years  of 
sustained  intense  grazing)  then  the  soil  surface  changes  and  infiltration  rate  declines  with 
the  result  that  less  water  enters  the  soil  but  proportionately  more  of  what  does  enter  may 
penetrate  to  the  subsoil.  The  combined  result  is  that  woody  vegetation  may  have  more 
water  available  to  it,  so  the  biomass  of  woody  vegetation  increases.  These  two 
factors — reduced  infiltration  and  greater  biomass  of  woody  vegetation — combine  to 
prevent  the  re-establishment  of  grasses  even  if  the  grazing  pressure  is  reduced.  Only  if  one 
or  both  of  these  factors  are  altered  to  allow  the  grass  to  develop  to  above  the  critical, 
unstable  equilibrium  (G  in  Fig.  5),  will  the  system  revert  to  its  original  structure.  It  may 
require  more  than  a  single  reduction  in  above-ground  woody  vegetation  to  achieve  this,  as 
the  root  systems  of  savanna  trees  are  remarkably  extensive  and  persistent.  The  competitive 
effect  of  the  woody  plants  must  be  removed  for  long  enough  to  allow  the  grass  root  system 
to  develop  to  the  point  where  it  can  once  again  take  up  more  water  than  the  tree  roots,  and 
thus  allow  the  grass  cover  to  increase  and,  eventually,  infiltration  (to  the  topsoil)  to 

increase. 

This  graphical  analysis  has  proved  useful  in  providing  a  conceptual  understanding  of  the 
dynamics,  but  further  progress  is  not  possible  without  a  more  detailed  quantitative 
understanding  of  the  response  of  the  system  to  changes  in  particular  parameters.  We 
proceed  therefore,  to  a  more  detailed  analysis  of  these  ideas.  The  parameters  introduced 
here  are  summarized  in  Table  1  (p.  490). 

QUANTITATIVE  ANALYSIS 

The  interpretation  of  equations  and  their  parameters  is  greatly  simplified  if  a  consistent 
system  of  units  is  adopted.  In  this  analysis,  the  use  of  []  indicates  the  units,  and  hence  the 
dimensions,  of  the  variable  or  parameter  they  enclose.  This  convention  is  superior,  for  our 
purposes,  to  the  use  of  simple  physical  dimensions  of  mass,  length  and  time  because 
although  the  biomass  of  grasses  and  herbivores  (for  example)  have  the  same  physical 
dimensions  it  is  not  usually  valid  to  add  them.  In  places  where  the  dimensions  are  obvious 
we  omit  specific  definitions. 

Grass  dynamics 
Change  in  the  biomass  of  grass  is  determined  by  five  factors:  the  amount  of  water  in  the 
soil,  the  proportion  of  this  amount  which  is  taken  up  by  grass,  the  efficiency  of  conversion 
of  the  water  that  has  been  taken  up  into  grass  biomass,  the  loss  of  biomass  by  respiration 
and  death,  and  consumption  by  herbivores. 

Amount  of  soil  water 

As  described  previously,  if  we  assume  an  average  annual  rainfall  the  amount  of  soil 
water  is  a  function  of  the  rate  of  infiltration,  /,  which  we  define  in  relative  terms  as  the 
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amount  of  water  entering  the  soil  expressed  as  a  proportion  of  the  amount  which  enters 
when  G  is  at  a  maximum  (Fig.  1).  Therefore,  /  is  dimensionless  and  has  no  units. 

The  form  of  the  relationship  between  /  and  G,  described  by  Fig.  1,  will  be  developed  in 
several  steps.  The  saturation  effect,  modelled  by  a  hyperbolic  relationship,  may  be 
expressed  in  Michaelis-Menten  form  as 

a 
1  = 


G  +  d, 

where  d,  is  the  value  of  G  where  7  reaches  half  its  maximum,  i.e.  7  =  0-5.  However,  we 
require  a  non-zero  infiltration  fraction  when  G  =  0,  so  we  displace  the  origin  of  the  G  axis 
byd2. 

G  +  d, 

/  = 


G  +  d2  +  dj 


When  G  =  0, 7  =  i0  is  given  by 


Then 


d2 


in  = 


d,  +  d, 


/  =  £+cijo<  (1) 

G  +  c, 

Water  uptake 

We  define  a  parameter  UG,  which  is  the  relative  utilization  of  topsoil  water  per  unit  of 
grass  biomass,  such  that  Ua  G  gives  the  utilization  of  water  by  grass,  expressed  as  a 
proportion  of  the  total  amount  of  water  available  (Walker  &  Noy-Meir  1981).  It  is  defined 
by 

U       ?2 ;     [UQ]-l/lGl  (2) 

G      8G.G  +  6w.W+e 

where  8Q  and  0W  are  the  rate  of  water  uptake  (as  a  depth  of  water  per  unit  time)  per  unit 
biomass  of  grass  and  woody  vegetation  respectively,  and  e  is  the  loss  of  water  (as  depth  per 
unit  time)  by  evaporation.  If  /  is  the  rainfall  in  a  given  time,  T,  then  [6G]  =  [/]/[G][T]; 
[0W]  =  l/]/[*FHT]  and  [e]  =  Ul/lT].  The  evaporation  loss  is  generally  insignificant  and  has 
little  effect  on  the  dynamics  of  the  system,  but  is  included  for  the  sake  of  completeness. 
Ua  is  therefore  a  term  which  reflects  the  competitive  effect  of  woody  vegetation  in  reducing 
the  amount  of  water  taken  up  by  grasses. 

Conversion  efficiency  of  water  into  grass  biomass 

This  can  be  considered  as  the  growth  rate  (productivity)  of  grass,  denoted  by  rG;  [rG]  = 
[G]/IT].  It  is  a  measure  of  the  increase  in  the  biomass  of  grass  during  a  given  time  which 
results  from  the  total  amount  of  water  taken  up. 

Losses  by  respiration  and  death 

Let  LG  denote  the  proportional  rate  of  loss  of  grass  biomass;  [LG1  =  1/[T).  Then  G  LG  is 
a  measure  of  the  grass  biomass  lost  per  unit  time. 

. ,  r\  n  7  ° 
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Consumption  by  herbivores 

The  stimulating  effect  on  grass  growth,  described  by  Fig.  7,  is  not  included  here  for  two 
reasons.  First,  it  affects  only  the  upper  parts  of  the  grass  isocline  (when  the  biomass  of 
grasses  is  large),  and  therefore  has  no  effect  on  the  qualitative  behaviour  of  the  system, 
which  is  determined  by  the  change  from  stable  to  unstable  equilibrium  when  grass  biomass 
is  small  Secondly,  as  already  described,  it  is  difficult  to  quantify  the  stimulating  effect  in 
terms  of  individual  plant  effects  (as  opposed  to  the  accumulated  effects  on  the  sward  as 
a  whole). 

Losses  due  to  grazing  must  incorporate  the  two  features  shown  in  Fig.  1 1 ,  where  the 
relation  between  total  consumption,  C,  at  a  fixed  level  of  grazing,  and  the  grass  biomass,  G, 
is  shown.  The  consumption  asymptotically  approaches  a  satiation  level  as  G  increases  (cf. 
Morley  1966;  Noy-Meir  1975).  Below  some  minimum  value  of  G  no  grazing  occurs,  either 
because  it  is  impossible  (the  animals  cannot  get  at  the  grass)  or  because  food  quality  or 
some  other  feature  of  the  grass-herbivore  interaction  renders  the  grass  unacceptable  to  the 
herbivores.  When  considered  over  a  short  time  (a  few  seasons),  recruitment  of  grass  is  a 
function  of  total  biomass,  above  and  below  ground.  Over  50%  of  the  biomass  of  many 
savanna  grasses  may  be  below  ground  (e.g.  Rushworth  1975;  Kelly  &  Walker  1976).  Even 
if  all  of  the  above-ground  biomass  is  removed,  re-growth  will  occur  the  following  season. 
Defoliation  at  this  intensity,  if  continued,  will  kill  the  plants,  but  in  the  short  term  the 
stolons  and  rhizomes  supply  a  refuge.  A  further  mechanism  that  allows  some  grass  to 
survive  above  ground  is  the  development  of  a  creeping  growth  habit.  It  is  a  commonly 
observed  phenomenon  that  grasses  with  this  ability  (e.g.  Cynodon  spp.*)  tend  to  take  over 
from  erect  bunch  grasses  (e.g.  species  of  Brachiaria,  Setaria  and  Themeda)  under  heavy 
grazing,  and  McNaughton  (1980),  working  in  the  Serengeti,  describes  how  some  grasses 
change  their  growth  habit  under  different  grazing  pressures. 

There  are  a  variety  of  reasons  why  herbivores  will  not  eat  certain  kinds  or  quantities  of 
grasses.  Wilson  &  Graetz  (1979)  described  an  Australian  grassland  where  the  nutritive 
value  of  all  but  a  small  percentage  of  the  plant  parts  is  too  low  to  be  of  any  use  to  domestic 
livestock.  Similar  types  (e.g.  Eragrostis  pallens)  occur  in  southern  Africa.  Some  grasses  are 
unpalatable  owing  to  taste,  e.g.  species  of  the  genera  Cymbopogon,  Bothriochloa, 
Trachypogon,  Elyonurus  and  Urelytrum,  and  others  because  they  are  too  tough  and 
fibrous,  e.g.  Eleusine  indica,  Aristida  junciformis,  some  Sporobolus  spp.,  and  so  on. 


Fig.  11.  Hypothesized  relation  between  tota!  consumption  (C)  of  grass  on  semi-arid  savanna  at 
a  fixed  level  of  herbivore  biomass,  and  the  biomass  of  grass  ((?).  Below  GQ  no  grazing  is 

possible. 


'  Nomenclature  follows  Simon  (1971). 
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The  combined  effect  of  these  different  mechanisms  is  the  grazing  cut-off,  represented  by 
G0  in  Figs  11  and  12.  The  value  of  G0  can  vary  considerably. 

We  may  incorporate  the  saturation  and  cut-off  effects  by  following  the  procedure  which 
led  to  Eqn  (1).  The  rate  of  consumption  of  grass,  Q. ,  can  be  described  by 


/c2//(G-G0) 

Q  =  max ,  0 

G  +  c, 


(3) 


where  c2  =  rate  of  consumption  of  grass  biomass  per  unit  of  herbivore  biomass;  [c2]  = 
[G]/([H][7]);  and  d3  =  the  value  of  G  at  which  the  consumption  rate  is  half  maximum,  and 
c3  =  d3  —  G0. 

Now  we  have  terms  which  correspond  to  each  of  the  factors  which  affects  the  change  in 
the  biomass  of  grass.  The  combined  equation  is 


—  =  rGJC/GG-LGG-//eG. 

at 


Here  we  have  set 


UG  = 


I 


where  a=6v/8G;P  =  e/0G ;  and 


2  =  max 


G  +  aW  +  0 
c2(l -G0/G) 


(4) 


(5) 


G  +  c, 


,0 


The  zero  isocline  for  G  (that  is  the  equilibrium  line  shown  in  Fig.  5)  is  the  curve  where 
dG/dt  =  0.  From  Eqn  (4),  we  see  that  this  curve  has  two  branches:  G  =  0  (a  trivial  result) 
and 


rGJ 


LG  +  OH 


G+aW+p 
After  substitution  of  Eqn  (1),  this  relation  may  be  solved  for  W: 


W  =  -\-p-G  +  - ~ 


1- 


c,(l-  i0) 


(6) 


G  +  Cj 
The  intercept  of  this  isocline  with  the  W  axis  is  obtained  by  setting  G  =  0  in  Eqn  (6)  and  4s 


wBJ-[r^-PY 


(7) 


The  value  of  W0  is  independent  of  if  (provided  G0  >  0)  because  of  the  grazing  cut-off.  The 
slope  of  the  isocline  may  be  obtained  by  differentiating  Eqn  (6).  If  G  is  large,  it  can  be 
shown  that  this  slope  is  negative.  Where  G  =  0,  we  have 

rG(l-io) 


qW 
dG 


0=0 


-1  + 


Loci 


(8) 


This  expression  will  generally  be  positive.  Thus  Eqn  (6)  when  differentiated  will  have  a 
root:  the  corresponding  point  {W,  G;  Fig.  5)  is  a  'turning  point'  for  the  isocline.  The 
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unstable  section  at  low  values  of  G  is  due  to  reduced  infiltration  when  the  biomass  of  grass 

is  small. 

Now  we  consider  the  effect  of  the  herbivore  biomass  upon  the  isocline.  If  G  >  G0,  then 
Q  >  0.  As  H  in  Eqn  (6)  increases,  the  factor  rG/(Lu  +  QH)  decreases.  The  quantity  in 
brackets  which  multiplies  it  is  positive.  Hence  the  isocline  moves  to  the  left  as  H  increases. 
Finally,  for  very  large  values  of  H,  the  isocline  will  coincide  with  a  segment  where  G  -  G0. 
If  G  <  G0,  the  isocline  does  not  depend  on  H.  The  net  effect  is  shown  in  Fig.  12. 


Woody  vegetation  dynamics 

The  equation  for  the  change  in  woody  vegetation  will  be  developed  in  a  similar  way  to 
that  for  grass,  except  that  the  woody  vegetation  has  access  to  subsoil  water,  which  is  not 
available  to  the  grass.  The  Michaelis-Menten  form  of  the  hyperbola  is  adequate:  the  extra 
term  has  the  form  reW/(W  +  p).  Here  r5  =  increase  in  biomass  of  woody  vegetation,  during 
a  given  time,  for  the  subsoil  water  taken  up,  [rs]  =  [W]/[T];  and  p  =  biomass  of  woody 
vegetation  which  might  have  been  formed  in  response  to  water  lost  by  percolation  beyond 
the  roots  in  the  subsoil  (i.e.  a  measure  of  the  inefficiency  of  subsoil  water  use;  [p]  =  [  W\. 

There  is  one  further  modification.  If  G  decreases,  the  total  amount  of  water  entering  the 
soil  decreases  but,  as  there  is  very  little  grass  to  take  it  up,  there  may  be  more  water  passing 
through  to  the  subsoil  (depending  on  soil  type  and  the  value  of  1$).  This  means  that  rs  is  not 
constant,  but  is  a  function  of  G,  in  the  opposite  sense  to  I.  We  therefore  set 

S  =  rs[l-G/(K(G  +  d«))] 

where  rs  =  the  increase  in  biomass  of  woody  vegetation  at  G  =  0,  during  a  given  time,  for 
the  subsoil  water  taken  up;  K  determines  the  proportional  reduction  from  rs  to  the  lower 
asymptotic  value  of  S  =  rs(l  -  1/K)  at  high  values  of  G;  and  d„  =  the  half  saturation  value 
of  G.  Usually,  K  will  be  greater  than  or  at  least  equal  to  1,  but  if  i0  is  very  low,  it  could  drop 
below  1. 

By  analogy  with  Eqn  (5),  we  set 
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Fig.  12.  Equilibrium  conditions  (AG  =  0)  for  grass  biomass  in  relation  to  grass  biomass  (G)  and 
woody  biomass  (W)  at  various  grazing  levels  (herbivore  biomass.  H).  The  cut-off  value  G0  is  the 
biomass  of  grass  below  which  stock  do  not  eat  the  grass.  See  text  and  Fig.  5  for  further 

explanation. 
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where  y  =  0G/0W  =  Ma  and  S  =  e/6>w  =  fila.  Dimensions  are  defined  by  |t/wl  =  1/1  W); 
\y]  =  [W]/lGl;and  \S]  =  \  W\.  Thus  the  dynamic  equation  for  W  has  the  form 

dW 

=  rw/l/w  W+SW/(W+p)-LvW  (10) 

dt 

where  rw  and  Lw  (analogous  to  rG  and  LG)  represent  the  growth  rate  due  to  topsoil  water 
uptake,  and  the  proportional  rate  of  loss  (by  respiration  and  death),  respectively.  A 
consumption  term  is  omitted  from  Eqn  (10),  because  virtually  no  woody  vegetation  is 
eaten  by  stock. 

The  zero  isocline  for  W  (the  equilibrium  line  shown  in  Fig.  10)  is  obtained  by  setting 
dW/dt  =  0.  The  two  branches  are  given  by  W=0,  and 

rw/£/w  + =  LW.  (11) 

w      w      W  +  p 

The  intercept  Wy  of  this  isocline  with  the  Waxis  is  obtained  by  setting  G  =  0  in  Eqn  (11). 

If  G  is  large,  C/w  is  small  and  the  first  term  in  Eqn  (1 1)  is  negligible.  Thus  the  isocline  has 
a  vertical  asymptote:  W=  Wa.  From  Eqn  (1 1)  it  follows  that 

Lw 

Now  we  turn  to  intermediate  values  of  G.  As  a  first  step,  the  effect  of  G  upon  infiltration 
may  be  ignored;  thus  we  set  J  =  1  in  Eqn  (1 1).  The  equation  of  the  isocline  then  takes  the 
approximate  form 

=  LW.  (13) 


W+yG  +  6      W+p 

If  G  increases,  the  first  term  in  Eqn  (13)  decreases.  The  balance  can  only  be  made  up  by  a 
decrease  in  W.  Therefore  the  biomass  of  woody  vegetation  decreases  as  that  of  grass 
increases.  This  relationship  may  not  hold  when  there  is  very  little  grass,  because  of  the 
effect  of  grass  upon  infiltration.  The  effect  upon  the  W  isocline  will  depend  upon  whether 
the  benefit  to  woody  vegetation  of  a  reduced  uptake  of  water  by  grass  is  offset  by  the 
reduced  infiltration.  The  general  form  of  the  isocline  is  shown  in  Fig.  10. 

Combined  dynamics 

When  the  zero  isoclines  (equilibrium  lines)  for  G  and  W  are  combined,  the  result  can  be 
as  shown  in  Fig.  13.  The  details  may  change  as  parameter  values  are  changed.  There  may 
be  two  stable  equilibria  in  the  absence  of  herbivores.  The  equilibrium  point  with  large 
biomass  of  grass  and  small  biomass  of  woody  vegetation  may  be  regarded  as  the  usual 
condition.  The  two  stable  equilibria  are  separated  by  a  saddle  point,  which  is  an  unstable 
equilibrium.  The  whole  area  of  the  grass  biomass-woody  vegetation  biomass  graph  may  be 
divided  into  two  parts,  each  containing  a  stable  point.  From  anywhere  within  one  area  the 
vegetation  will  tend  to  move  toward  the  stable  point.  The  line  defining  these  areas  is  the 
separatrix  (Fig.  13),  and  its  position  is  determined  by  the  relative  speeds  with  which  G 
and  JF  change,  not  merely  by  the  the  positions  of  the  isoclines. 

Varying  the  parameters  in  the  model,  in  various  combinations,  leads  to  two  main 
changes  in  the  relative  positions  of  the  isoclines.  First,  the  W  zero-isocline  moves  to  the 
right  if  5  (the  rate  of  growth  of  W  due  to  subsoil  water)  increases  or  if  Lw  (metabolic  loss 
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w0 

Fig.  13.  Probable  relative  positions  of  the  C  and  W  zero-isoclines  Oocus  of  points  of 
equilibrium  where,  for  X  =  G  or  fY,  X^JX%  =  I  -0;  AX  —  0)  in  relation  to  biomass  of  grass  and 
of  woody  vegetation.  There  are  two  stable  and  two  unstable  equilibria.  The  stable  equilibria  are 
separated  by  an  unstable  equilibrium,  or  saddle  point.  The  dashed  lines  show  the  directions  of 
movement  of  a  system  not  at  equilibrium.  The  interrupted  line  (separatrix)  defines  areas  from 
within  which  a  system  would  move  to  one  or  other  of  the  stable  positions. 


rate  of  W)  or  p  (percolation  loss)  decrease  (Fig.  14).  Eventually  the  two  intersections  of  the 
xero-isoclines  are  lost,  and  woody  vegetation  dominates.  In  the  reverse  situation,  if  S 
decreases  or  if  Lw  or  p  increase,  the  W  zero-isocline  moves  to  the  left,  with  two  possible 


Fig.  14.  Effect  on  the  woody  vegetation  zero-isocline  (AW  =  0)  of  increasing  5  (the  rate  of 

growth  in  response  to  subsoil  water)  or  of  decreasing  Lw  (respiratory  loss)  or  p  (percolation  loss). 

The  grass  zero-isocline  is  in  the  same  position  as  it  is  in  Fig.  1 3,  and  does  not  intersect  the  A  If  =  0 

line.  All  mixtures  move  toward  domination  by  woody  vegetation. 
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Fig.  15.  As  in  Fig.  1 4  but  decreasing  S  or  increasing  L,,  or  p.  The  A  W  =  0  line  moves  to  the  left. 
In  (a)  there  are  two  stable  points.  In  (b),  where  the  zero-isocline  has  moved  even  further  left, 

grassland  alone  is  stable. 


Fig.  16.  Effect  on  the  grass  zero-isocline  (AG  =  0)  of  increasing  i0,  rc  or  G„,  or  of  decreasing  cr, 

LG  or  H  (see  Table  1  for  definition  of  symbols).  The  woody  vegetation  zero-isocline  is  in  the 

same  position  as  it  is  in  Fig.  1 3.  All  mixtures  move  towards  a  single  stable  mixture. 

outcomes.  In  Fig.  15(a)  the  upper  intersection  is  lost  and  the  result  is  that  either  grass  or 
woods  vegetation  dominate,  depending  on  which  side  of  the  separatrix  the  system  starts 
from.  In  Fig.  15(b)  both  intersections  are  lost  and  grass  comes  to  dominate.  Secondly,  the 
G  zero-isocline  may  be  moved:  to  the  right  if  i0  or  rG  or  GQ  increases,  or  if  or  or  LG  or  H 
decrease.  The  lower  intersection  eventually  disappears,  leading  to  a  single  stable 
equilibrium  and  co-existence  of  grass  and  woody  vegetation  (Fig.  16).  If  the  parameters 
change  in  the  opposite  direction,  the  G  zero-isocline  moves  to  the  left  until  eventually  the 
same  situation  occurs  as  in  Fig.  14,  and  grass  is  eliminated. 

Example  from  an  'average'  savanna 
As  a  final  step  in  our  analysis  it  was  necessary  to  establish  whether  the  basic  situation  as 
depicted  in  Fig.  13  can,  or  actually  does,  occur.  To  do  this  we  need  to  quantify  the 
variables  and  parameters  of  the  model,  and  to  test  the  effects  of  changing  them.  A  complete 
set  of  values  does  not  exist  for  any  one  site,  but  an  'average'  semi-arid  savanna  was 
synthesized  from  various  sources  (Kelly  1973;  Rutherford  1979;  Huntley  &  Morris  1978; 
Walker  &  Noy-Meir  1981).  Where  no  data  were  available,  such  as  for  the  partition  of  total 
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Table    1.  The  parameters  used  in  the  model  of  grass  and  woody  vegetation 
dynamics.  The  values  in  parentheses  are  those  used  for  the  standard  run  of  the 

model,  shown  in  Fig.  17. 

Parameters  for  grass: 

rG         rate  of  growth  due  to  topsoil  water  (4000  kg  ha" '  yr" ' ) 

Cj  rate  of  consumption  of  grass  by  herbivores  (315  kg  kg"'  yr';at  H=  lOOkgha"') 

G0        the  grazing  threshold  (500  kg  ha"') 

d3  half-saturation  level  for  grazing  (600  kg  ha"1) 

i0  minimum  infiltration  at  low  values  of  G  (0- 1) 

c,         the  rate  at  which  infiltration  increases  with  G 

c,  =  d,  +  dj;     d,  =  500  kg  ha"1;     d2  is  derived  from  other  parameters,  see  p.  483. 
a  the  quotient  of  the  efficiency  of  uptake  of  water  from  topsoil  by  woody  vegetation  to  that  by  grasses 

(=eje0\6Q  =  0-21  mm  kg"1  ha"'  yr"1;  9m  =  0- 122  mm  kg"1  ha"1  yr1) 
LG        the  proportional  rate  of  loss  of  grass  by  respiration  and  death  (1-0  yr"1) 

Parameters  for  woody  vegetation: 
r. 


rate  of  growth  due  to  topsoil  water  (3000  kg  ha"1  yr"1) 

rate  of  increase  of  biomass  of  woody  vegetation  due  to  subsoil  water  taken 

upatG  =  0(4000  kg  ha"1  yr"1) 
dt         half-saturation  value  of  G  for  the  decrease  in  S  with  increasing  G  (500  kg  ha-1) 
K         proportional  reduction  from  rs  to  lower  asymptote  of  S  at  high  values  of  G  (2-0) 
p  biomass  of  woody  vegetation  which  might  have  been  formed  by  water  lost  through  the 

subsoil  by  percolation  (inefficiency  of  sub-soil  water  use)  (0- 1  kg  ha"1) 
Lw        the  proportional  rate  of  loss  of  woody  vegetation  by  respiration  and  death  (1-0  yr"') 


woody  vegetation  growth  into  that  due  to  topsoil-water  and  that  due  to  subsoil-water,  (rw 
and  S),  the  estimated  minimum  and  maximum  values  were  tested.  In  the  absence  of  better 
information,  losses  from  respiration  and  death  (LG  and  Lw)  were  made  equal,  and  subsoil 
water  loss  (p)  and  evaporation  loss  (e)  were  kept  at  very  low  values.  The  standard  set  of 
values  used  for  what  was  considered  to  be  an  average  semi-arid  savanna  with  soils  prone  to 
capping,  under  light  grazing,  are  shown  in  Table  1.  The  position  of  the  separatrix  was 
obtained  by  plotting  the  trajectories  of  G  and  W,  starting  at  various  positions  along  the  W 
axis  and  increasing  G  until  the  direction  of  the  trajectory  changed.  The  results  of  the 
standard  run  are  given  in  Fig.  17. 

With  the  parameter  values  given  in  Table  1,  the  multiple  equilibrium  situation  does  arise, 
and  changes  in  the  parameters  lead  to  the  expected  relative  changes  in  the  position  of  the 


All  units 
(kg  ho"1) 


2000     3000      4000   W 
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Fig.  17.  Zero-isoclines  for  woody  vegetation  (AW  =  0)  and  for  grass  (AG  =  0)  at  different 

herbivore  biomass  (H)  using  parameter  values  given  in  Table  1.  In  (a)  with  grazing  threshold 

500  kg  ha"1,  the  upper  equilibrium  exists  for  all  values  of  H.  In  (b)  with  grazing  threshold 

10  kg  ha"1,  grass  is  eliminated  when  H  is  slightly  less  than  100  kg  ha"1. 
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grass  zero-isocline.  A  number  of  interesting  consequences  emerge  from  the  results  obtained 
by  varying  the  parameter  values,  as  described  in  the  next  section.  At  this  stage,  one 
important  result  is  worth  noting.  With  a  value  of  i0  (infiltration  minimum)  greater  than  0-5 
and  with  a  low  value  of  K  (proportional  increase  in  subsoil  water  as  grass  biomass 
decreases)  there  are  only  two  equilibria:  the  single  stable  joint  equilibrium  (with  grass  and 
woody  vegetation)  and  woody  vegetation  alone.  Because  i0  has  such  a  powerful  effect,  an 
unstable  threshold  does  not  occur  on  sandy  soils  (for  which  i0  is  large)  and  the  system  will 
always  be  able  to  return  to  a  joint  equilibrium  (provided  grass  is  not  totally  eliminated). 
Furthermore,  if  K  <  0  then  even  at  very  low  values  of  i0  the  unstable  joint  equilibrium  either 
does  not  exist  or  it  occurs  very  close  to  the  If-axis,  virtually  at  G  =  0.  This  highlights  the 
importance  of  changes  in  the  absolute  amount  of  water  passing  through  to  the  subsoil  as 
the  biomass  of  grass  decreases.  There  are,  as  yet,  no  empirical  tests  of  this  prediction,  but 
two  studies  have  been  initiated  in  South  African  savanna  to  obtain  estimates  of  K  and  d4. 


EFFECTS  OF  GRAZING  AND  MANAGEMENT 

We  return  now  to  the  theme  of  the  introduction,  using  the  stability  analysis  to  suggest 
explanations  for  the  observed  impact  of  different  patterns  of  grazing  on  semi-arid 
savannas.  We  describe  four  different  patterns  emerging  from  inferred  differences  in  the 
form  and  relative  positions  of  the  isoclines.  These  four  examples  are  then  used  to  define  the 
responses  of  the  grass-wood  vegetation  systems  to  disturbance. 

Grazing  prior  to  1800 

Prior  to  the  main  impact  of  European  settlement,  the  savannas  of  southern  and  eastern 
Africa  were  grazed  by  many  species  of  native  ungulate  and  by  the  cattle  of  semi-nomadic 
herdsmen.  Human  population  density  was  low,  especially  in  the  drier  regions  where,  apart 
from  the  permanently  watered  areas  along  rivers,  grazing  by  cattle  was  irregular  and 
opportunistic,  depending  on  the  rainfall.  Because  of  the  migratory  behaviour  of  the 
native  ungulates  and  the  movements  of  the  herdsmen  any  one  area  was  grazed  fairly  lightly 
for  most  of  the  time,  but  periodically  it  was  grazed  heavily,  followed  by  a  period  of 
recovery.  A  number  of  perennial  grasses,  including  many  early  successional  species  (here- 
after referred  to  as  'resistant'  grasses),  can  withstand  intense  grazing  because  a  significant 
proportion  of  their  biomass  is  underground.  In  addition,  some  develop  a  recumbent  habit 
when  grazed  and  others  are  unpalatable  and  of  low  value  as  food.  At  the  other  end  of 
the  spectrum  are  the  palatable,  leafy  perennials,  regarded  as  desirable  by  range  managers. 
These  grasses,  while  more  productive,  are  more  'sensitive'  because  they  lack  the  under- 
ground refuge  or  recumbent  habit  and  are  the  first  to  be  grazed.  Hence  the  particular 
balance  achieved  between  the  two  types  is  the  result  of  the  shifting  dynamics  of  competition 
between  them.  During  periods  of  intense  grazing,  resistant  grasses  are  favoured;  during 
periods  of  low  or  moderate  grazing,  the  sensitive  grasses  are  favoured. 

The  resulting  zero-isoclines  for  pre- 1900  are  suggested  in  Fig.  18(a),  together  with  a 
typical  trajectory.  Two  equilibrium  conditions  are  possible — one  allowing  coexistence  of 
grasses  and  woody  vegetation  and  the  other  in  which  woody  vegetation  predominates.  On 
the  one  hand,  overall  grass  productivity  is  fairly  high  because  the  productive  sensitive 
grasses  are  maintained  in  the  system.  On  the  other,  persistence  at  the  upper  equilibrium, 
with  grasses  dominant,  is  encouraged  because  of  the  maintenance  of  resistant  grasses 
(those  with  a  high  C0).  In  these  conditions,  the  system  is  insensitive  to  even  unusually  large 
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Fig,  18.  Zero-isoclines,  as  in  Fig.  17,  for  different  herbivore  biomass  (H)  in  four  situations. 
Large  filled  circles  on  the  W  axis  and  thick  lines  represent  stable  equilibria  with  grass  and  woody 
vegetation  coexisting.  Unfilled  circles  indicate  unstable  equilibria.  Arrowed  ellipses  represent  the 
normal  variability  in  G  and  W  under  the  different  forms  of  land-use:  (a)  southern  Africa, 
pre- 1800:  (b)  southern  Africa,  settled  peasant  fanners,  stock  grazing:  (c)  southern  Africa, 
ranching,  fixed  moderate  level  of  stocking:  (d)  south-western  U.S.A.,  mesquite  grassland. 

changes  of  the  variables.  In  this  simplified  caricature,  herbivore  density  could  even  rise  to 
infinity,  and  the  available  grass  biomass  could  go  to  zero,  without  the  system  moving  to  the 
lower  equilibrium.  In  addition,  the  distance  between  the  upper  set  of  equilibrium  points 
(Fig.  18(a),  A  and  B)  and  the  separatrix  that  runs  through  the  saddle  point  (D)  is  large 
enough  that  it  is  difficult  to  imagine  a  mechanism  that  would  move  shrub  density  into  the 
lower  equilibrium  position.  Although  this  argument  is  appropriate  for  this  simplified 
caricature,  the  reality  would  not  be  as  extreme.  In  particular,  a  persistent,  rather  than 
transient,  large  increase  of  H  or  decrease  of  G  would  lead  to  elimination  of  sensitive  grasses 
and  the  movement  of  the  G  isocline  to  the  left,  perhaps  to  the  point  where  the  separatrix 
was  crossed.  Drought  could  produce  the  same  result. 

Settled  peasant  farmers,  stock  grazing 

Those  grazing  areas  managed  by  peasant  farmers  in  southern  Africa  are  now 
characterized  by  higher  human  density,  fixed  settlements,  and  intense  and  persistent  (i.e. 
unvarying)  grazing  by  cattle.  This  persistent  and  intense  grazing  results  in  areas  with 
significant  amounts  of  shrubs,  little  available  grass,  and  a  high  density  of  emaciated  cattle. 
Although  such  areas  appear  to  be  on  the  edge  of  collapse  (e.g.  Anon.  1920;  West  1947, 
1967;  Allan  1965;  Cleghorn  1966),  they  typically  persist  on  that  side  of  the  separatrix 
where  both  grass  and  woody  vegetation  coexist. 

These  effects  on  the  isoclines  are  suggested  in  Fig.  18(b).  Under  such  prolonged,  heavy 
grazing,  the  biomass  of  productive  grasses  is  greatly  reduced,  water  run-off  is  greater,  and 
erosion  leads  to  a  decline  in  soil  productivity.  The  result  is  that  both  rG  and  rw  (the  growth 
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rates  of  grass  and  woody  vegetation)  are  reduced,  leading  to  a  shift  of  both  the  grass  and 
shrub  isoclines  to  the  left.  Moreover,  the  demand  for  firewood  can  shift  the  shrub  isocline 
even  further  left  than  that  shown,  potentially  eliminating  the  lower  equilibrium  where 
shrubs  can  predominate. 

Hence  the  reduced  variability  and  the  high  intensity  of  grazing  causes  parameters  of  the 
system  to  evolve,  producing  a  much  less  productive  savanna,  but  one  which  still  commonly 
persists  at  an  upper  equilibrium  point.  The  biomass  of  ungrazeable  grasses,  G0,  is 
maintained  at  a  high  value  as  in  the  previous  pattern,  and  the  area  on  the  G  v.  W  graph 
which  is  occupied  by  the  lower  equilibrium  could  well  be  reduced.  Although  this  system, 
like  the  previous  one,  is  insensitive  to  large  changes  in  value  of  H  and  G,  productivity  could 
decrease  so  much  that  it  lowers  human  carrying  capacity,  leading  to  a  human  density  too 
low  to  effectively  reduce  W.  Drought  would  further  contribute  to  such  a  collapse.  Hence 
drought-induced  disasters,  even  though  uncommon,  are  more  likely  to  occur  than  in  the 
previous  system. 

Ranching,  with  fixed  moderate  stocking  level 
The  types  of  grazing  management  on  the  larger  commercial  ranches  feature  either 
moderate  but  persistent  levels  of  stocking,  or  seasonal  rotational  grazing  and  provision  of 
drinking  water  supplies  for  cattle.  Both  the  variability  and  intensity  of  grazing  are 
considerably  less  than  in  the  peasant  farmer  system.  As  a  consequence,  the  balance  of 
vegetation  composition  shifts  towards  highly  productive  sensitive  grasses  so  that  the  land 
has  higher  rates  of  meat  production  than  either  of  the  two  previous  cases.  Experienced 
range-managers  note  that  the  effects  of  sudden  management  changes,  especially  if  coupled 
with  a  drought,  can  be  far  more  devastating  on  a  well-managed  ranch  than  on  areas  which 
have  been  'misused'  for  some  time.  It  is  a  system  of  high  productivity  but  one  that  is 
sensitive  to  change. 

Figure  18(c)  summarizes  these  effects  on  the  zero-isoclines.  That  of  woody  vegetation  is 
little  affected  but  the  grass  zero-isocline  shifts  to  the  right  because  the  rate  of  grass  growth 
(rG)  is  increased.  More  important,  the  ungrazeable  grass  refuge  (G0)  is  greatly  reduced,  as 
the  sensitive  grasses  become  dominant,  so  that  grazing  intensities  that  previously  could  be 
supported  (without  crossing  the  separatrix)  no  longer  can  be.  Under  reduced  but  persistent 
grazing,  G0  evolves  to  bring  the  separatrix  close  to  the  region  in  the  G  v.  W  graph  which  is 
occupied  by  the  system.  Hence  small  increases  in  stocking  rate  or  in  exogenous  events  (e.g. 
drought)  that  move  the  grass  isocline  leftward,  can  easily  lead  to  the  system  moving  to  the 
lower  equilibrium  and  a  collapse  of  the  productive  cattle  grazing  system. 

Mesquite-grass  systems  in  southwestern  U.SA. 

Prior  to  the  high  density  of  cattle  stocking  in  the  late  1800s,  the  ranges  of  the 
southwestern  United  States  had  very  little,  if  any  woody  vegetation  on  the  uplands. 
Mesquite  (Prosopis  glandulosa  Torr.  var.  glandulosa)  was  confined  to  drainage  lines 
(Buffington  &  Herbel  1965).  Subsequently,  intense  grazing  by  cattle  resulted  in  invasion 
and  predominance  of  mesquite  in  the  uplands,  with  much  less  grass.  Once  grass  had  gone, 
an  erosion  pavement  was  formed  that  reduced  infiltration  of  water  into  the  soil,  which, 
in  turn,  led  to  a  high  susceptibility  to  erosion  (Glendenning  1952).  Moreover,  removal  of 
cattle  at  this  point  was  not  followed  by  recovery  of  the  range  (Buffington  &  Herbel  1965). 

Three  conditions  distinguish  this  system  from  those  in  south-eastern  Africa.  First, 
mesquite  roots  apparently  do  not  have  access  to  a  reservoir  of  water  that  is  unavailable  to 
grass.  Competition  for  water  between  grass  and  shrubs  is  therefore  direct.  Secondly, 
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establishment  of  seedlings  of  mesquite  is  markedly  curtailed  by  perennial  grass  cover.  The 
number  of  summer-planted  seeds  required  to  produce  one  established  seedling  in  the 
autumn  is  16  times  higher  on  grassed  than  ungrassed  sites,  and  subsequent  losses  are  8 
times  higher  (Glendenning  1952).  Both  these  conditions  would  greatly  depress  the  shrub 
zero-isocline  at  low  H^and  high  G  (Fig.  18(d)).  In  addition,  the  direct  competition  for  water 
between  grass  and  mesquite  would  move  the  W  zero-isocline  to  the  right  and  the  grass 
zero-isocline  to  the  left  when  shrubs  were  abundant.  Thirdly,  seed  dispersal  of  mesquite  is 
greatly  facilitated  by  cattle.  High  stocking  levels  result  in  increased  rodent  populations,  and 
Glendenning  &  Paulsen  (1955)  have  recorded  the  resulting  increase  in  seed  dispersal  both 
by  rodents  (through  lost  caches)  and  cattle  (by  the  spreading  of  seeds  in  dung).  Although 
not  shown  in  Fig.  19(d),  this  would  result  in  the  W  zero-isocline  shifting  even  further  to  the 
right  as  grazing  intensity  increased.  The  W  zero-isocline  shown  in  the  figure  represents  that 
likely  at  intermediate  to  high  grazing  intensity. 

These  changes  in  isoclines  result,  in  only  two  stable  equilibrium  conditions-with  either 
grass  or  mesquite  predominating  (Fig.  18(d)).  Stable  coexistence  of  both  is  not  possible.  So 
long  as  grazing  is  of  moderate  intensity,  the  range  persists  as  a  grassland.  At  a  critical  point 
however,  further  increase  in  grazing  intensity  leads  to  the  vegetation  moving  towards 
dominance  by  mesquite.  The  grass  sward  does  not  recover  if  cattle  are  removed. 

DISCUSSION 

This  analysis  of  semi-arid  savanna  grazing  systems  has  largely  concerned  qualitative 
properties  of  stability.  It  was  developed  from  knowledge  of  the  existence  of  key  processes 
and  their  qualitative  form.  The  parameter  values  are  so  approximate  that  only  a  rough  idea 
of  magnitude  is  possible.  Moreover,  the  model  was  confined  to  the  dynamics  of  only  two 
variables — grass  and  woody  shrubs — under  different  levels  of  grazing.  Other  dynamic 
variables  such  as  herbivore  populations  and  the  proportion  of  sensitive  and  of  resistant 
grasses  are  likely  to  be  important,  but  their  dynamics  were  inferred  rather  than  explicitly 
represented.  Nevertheless,  this  simplified  representation  identifies  some  causes  and 
interrelations  between  three  important  properties  of  ecological  systems — productivity, 
resilience  and  variability.  The  relative  values  of  these  properties  in  the  three  grazing 
systems  of  southern  Africa  are  summarized  in  Table  2. 

The  differences  in  herbivore  productivity  were  inferred  from  movements  of  the 
grass  zero-isocline  only.  The  differences  in  resilience  were  inferred  from  the  relationships 
between  both  the  G  and  W  zero-isoclines.  Because  the  two  estimates  are  based  on  different 
information,  the  relationship  between  productivity  and  resilience  therefore  need  not  be 
simple  and  direct.  In  this  example,  high  resilience  is  associated  with  intermediate 
productivity. 

These  examples  illuminate  the  concept  of  resilience.  The  original  description  by  Holling 
(1973)  contrasted  resilience  with  the  notion  of  local  stability  near  an  equilibrium.  The  latter 
emphasizes  resistance  to  small  disturbances,  and  rapid  recovery  from  disturbance.  Thus 

Table  2.  Summary  of  properties  of  grazing  systems  in  southern  Africa. 

Herbivore  Variability 

Grazing  pattern  Resilience  productivity  in  time 

Nomadic  and  native  ungulate  grazing  High  Intermediate  High 

Settled  peasant  farmers,  stock  grazing  Intermediate  Low  Intermediate 

Ranching,  fixed  number  of  stock,  and  range 

management  Low  High  Low 
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local  stability  is  associated  with  low  variability.  The  notion  of  resilience  emphasizes  the 
existence  of  more  than  one  domain  of  attraction  (area  bounded  by  the  separatrix  on  the  G 
v.  W  graphs)  and  the  maintenance  of  that  global  structure  through  variability.  In  the 
present  situation,  reduction  of  variability  changes  competitive  balances,  so  that  parameters 
of  the  system  change.  Thus  grass  growth  rate,  rG,  decreases  in  the  settled  peasant  farmer 
systems,  primarily  because  of  reduction  in  the  biomass  of  sensitive  grasses;  the  biomass  of 
ungrazeable  grass,  G0,  decreases  in  constant  stocking  systems  because  of  the  reduction  in 
biomass  of  resistant  grasses. 

Similar  analyses  of  pest  control  and  forest  fires  (Wright  &  Heinselman  1973;  Ludwig, 
Jones  &  Holling  1978;  Holling  1979)  lead  to  the  same  conclusion.  Insecticide  spraying  of 
spruce  budworm  in  eastern  Canada  since  1950  has  eliminated  the  boom-and-bust  epidemic 
pattern,  leaving  more  constant  and  moderate  densities  of  budworm  over  larger  areas.  The 
result  has  been  an  accumulation  of  susceptible  foliage  over  those  same  areas.  Hence  any 
relaxation  of  control  would  result  in  a  more  intensive  and  extensive  outbreak  than  has  ever 
been  experienced.  Effective  control  and  suppression  of  fires  in  the  national  parks  of  the 
western  United  States  reduced  the  incidence  of  fires  but  led  to  an  accumulation  of  fuel  and 
a  more  homogenous  expanse  of  inflammable  canopy.  Again,  any  relaxation  of  control 
would  result  in  costly  consequences.  In  each  case,  reduced  variability  has  led  to  the 
presumed  parameters  changing  to  produce  a  less  resilient  system.  In  the  example  of  the 
fixed  level  of  stocking  on  savanna  (Fig.  18(c)),  the  ratio  of  sensitive  to  resistant  perennial 
grasses  shifts  to  favour  the  former  as  represented  by  decrease  of  the  parameter  G0.  In  the 
budworm  example,  the  spatial  distribution  of  susceptible  foliage  becomes  more 
homogenous.  In  the  forest  fire  example,  the  physical  structure  of  the  forest  shifts  to  allow 
continuity  of  fuel  from  the  ground  to,  and  through,  the  canopy.  Each  system  becomes 
much  more  sensitive  to  errors  in  management  and  to  exogenous  events. 

These  examples  suggest  a  more  precise  definition  of  resilience:  resilience  is  the  ability  to 
adapt  to  change  by  exploiting  instabilities,  rather  than  the  ability  to  absorb  disturbance  by 
returning  to  a  steady  state  after  being  disturbed. 

The  four  examples  shown  in  Fig.  18  represent  different  combinations  of  parameters.  The 
combinations  were  chosen  to  represent  the  particular  parameter  values  that  would  evolve 
under  grazing  regimes  that  differed  in  the  variability  and  intensity  of  grazing.  There  are  two 
questions  concerning  resilience.  The  first  is  how  to  compare  the  resilience  of  the  different 
grazing  systems,  given  their  parameters.  The  second  is  how  to  predict  the  direction  and 
magnitude  of  parameter  change,  given  the  set  of  equations. 

The  earlier  description  of  differences  in  resilience  between  three  southern  African 
examples  related  to  the  size  of  the  upper  area  of  stability  on  the  G  v.  W  graphs.  But  the 
descriptions  were  meaningful  only  in  terms  of  kind  and  direction  of  disturbance.  The 
important  disturbances  were  increase  in  H,  decrease  in  G,  and  leftward  movement  of  the 
grass  zero-isocline  induced  by  drought.  Because  of  these  directions  of  disturbance,  a 
convenient  common  reference  point  would  be  the  lowest  point  in  the  set  of  upper  equilibria 
(B,  Fig.  18).  One  possible  set  of  measures  of  resilience  would  be  the  distances  in  H,  G,  W 
space  from  the  reference  point  to  the  separatrix.  Since  the  separatrix  is  complicated  to  find, 
surrogates  would  be  distances  to  the  saddle  point,  D,  through  which  the  separatrix  passes. 
Such  resilience  measures  would  decrease  from  (a)  to  (b)  to  (c)  in  Fig.  18. 

The  second  aspect  of  resilience  concerns  the  effects  on  parameters  of  controlling  the 
variability  and  magnitude  of  variables  (Peterman,  Clark  &  Holling  1979).  Many 
parameters  of  ecological  models  are,  in  reality,  variables.  They  are  chosen  to  be 
represented  as  constants  (or  perhaps  as  stochastic  variables)  because  they  are  likely  to 
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change  too  slowly  to  be  of  interest  (e.g.  geomorphological  changes),  because  of  convenience 
and  the  need  for  simplification  (e.g.  herbivores  in  the  model  presented  here)  or  because 
information  or  understanding  of  their  behaviour  is  lacking.  And  yet  the  capability  to  adapt 
to  disturbance  is  very  much  tied  to  the  stability  of  key  parameters  and  their  tendencies  to 
evolve  under  changed  circumstances. 

The  general  solution  is  not  simply  to  expand  the  model  to  represent,  where  possible, 
parameters  as  variables.  Such  an  expansion  itself  introduces  additional  parameters  that  are 
subject  to  the  same  uncertainty  concerning  their  stability.  If  anything,  the  problem  is 
magnified  because  the  number  of  parameters  increases  rapidly  as  a  model  is  expanded.  An 
explicit  effort  to  structure  models  as  hierarchies  can  help  to  formalize  the  problem.  That  is, 
the  problem  can  be  structured  as  a  hierarchy  and  defined  by  focusing  on  three  or  four 
levels  of  that  hierarchy  that  differ  in  their  speed  or  geographical  scale  or  both.  In  the 
problem  described  here  the  speed  of  variables  decreases  from  sensitive  grasses  to  resistant 
grasses,  to  woody  shrubs  to  herbivores.  Events  at  levels  below  this  hierarchy  will  operate 
so  rapidly  it  might  be  assumed  that  an  averaging  effect  allows  their  representation  as 
parameters.  Events  at  still  higher  levels  might  be  presumed  to  operate  so  slowly  that  they 
can  be  represented  as  constants.  Although  this  formalizes  the  problem,  it  does  not  resolve 
it  since  very  fast  variables  can  be  eliminated  by  man's  activity  (e.g.  soil  microfauna 
important  for  maintaining  soil  fertility)  and  normally  slow  variables  can  be  speeded  up  (e.g. 
hydrodynamic  processes  through  irrigation). 

At  best  we  can  suggest  only  a  modest  step.  Whatever  the  choice  of  variables,  explicit 
efforts  can  be  made  to  infer  the  consequences  on  what  was  left  out  of  the  model.  Often 
there  is  some  information  that  can  suggest  the  direction,  if  not  the  magnitude,  of  change  of 
key  parameters.  That,  in  essence,  was  done  in  constructing  the  four  examples  of  Fig.  18. 
Knowledge  of  the  change  in  numbers  of  grazers  and  of  the  competition  between  the  two 
kinds  of  grasses  led  to  inferences  about  the  way  parameters  would  change  under  different 
grazing  regimes,  and  that  generated  predictions  of  changes  in  productivity  and  resilience 
measures.  It  is  then  possible  to  design  the  key  field  experiments  that  can  best  test  these 
qualitative  predictions.  Where  these  can  be  designed  as  part  of  management,  additional 
information  on  economic  and  social  feasibility  can  be  produced. 

In  southern  Africa  one  such  key  management  experiment  would  be  a  pulsed  grazing 
scheme  within  subsections  of  a  typical  ranch.  Periods  of  intense  over-grazing  would 
alternate  with  periods  of  recovery  of  productive  grasses,  the  latter  accelerated  by  various 
experimental  techniques.  The  timing  of  the  grazing  pulses  would  rotate  among  the 
subsections  so  that  temporal  variability  for  the  ranch  as  a  whole  would  be  reduced  by 
transforming  it  to  spatial  variability.  Productivity  (biological  and  economic)  and  resilience 
changes  would  be  measured  until  they  had  stabilized.  That  would  be  followed  by  controlled 
application  of  different  kinds  and  intensities  of  stress  as  a  key  test  of  the  qualitative 
predictions  discussed  earlier.  Similar  tests,  applied  to  areas  managed  with  moderate  fixed 
stocking  strategies  and  to  those  subjected  to  intense  grazing  by  stock  kept  by  settled 
peasant  farmers,  would  provide  the  necessary  controls. 
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VARIABLE  DICTIONARY 


AHEN  -  variable;  total  number  of  adult  hens 

AL  -  constant;  the  ratio  of  water  uptake  by  shrubs  to   water 
uptake  by  grass 

ANEST  -  variable;  proportion  of  adult  females  which  nest   as 
a  function  of  grass  biomass 

BE  -  constant;  ratio  of  evaporation  loss  rate  to   water   up- 
take rate  by  grass 

CI  -  constant;  rate  at   which   infiltration   rate   increases 
with  increasing  grass  biomass  (kg/ha) 

C2  -  constant;  rate  of  consumption  of  grass  biomass  per  unit 
of  herbivore  biomass  (kg/kg) 

CHICKS  -  variable;  total  number  of  chicks  born 

CPHEN  -  variable;  number  of  chicks  produced  per  nesting   hen 
as  a  function  of  grass  biomass 

D4  -  constant;  grass  biomass  at  which  shrub  growth  rate  from 
subsoil  water  is  50%  of  maximum  (kg/ha) 

DE  -  constant;  ratio  of  evaporation  loss  rate  to   water   up- 
take rate  by  shrubs 

E  -  constant;  evaporation  loss  rate  (mm/year) 

G  -  variable;  grass  biomass  (kg/ha) 

GO  -  constant;  grass  biomass  below  which  herbivore   consump- 
tion does  not  occur  (kg/ha) 

GA  -  constant;  the  ratio  of  water  uptake  by  grass   to   water 
uptake  by  shrubs 

H  -  constant;  herbivore  biomass  (kg/ha) 

10  -   constant;   minimum   infiltration   rate   at   low   grass 
biomass  levels 

K  -  constant;  proportional  reduction  in   shrub   growth   rate 
from  RS  to  S  when  grass  biomass  is  high 

LG  -  constant;  the  proportional  rate  of  loss  of   grass   from 
respiration  and  death  (/year) 

LW  -  constant;  the  proportional  rate  of  loss  of   grass   from 
respiration  and  death  i./yesr) 


MOD  -  variable;  hen  survival  rate  from  shrub  cover  effects 

NS  -  variable;  variable  which  points  to  proper  set  of  X   and 
Y  points  to  use  in  SLP  interpolation 

P  -  constant;  inefficiency  of  subsoil  water  use  by  shrubs 
(kg/ha) 

POOLER  -  variable;  percent  shrub  cover  as  a  function  of 
shrub  biomass 

PFEM  -  constant;  proportion  of  population  which  is  female 

Q  -  the  rate  of  grass  consumption  by  herbivores  (kg/ha/year) 

RQ  -  constant;  grass  growth  rate  due  to  topsail  water 
(kg/ha/year) 

RS  -  constant;  the  increase  in  biomass  of  shrubs  due  to  sub- 
soil water  when  there  is  no  grass  biomass 
(kg/ha/year) 

RN  -  constant;  shrub  growth  rate  from  topsoil  water 
C kg/ha/year ) 

S  -  variable;  actual  shrub  growth  rate  from  subsoil  water 
(kg/ha/year) 

SHEN  -  variable;  hen  survival  rate  as  a  function  of  grass 
bi  omass 

SMALE  -  variable;  male  survival  rate  from  shrub  cover  ef- 
fects 

TG  -  constant;  grass  water  uptake  rate  (mm/year) 

THN  -  variable;  total  hen  population  (yearlings  plus  adults) 

TMALE  -  variable;  total  male  population  (yearlings  plus 
adults) 

TT  -  variable;  total  number  of  hens  nesting 

TW  -  constant;  water  uptake  by  shrubs  (mm/year) 

UG  -  variable;  relative  utilization  of  topsoil  water  by 
grass  per  unit  of  grass  biomass  (mm/kg) 

UW  -  variable;  relative  utilization  of  topsoil  water  by 
shrubs  per  unit  of  shrub  biomass  (mm/kg) 

W  -  variable;  shrub  biomass  (kg/ha) 

NHSUR  -  constant;  hen  winter  survival 
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WMSUR  -  constant;  male  winter  survival 

XI  -  variable;  infiltration  rate  (proportion  of  precipita- 
tion which  enters  soil  when  grass  biomass  is  at 
maximum) 


XN  -  variable;  number  of  interpolation  points  for  SLP  to   be 
evaluated 

XX(5,6)  -  constant;  the  array  of  X  points  for  the  SLP   func- 
tions 

YHEN  -  variable;  total  number  of  yearling  hens 

YNEST  -  variable;  proportion  of  yearling  females  which   nest 
as  a  function  of  grass  biomass 

YY(5,6)  -  constant;  the  array    of  Y  points  for  the  SLP   func- 
tions 

ZR  -  variable;  normal  random  variate  used  to  vary       infiltra- 
tion rate 
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Z  Array  Storage  Variables 


1  -  XI 

2  -  UG 

3  -  Q 

4  -  G 

5  -  UW 

6  -  S 

7  -   \4 

8  -  H 

9  -  AHEN 

10  -  AMALE 

11  -  TMALE 

12  -  THN 


SLP  Functions 
(refer  to  first  AEA  Inc  Saval  Report  for  Figures) 

1  -  Adult  nesting  success  Mrs  grass  biomass  (Fig.  £*€) 

2  -   yearling  nesting  success  vrs  grass  biomass  (Fig.  6.6) 

3  -  Chicks  born  per  hen  vrs  grass  biomass  (Fig.  6.7) 

4  -  Hen  survival  vrs  grass  biomass  (Fig.  6.8) 

5  -  Hen  survival  vrs  percent  cover  (Fig.  6.9) 

6  -  Male  survival  vrs  percent  cover  (Fig.  6.10) 
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PROGRAM  LISTING 


I  DIM  2(20, 100), 2M!20) 

10  FOR  1  =  1  TO  20-.ZMII)  ■  0:  NEXT  :NV  =  20 

II  IF  PEEK  1104)  <  >  64  THEN  POKE  104,64:  POKE  103,1:  POKE  16384,0:  PRINT  CHR*  (4);"RUN  SIHCON" 

14  INPUT  'SIM  FROM  'iA$:  IF  LEFT!  (A*,21  «  "BR"  THEN  10100 

15  ZS  =  VAL  (  LEFT*  (A$,2)):NT  =  VAL  (  RI6HTJ  (A$,31) 

16  IF  NT  (  1  THEN  NT  =  VAL  (  RIGHT*  (At, 2)) 
18  FOR  TIME  =  ZS  TO  NT 

100  IF  TI  >  0  THEN  1000 

1 10  H  =  0:6  =  200:LS  =  1:80  =  500 

120  10  =  O.l'.TG  =  0.2i:TK  «  0.122:RG  =  4000 

130  CI  =  600:C2  *  3.15:C3  =  0.0:E  =  0 

131  AL  ~  TH  /  TG:BE  =  E  /  TG 

135  DATA   0,0, .35, 150, .8, 500, 0,0, 0,0 

136  DATA   0,0, .2, 150, .65, 500, 0,0, 0,0 
138  DATA   0,0,1,100,5,500,0,0,0,0 
140  DATA   .6,0, .6, 100, .9, 500, 0,0, 0,0 
142  DATA   .25,0, .25, 5, 1,25, .5, 50, .5, 100 
144  DATA   .35,0, .35, 5, .65, 15, .65, 30, .5, 50 

146  FOR  IN  =  1  TO  6 

147  FOR  I  =  1  TO  5:  READ  YY(I,IN!,XX(I,IN):  NEXT  I 

149  NEXT  IN 

150  WHSUR  -   .fiNKSUR  =  .8 

152  AKALE  =  100-.AHEN  =  100:YHALE  =  100:YHEN  =  100 

153  PFEM  =  .6 
200  Kb  0 

210  RH  =  3000:RS  =  4000:LK  =  1:K  =  1:D4  =  500:P  =  .1 

99e  "mLg  =1oD00=  BE  7  ^  Equation   #s   in  Walker  et.    al    (1981) 

999  PRINT  'ENTER  COMMANDS, THEN  CONT":  END 

1000  XI  =  (G  +  CI  ♦  10)  /  (G  +  CD  -------------  (1) 

1005  BOSUB  12000: XI  ■  XI  *  EXP  (.25  *  ZR) 

1010  US  =  1  /  fB  +  ftL  *  H  +  BE)  -  -  -  -  -  -  -  -  -  -  -  -  -  -  (5) 

1020  0  -   (C2  *  (1  -  GO  /  6!)  /  (6  +  C3>:  IF  Q  <  0  THEN  Q  =  0   -----  (3) 

1030  6  =  6  +  RS*XI*US*6-L6*G-H*Q*G   --------  (4) 

1040  IF  6  <  IE  -  5  THEN  G  =  IE  -  5 

2000  UN  =  1  /  (M  +  GA  *  G  +  DE)  -  -  -  -  -  -  -  -  -  -  -  -  -  -  (9) 

2010  S  =  RS  *  (1  -  6  /  (K  *  (G  +  D4))l  -  -  -  -  -  -  -  -  -  -  -  -  before  (9)  [unnumbered] 

2020  K  =  N  +  N  *  RH  *  II  *  UK  +  (S  *  H  /  (H  +  PI)  -  LK  *  H  -  -  -  -  -  -  (10) 

2030  Z(5,TI)  =  UK:Z(6,TI)  =  S:Z(7,TI)  ■  H 

3000  REM    NESTING  SUCCESS  ADULTS 

3010  XJ  =  6 

3015  NSLP  =  1 

3020  XN  s  3 

3025  BOSUB  13000 

3030  ANEST  -  SLP 

3035  RED   NESTING  SUCCESS  YEARLINGS 

3040  XJ  s  G 

3045  NSLP  =  2 

3050  IN  =  3 

3055  60SUB  13000 

3057  YNEST  =  SLP 

3060  REM    CHICKS  PER  HEN 

3062  XJ  =  G:NSLP  =  3:XN  =  3  "  ^  *  9  * 

3063  GOSUB  13000 
3065  CPHEN  =  SLP 

3070  REM   I  OF  CHICKS 

3078  TT  b  AN  *  AHEN  +  YN  «  YHEN:CH  =  CPHEN  *  TT  *  (1  -  TT  /  5000) 

30BO  REM    HEN  SURVIVAL 


3015  W  *  6 
30V0  NSLP  *  4 
3095  KN  *  I 
3100  BOSUB  13000 
3105  SHEN  *  SLP 

3110  REM    MODIFY  MITH  COVER 

3111  PCOVER  =  (N  /  4000)  «  100 

3112  IF  PCOVER  >  100  THEN  PCOVER  =  100 
3115  11   =  PCOVER 

3120  NSLP  =  5 

3125  m  -  4 

3130  60SUB  13000 

3135  MOD  a  SLP 

3140  REM   MALE  SURVIVAL 

3145  U  -  PCOVER 

3150  IN  ~  5 

3155  NSLP  -   6 

3160  BOSUB  13000 

3165  SHALE  ■  SLP 

3170  AHEN  «  (AHEN  +  YHEN)  *  SHEN  *  MOD 

3175  AMALE  -   (AMALE  +  YMALE)  *  SMALE 

3180  YHEN  -  CHICKS  *  PFEM 

3185  YMALE  «  CHICKS  *   (1  -  PFEM) 

3190  REM   ADD  SN  WINTER  MORTALITY 

3195  AHEN  *  AHEN  *  HHSUR 

3200  AMALE  *  AMALE  *  HMSUR 

3205  YHEN  -  YHEN  *  HHSOR 

3210  YMALE  =  YMALE  *  MSliR 

3215  THN  =  AHEN  +  YHENsTMALE  =  AMALE  *  YMALE 

3220  PRINT 

3300  PRINT  ■««"»■(  AN,"  YN='YN 

3310  PRINT  aCP=";CP,BSH=BsSH 

3315  PRINT  "PC=';PC,"MQD=B?HOD 

3320  PRINT  "SM=*;SM 

5000  Z(iJI)  =  H:Z(2,TI)  *  UG:Z(3,TU  *  llZH.TH  2  B 

5005  Z(8,TI>  =  H 

5006  Z19,TI>  =  AHENsZQOJD  -  AHALEiZMJI)  ■  THALE:Z(12,TI)  »  THN 
10000  FOR  0  =  1  TO  NV 
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I.   INTRODUCTION 

This  model  is  a  simplified  version  of  the  hydrology  and 
vegetation  submodels  described  in  Sonntag  et  al.  (1982)  and  the 
first  section  of  this  report. 

The  model  considers  only  one  range  site  on  one  acre  of  a 
pasture.   However,  the  time  step  is  flexible;  it  is  determined 
by  setting  the  variable  ND  (e.g.,  7  for  a  weekly  time  step).   The 
basic  structure  of  the  model  is  as  follows. 

Changes  in  soil  moisture  in  two  soil  layers  (0  -  20", 
20  -  40")  are  computed  based  on  precipitation,  infiltration  and 
evapotranspiration.   The  actual  above  ground  biomass  and  potential 
above  ground  biomass  stored  in  roots  are  modelled  for  both  grasses 
and  shrubs.   Grasses  draw  water  only  from  the  top  soil  layer;  shrubs 
draw  water  from  whichever  soil  layer  has  more  water.   Growth  rates 
of  vegetation  depend  on  soil  moisture  levels,  a  phenology/temperature 
factor,  the  vegetation  type  (grasses  or  shrubs) ,  a  competition 
factor,  and  translocation  between  above  and  below  ground  biomass. 

Cattle  grazing  can  be  set  to  occur  in  any  period  of  one 
week  or  longer;  the  user  sets  the  number  of  acres  per  cow,  and  the 
estimated  weekly  biomass  of  forage  consumed  per  cow. 
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VARIABLE    DICTIONARY 
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AB(2)  actual  biomass  of  grasses  &  shrubs  (green  and  non-green) 

AC  #  acres /cow 

ANS  temporary  variable  (for  growth  rate) 

AW  available  water  capacity  ("/") 

B2  green  biomass  "mortality"  when   RG  =  0.0 

BC(2)  base  cover  (%  cover  with  zero  biomass  of  grass  &  shrubs 

■  respectively) 

BG  fraction  bare  ground  (e.g.  0.2) 

CC  cattle  consumption  (lbs/week/cow) 

DD(2)  death  rate/week  of  grasses  and  shrubs 

DP (2)  depths  of  2  soil  layers 

DS  week  of  seed  set 

I       EF  growth  reduction  factor  for  competition  for  water 

ET  actual  evapotranspiration/week  (") 

FD  week  of  flowering 

GB(2)  above  ground  green  biomass  (lb/acre)  of  grasses  (1)  & 
shrubs 

IW  week  #  in  year 

KB (2)  half  saturation  constant  for  %  cover  relationship 

KQ  half  saturation  constant  for  soil  moist. ->  plant  growth 
equation 
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MB (2)  maximum  biomass  (100%  cover) 

MR (2)  reduction  in  growth  rate  due  to  moisture 

MX (2)  maximum  growth/week  of  grasses  &  shrubs  respectively 

ND  #  days/timestep 

[J       NE  necessary  evapotranspiration  (in/week  per  lb/a  of  biomass) 

NQ  equation  #  (pointer) 


0 

I 
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old  green  biomass 

%  cover  of  grasses;  shrubs  (grasses  &  shrubs)  respectively 

phenology /temperature  factor 

total  potential  water  transpiration 

random  number 

radiation  (langleys) 

real  growth  rate  of  plant 

runoff  (inches) 

r  parameter  for  soil  moisture  -*■  plant  growth  equation 

Pr  (no  rain) 

translocation  from  root  to  shoot  (fraction  of  root  storage/ 
week) 

root  biomass  (lb/a)  (potential  above  ground  green  biomass) 

actual  retention  parameter 

slope  of  mortality-growth  line 

maximum  retention  parameter 

%  soil  moisture  by  layer 

translocation  from  shoot  to  root  (fraction  of  above  ground 
green  biomass/week 

SS(1)    saturation  parameter  for  soil  moisture  ■+■  plant  growth 
equation 

TB  total  above  ground  biomass 

TE  temperature  (°F) 

TI  time  (#  weeks) 

TP  cumulative  precipitation 

WD (2)  water  density  in  wettest  of  two  layers  ("/") 

WI  total  infiltration 

WN  potential  evapotranspiration  p  r\  p.  C)  £ 


0G(2) 

PC(3) 

PF 

PW 

R 

RD 

RG 

RO 

RQ 

RR 

RS(2) 

RT(2) 

s 

S2 

SM 

SO  (2) 

SR(2) 
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WT(2)    water  transpired  from  each  soil  layer  (") 
WW (2)    amount  of  water  in  each  soil  layer  (") 
XX      temporary  variable  for  soil  moisture 
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105  - 

199 

210  - 

240 

270 

300 

320 

330 

350  - 

430 

470  - 

530 

550  - 

590 

600  - 

710 

1000 

-  1020 

1040 

-  1210 

1260 

-  1300 

1310 

-  1320 

1330 

-  1360 

1370 

-  1380 

1384 

1390 

-  1430 

1440 

-  1560 

Initialization. 

Compute  temperature  and  radiation  from  sin  curves 

Compute  potential  evapotranspiration 

Allocate  potential  evapotranspiration  to  each  soil 
layer,  accounting  for  %  runoff  and  %  green 

Choose  wettest  layer  for  shrubs 

Compute  water  transpired  from  soil  layer  if  layer 
filled  to  AWC  (water  not  limiting) 

Take  water  out  of  soil 

Compute  rainfall  (sin  curve  for  probability  of  no 
rain;  random  number  for  actual  rain) 

Compute  curve  numbers,  retention  parameter 

Fill  up  water  layers  with  infiltration 

Compute  phenology/temperature  factor 

Compute  soil  moisture  constraint  (using  logistic 
equation) 

Compute  evapotranspiration  limitation  on  growth 
Compute  real  growth  rate  of  plant 

Death  rate  (proportional  to  phenology  factor  if 
early  spring  or  late  fall) 

Growth  and  death  of  green  biomass 
Cow  grazing 

Translocation  (root  storage  after  flowering  removes 
green  biomass  but  not  total  biomass) 

Compute  %  cover  based  on  biomass 
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1  DIM  2(30,52),ZM<30) 

2  OIK  S0(2),HR(2),Hn2),SS(l),RO(l),Ke(l),6E(2),RT(2),O6(2),R6(2),KB(2),DP!2),MW(2),«U(2),KT(2),RSl2),SR(2),BC12),DI)t2),S2(2) 

10  FOR  I  =  1  TO  SOlZHtl)  ■  0:  NEXT  :NV  =  30 

11  IF  PEEK  (104)  <  >  64  THEN  POKE  104,64:  POKE  103,1:  POKE  16384,0:  PRINT  CHR$  (4)|*MM  SIMCON" 

12  PRINT  'ENTER  COMMANDS, THEN  CONT":  END 

14  INPUT  'SIM  FROM  ';A$:  IF  LEFTS  (A$,2)  =  "BR"  THEN  10100 

15  ZS  =  VAL  (  LEFT*  (A$,2)):NT  =  VAL  (  RIGHT*  (A$,3>) 

16  IF  NT  <  1  THEN  NT  =  VAL  (  RI6HT*  (A$,2)) 
18  FOR  TIME  =  ZS  TO  NT 

100  IF  TI  >  0  THEN  GOTO  200 
105  REM   INITIAL  CONDITIONS 
107  AW  =  0.14 


110  DF(1) 

=  20:DP 

(2 

=  2 

5 

115  m\\ 

=  DP(l) 

1 

AK  * 

0.75 

117  HH12I 

=  DP  (2) 

1 

AH  * 

0.25 

120  R  * 

RND  (  - 

i) 

125  TP  = 

0 

126  CC  = 

200 

127  AC  = 

20 

128  PF  = 

1 

i 

0 


130  SS(1)  =  1 
135  RGil!  =  0.8 
140  KQ(1)  =  15. 

145  KM)  =  1.2 

146  6S  =  23 

147  SF  =  26 

150  HX(2)  =  0.6 

151  RS(1)  =  -  0.1:RS(2)  =  -  0.1 

152  SR(1)  =  0.05:SR(2)  =  0.05 
155  IN  ■  0 

160  NE  =  0.01 

16!  BC(l)  "  0.1:BC(2)  =  0.1 

162  KB(1)  =  1000:KB(2)  =  1000. 

165  BS  =  0.2 

165  B2  =  0.10 

170  S2U1  =  0.07:52(2)  =  0.07 

172  PI  =  3.1*159:ND  =  7 

175  FD  =  26 

177  B3  =  0.1  ♦  FD  /  (FD  -  8) 

178  S3  =  B3  /  FD 

190  Bid)  =  50.:AB(1)  =  50 
195  BB(2)  =  50.:AB(1)  =  50 

197  RT(1)  =  200. 

198  RT(2)  =  200. 

199  PRINT  'ANY  CHANGES?  THEN  CONT":  END 

200  IK  =  IH  +  1 

205  IF  IN  *  (ND  /  7)  >  52  THEN  IH  =  0 

207  IF  IB  «  0  THEN  TP  =  0:CR  =  0 

210    REM       TEMPERATURE*«tmm 

220  TE  =  10  +  60  *    SIN  ((IN  *  (ND  /  7)  /  52)  *  PI) 

240  REM   RADIATION  *********** 

250  RD  =  200  +  600  *  SIN  ((IN  *  (ND  /  7)  /  52)  *  PI)  n  A  O  Q  O 

260  REM   POTENTIAL  ET  ********  l ;  v-  U 

270  HU  =  ND  *  (0.014  *  TE  -  0.37)  *  RD  /  580. 

280  IF  U  <  0  THEN  HU  =  0 

285  PW  =  HU 

290  FOR  I  =  1  TO  2 

300  ITU)  =  HU  *  0.4  *  PF  *  PCU) 


510  KM  KM1)    /   L>     ■  > 

320    if  i  --  :  an;-  i    i:»  /  Dn?n  >   «  kd  then  up  *  W(2)  /  DP(2) 
330  men  *  Min  *  k  »  hi  /  ak 

340     N[)l 

3'jO  KU  «  11(11   +  tftZl 

360  El  =  0 

370     IF  mil)  <  HKUI  THEN    GOTO  390 

375  ET  =  ET  ♦  WW  ( 1 ) 

380  KKtl)  =  0 

3E5    6CTO  100 

390  HM(l)  --  Mf(l)  -  HTClllET  «  ET  *  mill 

400  HX  '-  2:   IF  HH(l)   >  HK(2)  THEN  HX  =  1 

410  IF  NT  12)  <  HH(NX)  THEN  60T0  430 

415  ET  -  ET  ♦  KH(NX):KN(NX)  =  0 

425  GOTO  460 

430  mm)  *  mm  -  htuuet  =  et  +  mm 

460  PRINT  "     «*«*»*■ 

46!  PRINT  'KEEK  \IN  *  (NE  /  7);"  TE="|TEi"  RD=';RD 

462  PRINT  ■  AFTER  EVAPOTRANSPIRATION:" 

463  PRINT  -WK(1)=-;WW(1>;-  Wi(2)«MW(2):  PRINT  •HT<l)=8;KTU)s'  HT (2)«*;HT  12) s  PRINT  'PKVjPH 
470  R  =  RND  (1) 

4e0  RF;  =  SIN  ((IK  t  (ND  /  7)  /  52)  *  PI)  *  ©.72 

490  IF  R  <  RR  THEN  RN  =  O:  SOTO  550 

520  RN  *  (  RND  (1)  *  0.5  +  0.5)  *  (ND  /  7) 

530  TF  =  TP  +  RN 

550  CN  •  87  -  0.18  *  PC (3) 

540  Sn  «  (1000.  /  CM  -  10. 

576  S  «  St!  *  (AK  -  (HH(1)  ♦  IW(2»  /  CDPtlJ  ♦  DP(2>»)  /  AH 

575  RO  =  0 

580  IF  RN  >  0.2  *  S  THEN  RO  =  ((RN  -  0.2  «  S)  *  2)  /  (RN  ♦  0.8  *  S) 

590  IF  RN  <  (0.2  *  SI  THEN  RO  »  0. 

595  CR  =  CR  ♦  RO 

600  REN   INFILTRATION  *«**»»* 

610  MI  -  RN  -  RO 

630  IF  HI  <  0.  THEN  HI  -  0. 

635  NF  =  HI 

640  FOR  M  1  TO  2 

650  IF  (HI  *  HHiD)  <  (DF(I)  *  AK)  GOTO  690 

660  WI  =  Ml  -  (DP(I)  *  AH  -  Witt)) 

670  KH(I)  =  DPII)  *  AW 

680  60T0  710 

690  HHil)  i  MH(I)  +  HI 

700  6QT0  720 

710  NEXT  1 

720  PRINT  '  AFTER  INFILTRATION;':  PRINT  •NK(1)=';HH(1)?'  «K(2)sB?NK(2k  PRINT  'RN=S;RN 

740  TB  =  68(1)  +  61(2) 

1000  REM   PHENOLE;  I  TEHP.  «tt 

1020  PF  =  SIN  (IN  «  (ND  /  7)  /  52  «  PI) 

1040  REH   SOIL  BOISTURE  CONSTRAINT  «* 

1050  TB  =  BB(1)  ♦  86(2) 

1060  FOR  1  =  1  TO  2 

10B0  SO(I)  =  ((KH(I)  /  DP(D)  /  AH)  *  100. 

1090  NO.  =  1 

1100  U  -  S0(1) 

1110  IF  1  =  2  AND  80(2)  >  S0(1)  THEN  H  =  50(2) 

1180  60SUB  1690 

1200  HRU)  =  AN5 

1210  IF  Mtd)  )  1  THEN  HR(I)  =  1 

1240  REH   EVAP0TRANSPIRAT10N  »«i  ^ 

1260  EF  '-  (ET  /  TB)  /  (NE  »  (ND  /  7))  ,  J  ]   (J  'J 

1280    IF  EF  >  1.0  THEN  EF  -  1.0 

1300    IF  EF  (  0.0  THEN  EF  =  0.0 

1310  R6II)  '-  fix  (I)  «  PF  ♦  HR(I)  *  & 

1315     IF  1H  «   (ND  /  7)   <  8  THEN  RGII)  =  0. 

1320     IF   IN  ♦   (NO  /  7)    >  40  THEN  RGd)  =  0 

1330    REH         H0RTAL1TV  §♦♦»♦♦♦ 
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D 
D 
D 

G 

D 
D 
0 
0 
D 
D 
I 
D 
D 
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srii: 


=  RS(I)  *  RT(I)lGBU) 
*  6B(I):6B(I)  =  GB(i; 


=  6B(I> 

+  DR 


+  DR:AB(1)  =  AB(I)  +  DR 


u&  duu)  -  y 

1340  IF  IK  »  (ND  /  7)  <  B  OR  1M  >  FD  THEN  DD(I)  *  B2  -  S2(I)  *  PF 

12.60  IF  DD(I)  <  0.  THEN  DDf I>  =  0. 

1370  OB(I)  «  6B(I) 

1380  BBU)  --   BB(I)  *  (ND  /  7)  *  (RE(I)  +  1.  -  DD(I)):AB(I)  =  AE(I)  *  (ND  /  7)  *  (RB(I)  +  1  -  DDU!) 

13B2  REM  COW  EATS*"* 

1384  IF  IK  *  (ND  /  7)  >  GS  AND  IK  *  (ND  /  7)  <  GF  AND  I  =  1  THEN  GB(1)  «  SB ( 1 )  -  CC  *  (ND  /  7)  /  AC:AB(1)  =  AB(1)  -  CC  *  (ND  /  7)  / 

AC 

1390  REM 

1392  RS(I)  =  B3  -  S3 

1395  DR  =  0 

1400  IF  IK 

1410  IF  IK 

1430  RT(I)  =  RT(I) 

1440  NEXT 

1450  REM   I  COVER  ***** 
1520  FOR  I  =  1  TO  2 

1530  PC(I)  =  (ABU)  *  (t.  -  8C(I)  -  BG))  /  (KB(I)  +  AB(D)  +  BCCI1 
1540  IF  PC(I)  >  1  THEN  PC(I)  =  1.0 
1550  NEXT  I 

1555  PC (3)  =  PCU1  +  PCI21 

1560  IF  PC(3)  >  1  THEN  PC(1)  =  PC(1)  /  PC(3):PC(2) 
1565  2(21, TI)  =  ABU) 
AB(2) 
=  ET 


TRANSLOCATION  «*** 

*  IK  *  (ND  /  71 

(ND  /  7)  <  FD  AND  IK  *  (ND  /  7)  >  B  THEN  DR 
(ND  /  7)  >  FD  AND  IK  <  40  THEN  DR  i 
DR 


PC(2)  /  PC(3):PC(3)  =  PC(1)  +  PC(2) 


1567  Z(22,TI) 
1569  2(23, TI) 
1571  Z(24,TI>  =  RD 
1573  2(25, TI)  =  CR 
15B0  2(1, TI)  =  61(11 
1590  2(2,TI)  =  68(2) 
160O  2(3, TI)  =  50(1) 
1610  2(4, TI)  =  50(2) 
1620  2(5, TI)  =  RN 
1630  2(6, TI)  =  PC(1) 
1640  2(7,TI)  =  PC(2! 
1650  2(8, TI!  =  PF 
1660  Z(9,TI)  =  RT11J 

1670  ZtlOJI)  =  RT(2) 

1671  2(11, TI)  =  RS(1) 

1672  2(12,TI)  =  RS(2) 

1673  2(13, TI)  =  TP 

1674  2(14, TI)  =  R6111 

1675  2U5.TI)  »  RG(2) 

1676  2(16, TI!  =  RO 

1677  ZU7JI!  =  NF 
167B  2(18, TI)  =  KT(1) 

1679  2(20, TI)  =  MRU) 

1680  GOTO  2000 

16B9  REM  LOG  IST1CEQN  ***** 

1690  AN5  -   SS(NQ)  /  (1.0  ♦  EXP  (RQ(NQ)  *  (KQ(NQ)  -  XX))) 

1700  RETURN 

2000  REM   ***END*»* 

1O000  FOR  J  =  1  TO  NV 


"i  1 


n  i 


cd 


*,B* 


•■; 
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